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ABSTRACT 


The  purpose  of  this  investigation  is  to  study  the  noise  transport 
and  reduction  phenomena  in  multivelocity  electron  beams  with  significant 
potential  variations  in  the  transverse  direction.  Correlations  between 
the  noise  reduction  mechanism  and  the  multi -dimensional  nature  of  the 
electron  beams  are  sought.  The  Investigation  was  carried  out  by  means 
of  a  high-speed  digital  computer  employing  two  different  numerical 
methods.  The  mathematical  model  assumes  an  infinite  confining  magnetic 
field  so  that  transverse  motions  are  not  allowed.  In  addition  a  one- 
dlmensional  half -Maxwellian  velocity  distribution  function  is  considered 
for  the  electrons  at  the  potential  minimum. 

Details  of  the  derivation  of  a  small-signal  linearized  Boltzmann 
equation  in  terms  of  a  statistical  electron  density  function  and  the 
electric  field  in  a  two-dimensional  space-charge-limited  diode  are 
presented,  together  with  some  discussion  on  the  singularities  in  the 
boundary  conditions  and  the  formulation  of  the  a-c  electric  fields.  To 
facilitate  a  solution  of  the  transport  equation,  velocity  classes  are 
introduced  as  an  approximation  so  that  digital  techniques  become  appli¬ 
cable.  Haus's  noiue  theory  for  small-signal  longitudinal  beam  devices 
is  employed  to  evalviate  the  noise  transport  characteristics  of  the 
two-dimensional  electron  beam. 

In  parallel  with  the  density  function  analysis,  the  transport 
phenomena  in  a  two-dimensional  diode  are  treated  with  the  Monte  Carlo 
technique  in  which  random  numbers  are  generated  to  simulate  the  stochas¬ 
tic  emission  process  and  each  electron  emitted  from  the  thermionic  cathode 
is  traced  in  phase  space  by  nvimerical  integration  methods.  The  methods 
employed  to  evaluate  the  d-c  and  a-c  electric  fields  are  discussed, 
especially  the  implications  of  the  widely  used  open-circuit  assumption. 

A  brief  description  of  the  generalized  harmonic  analysis  is  included  and 
Haus's  noise  parameters  are  estimated  statistically  by  means  of 
sampling  techniques.  Because  of  the  thin  beam  assumption,  the  axial  a-c 
electric  field  is  assumed  to  be  uniform  in  the  transverse  direction. 

Only  the  d-c  multi-dimensional  effects  are  actually  included  in  the 
Monte  Carlo  calculations. 

Results  from  the  density  function  calculations  are  presented. 

No  reduction  from  the  thermionic  emission  shot  noise  is  observed,  nor 
is  there  any  evidence  of  the  multi -dimensional  noise  reduction  effects 
expected.  Re-examination  of  the  coupling  mechanism  employed  reveals  that 
the  z-directed  a-c  electric  field  decays  too  rapidly  in  the  transverse 
direction  to  bring  forth  any  signlficeut  two-dimensional  effects.  Because 
of  the  velocity  dependent  characteristics  of  the  coupling  fields,  the 
heavy  emphasis  placed  on  the  fluctuations  associated  with  the  slow  elec¬ 
trons  results  in  a  net  reduction  of  the  kinetic  power  carried  by  the 
electron  beam,  causing  deterioration  in  noise  performance. 
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Substantial  noise  reduction  at  low  frequencies  for  the  two- 
dimensional  diode  is  found  by  means  of  the  Monte  Carlo  calculations. 

A  minimum  noise  figure  of  less  than  2  db  is  predicted  for  0-type  devices. 
The  reduction  in  beam  noise  is  believed  to  be  a  result  of  the  large  II/S 
ratio  developed  in  the  extended  multivelocity  region  created  by  the  two- 
dimensional  space-charge  effects.  Comijarisons  are  made  between  the 
resxxlts  obtained  here  and  previous  res\ilts  from  the  analyses  employing 
one-dimensional  diodes  of  similar  d-c  configurations. 

The  large  discrepancy  between  the  results  from  the  two  numerical 
methods  of  noise  transport  calculations  in  two-dimensional  diodes  of 
identical  d-c  characteristics  pointed  out  the  extreme  importance  of  the 
coupling  mechanism  approximations.  Further  improvement  of  the  electron 
beam  noise  transport  theory  would  depend  on  a  thorough  understanding  of 
the  actual  coupling  mechanism  in  a  multi-dimensional  realistic  electron 
beam  model. 
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CHAFEER  I.  imODUCnON 


1.1  Noise  in  Longitudinal  Beam-Type  AnuHfiers 

Ever  since  the  Invention  of  such  electronic  communication  systems 
as  telegraphy,  telephone,  radio,  etc.,  much  effort  has  lieen  exerted  on 
the  Improvement  of  the  reliability  emd  sensitivity  of  the  receivers. 

For  Ideal  systems,  a  receiver  consists  of  one  or  more  amplifiers  and  a 
detector  designed. to  recover  the  original  signal.  Unfortunately,  It  Is 
Impossible  to  avoid  a  certain  amount  of  uncertainty  introduced,  by  those 
^perfect  amplifiers  or  detectors,  and  -^e  undesirable  fluctuations  thus 
Introduced  will' reduce  the  reliability  suad  especially  the '  sensitivity  of 
the  entire  receiving. system.  ‘ 

.  '  ■  \ 

Prior  to  the  discovery  of  such  ultra-low-noise  devices  as  masers, 
parametric  amplifiers  and  other  solid-state  amplifiers-,  the  longitudinal 
beam-type  device  has  been  -the  only  device  with  which  appreciable  ampli¬ 
fication  of  an  electronic  signal  at  microwave  frequencies  could  be- 
accomplished  wi-thout  Introducing  an  excessive  amount  of  xmcertalnty  In 
■the  signal.  Ho-wever,  because  of  -the  presence  of  velocity  and  current 
fluctuations  in  -the  electron  beam,  -the  kinetic-energy  source,  a  cer-taln 
amount  of  noise*  will  be  de-tected  at  -the  output  superposed  on  "the  ampli¬ 
fied  signal.  In  small-signal  de-tection  -the  sensltivl-ty  of  a  longitudinal 
beam  amplifier  Is  limited  by  -the  shot  noise  current  and  -the  velocity  - 
f luctua-tlons  originating  at  the  -thermionic  cathode . 


*  Noise  is  generally  considered  -to  be  the  amount  of  uncertednty  existing 
in  a  physlcELl  process;  it  may  be  represen-ted  by  a  chance  variable 
^ose  ins-tan-taneous  behavior  is  not  predictable  fixm  its  past  hls-tory. 
In  -this  dlsser-tatlon,  -the  represen-tation  of  noise  Is  limited  -to 
s-tatloneuy  ergodlc  randcm  processes  only. 
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Vco'lous  Investigations  on  noise  transport  phenonena  have  been 
caurried  out  in  the  past  starting  with  lov-frequeney  diodes  and  carrying 
to  accelerated  electron  beans.  The  one-dlnenslonal  single-velocity 
model  gave  a  beam  noise  conservation  theory  predicting  a  theoretical 
minimum  noise  figure  of  longlt\illnal  beam-type  devices  of  6-1/2  db.  On 
the  other  hand,  with  the  distribution  of  velocities  in  the  electron  beam 
taJcen  Into  account,  nearly  a  3  db  reduction  in  noise  figure  was  predicted 
by  Watkins'  multivelocity  auialysis.  Currently  the  best  experimental 
noise .performance  of  2-1/2  db  noise  figure  at  frequencies  near  3  Qc  Is 

.•'lA, 

still  unaccounted  for  by.Jlie  existing  one-dlmensloi^.  tioise  theories.  In 
view  of  the  fact  that  a  substantial  reduction  in  noise  temperature  has 
been  observed  in  amplifiers  with  a  controlled  potential  distribution  near 
the  cathode  and  that  hollow  electron  beams  are  employed  in  the  ultra-lcw- 
noise  amplifiers,  a  question  is  raised  as  to  how  the  noise  characteristics 
of  an  electron  beam  will  be  affected  by  the  multl-dlmensloncLL  potentleLL 
distribution. 

1.2  Thermionic  Emission  and  the  Origins  of  Noise  in  Electron  Beans 

Thermionic  emission  is  defined  as  the  phenonenon  by  which  some 
electrons  in  a  metsil  or  semiconductor  materleil  held  at  an  elevated  tem¬ 
perature  escape  from  the  surface  into  the  neighboring  space.  The 
emission  mechanism  can  be  represented  as  statloQEOT’  random  processes. 

At  the  cathode  plane,  two  kinds  of  fluctuations  are  observed,  namely, 
the  fluctuations  of  the  emission  current  and  the  Initial  velocities  of 
the  electrcms.  In  noise  transport  analyses  the  Interest  Is  primarily 
directed  to  the  frequency  domain,  where  the  noise  properties  are  conven¬ 
tionally  represented  and  the  spectral  densities  of  the  two  variables 
are  the  significant  characteristics. 
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At  low  freqvieacles,  where  the  transit  angles  of  the  electrons 
across  a  diode  are  negligible  compared  to  2k,  Schottky^  obtained  a 
formula  for  the  self -power  spectral  density  of  the  current  fluctuations, 
under  temperature-limited  conditions,  which  are  referred  to  as  shot  noise. 

S^(f)  •  2el^  for  f  -»0  ,  (l.l) 

where  Is  the  average  current.  Meanwhile,  an  expression  for  the  power 
spectral  density  of  the  velocity  fluctuations  created  by  shot  noise  in 
the  emission  current  was  derived  by  Rack^. 

|e|kT 

S^(f)  -  —55-^  (*♦-«)  for  f  -»0  ,  (1.2) 

o 

where  m  =  mass  of  an  electron, 

T  ■  temperature  of  the  cathode  in  degrees  K, 

0 

k  3  Boltzmann's  constant. 

The  equipartition  of  energy  in  frequency  so  indicated  is  valid  only  for 
low  frequencies  (kT  »^(d).  For  the  noise  transport  investigation 
carried  out  in  this  report,  both  Eqs.  1.1  and  1.2  are  assumed  valid. 

1.5  Survey  of  PreviovLs  Holse  Transport  Studies 

Previous  noise  studies  were  carried  out  for  one-dlmenslc»al  elec¬ 
tron  beam  models  and  they  could  be  classed  generally  in  one  of  two 
categories.  In  one  category  the  single-velocity  assumption  was  employed 
in  a  small-signal  perturbation  aneJjrsls  for  diodes  or  electron  beams  In 
which  the  thermal  velocity  spread  was  negligible;  In  the  other  category, 
the  formation  of  a  potentleO.  minimum  away  from  the  surface  of  the  cathode 
constituted  a  multivelocity  electron  beam  In  which  a  portion  of  the 
electrons  emitted  would  eventually  be  returned  to  the  cathode  because  of 
the  negative  potential  gradient. 
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1.3«1  Space-Charge  Reduction  of  Noise  In  Diodes  of  flanll  Transit 
Angle.  The  i^ncnenaa  of  space-charge  reduction  of  shot  noise  In  diodes 
was  first  treated  hy  Rack®  based  on  the  LLewellyn-Peterson*#*  equations. 
At  low  frequencies,  where  the  transit  tine  effect  Is  small,  a  shot  noise 
reduction  factor  of  0.644  was  obtained  for  space-charge-llmlted  diodes 
operated  under  Child's  law. 

Pierce^  made  use  of  the  open-circuit  assumption  and  showed  that 
the  fluctuations  In  current  density  at  the  cathode  produce  no  significant 
effects  on  the  fluctuations  at  the  anode.  An  Identical  noise  smoothing 
factor  was  derived  In  his  simplified  single  velocity  model  diode  In  which 
similarity  between  the  potential  distribution  beyond  the  mlnlimim  Anrt 
Child's  law  was  assumed. 

The  mechanism  of  noise  reduction  by  space  charge  has  long  been 
recognized  to  be  the  gating  action  of  the  potential  minimum  In  regulating 
the  number  of  electrons  passing  through.  In  a  diode  operated  under  full 
space  charge,  any  Instantaneous  Increase  over  the  average  number  emitted 
from  the  cathode  would  result  In  lowering  the  potential  minimum  Immedi¬ 
ately  and  In  turn  decrease  the  percentage  of  electrons  psuislng  throu^. 
Noise  reduction  Is  accomplished  by  the  establishment  of  the  compensation 
currents  opposite  to  the  original  disturbance. 

Thomson,  North  and  Harris^  made  an  attempt  to  evaluate  the  smooth¬ 
ing  factor  for  a  multivelocity  model  diode  assuming  Independent 
fluctuations  for  each  velocity  class  at  the  cathode.  For  a  moderate 
anode  potential,  the  noise  reduction  factor  thus  obtained  Is  the  same  as 
the  single-velocity  approximation. 

1.5.2  Noise  Transport  In  Hl^-Frequency  Diodes.  At  high  frequen¬ 
cies  the  transit  time  effect  becomes  Inqportant  and  the  noise  reduction 
mechanism  must  be  treated  from  a  somewhat  different  approach. 
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Rack?  has  extended  his  single-velocity  analysis  to  Include  the 
effects  of  finite  transit  angles.  Later  on,  Peterscm'^  applied  the 
Llevellyn-Feterson  equations  to  evaluate  the  noise  reduction  In  a 
tetrode,  and  Pierce^  employed  the  same  method  to  estimate  the  noise  In  an 
electron  beam  with  zero  velocity  at  the  potential  minimum,  realizing  that 
'such  cm  assumption  would  lead  to  the  elimination  of  the  effects  prodticed 
by  current  fluctuations  at  the  cathode  plane.  However,  his  assuiiqptlan 
ccxitradlcts  the  fact  that  the  average  electron  velocity  at  the  potential 
minimum  Is  small  but  finite,  and  the  velocity  fluctuations  there'  are  not 
the  only  source  of  noise.  .  ■ 

Whlnnery®,  Slegman  and  Watklns^®#^^  attempted  to  establish  the 
boundary  conditions  at  the  potential  minimum  of  a  multlveloclty  flow 
model  for  both  zero  and  finite  transit  angles.  An  equivalent  noise  was 
obtained  from  a  macroscopic  point  of  view)  unfortunately.  It  Is  rather 
doubtful  that  the  reduced  noise  obtained  by  the  multivelocity  aiuOysls 
would  be  -applicable  to  the  single-velocity  theory. 

1.3.3  Reduction  of  Noise  In  Beam-Type  High-Frequency  Devices  and 
the  Minimum  Noise  Figure  in  the  Single-Velocity  Theory.  An  artificial 
method  for  noise  reduction  was  first  suggested  by  Watkins^  at  the  time 
when  velocity  fluctuations  were  still  considered  to  be  the  only  effective 
noise  source  at  the  potentled  minimum.  Becatise  of  the  Invariant  char¬ 
acteristic  of  the  kinetic  potential, 

u  V 

V  .  °  ^ 

I  h  ' 

where  h  Is  the  absolute  value  of  electron  charge-to-mass  ratio,  the 
magnitude  of  the  a-c  velocity  v  may  be  reduced  In  the  case  of  a  sudden 
Increase  In  d-c  po-tentlal  at  the  velocity  fluctuatlcm  maxima.  Following 
this  argument,  Watkins  proposed  that  periodically  repeated  sudden 
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acceleratlons  asd  decelerations  placed  a  quarter  plaana  wavelength  apart 
may  In  the  limit  reduce  the  velocity  fluctuations  to  zero,  l.e.,  the 
theoretical  minimum  nolae  figure  would  he  zero  according  to  this  theory. 

Peter  and  Bloom^®j^  derived  the  transmlsslon-llne  analog  for 
space-cheurge-wave  propagation  In  unlveloclty  accelerated  electron  beams. 
The  r-f  Impedance  of  a  thermionic  cathode  may  be  matched  to  the  electrcm 
beam  Impedance  at  the  Input  of  the  amplifier  through  one  or  more  exponen¬ 
tial  transformer,  sections.  Theoretically  this  Is  an  Ideal  method  of 
noise  reduction -because  of  the  smooth  transition  of  the  d-c  potential 
along  the  electron  bew  In  an  exponential  Impedance  transformation  region. 
The  minimum  achievable  nolee  figure  of  a  one-dlmenslonal  longitudinal 
beam-type  amplifier  as  predicted  by  the  preceding  transport  analysis  Is 
approximately  6-1/2  db  with  current  fluctuations  also  taken  Into  account. 

Meanwhile,  Haus  and  Robinson^ carried  out  a  general  study  of 
noise  transport  behavior  in  a  single-velocity  electron  beam  using  the 
kinetic  power  theorem.  A  noise  invariance  theory  was  established  together 
with  the  Introduction  of  some  characteristic  noise  parameters  represent¬ 
ing  the  self-power-spectral  density  euid  the  cross -power-spectral  density. 
Haus  concluded  from  his  theory  that  correlations  between  velocity  and 
current  fluctuatlcxis.  should  remain  unchanged  In  the  propeigatlai  direction 
of  the  electron  beam. 

1.3.4  Multivelocity  Theories .  Slegman,  Watkins  and  Hsleh^’^>“ 

In  a  series  of  two  papers  described  the  formulatl(»i  and  the  method  of 
solutlcn  of  the  noise  transport  problem  In  a  multivelocity  flow  diode 
utilizing  the  density  function  analysis.  A  set  of  linearized  small- 
signal  multi-coupled  dlfferentisQ.  equations  governing  the  behavior  of 
the  density  functions  was  solved  numerically  by  a  forward  integration 
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aethod  for  an  open-clrc\ilt  diode  operated  under  Child 'a  lair.  Contzaiy 
to  the  aingle-velocity  theoriea,  a  certain  amount  of  correlation  vaa 
developed  between  the  current  and  velocity  flxictuatlona  previously 
supposed  to  be  independent  in  the  reglcai  where  oultiveloclty  flow  pre¬ 
vailed.  Neither  the  noise  parameters  nor  the  theoreticcd  mlnlaua  noise 
figure  were  found  to  be  invariant  by  the  authors  and  the  question  of  the 
existence  of  a  theoretlceil  minimum  was  again  raised. 

Vivian performed  similar  conqputatlons  on  a  space-charge-limlted 
diode  with  a  potential  minimum,  knowing  that  linear  transformation  methods 
could  be  applied  to  evaluate  the  boundary  conditions  at  the  potential 
minimum  plane.  Cue  to  the  rough  approximations  employed  in  his  discrete 
bean  model,  the  inclusion  of  the  reflected  flows  failed  to  show  appre¬ 
ciable  effects  on  the  end  resiilts. 

Following  the  invention  of  the  ultra- low-noise  electron  gun  by 
Currie®®,  extensive  studies  of  noise  propegatlon  in  low-velocity  drifting 
electron  beams  were  carried  out  both  analytically  and  numerically  by 
various  authors.  The  nonconservatlve  nature  of  the  noise  peurameters  wm 
demonstrated  by  Berghammer  and  Bloom®^  in  a  low-potential  multivelocity 
drifting  electron  beam;  and  by  means  of  the  kinetic  power  theorem  in 
electron  beams  with  appreciable  velocity  spread,  Hok®®  showed  that  the 
real  part  of  the  kinetic  power  is  not  necessarily  an  invariant.  Similar 
results  were  also  obtained  by  Shaw®®  and  others  while  solving  the 
Boltzmann  transport  equation  for  a  drifting  beam  with  square  velocity 
distributions . 

The  Monte  Carlo  approach  to  the  noise  problem  in  electron  beams 
was  first  Investi^ted  by  Tien  and  Moshman®*  to  csdculate  the  shot  effect 
at  a  plane  Just  beyond  the  potential  minimum  of  a  short-circuited  diode. 
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A  direct  slaulatloa  of  the  ealssloc  process  and  the  nechanlsn  of 
establishing  the  self-consistent  electric  fields  In  a  one-dlnenslonal 
diode  vas  performed  by  means  of  a  hl£^-speed  electronic  caqiuter. 

Instead  of  studying  the  macroscopic  behavior  of  the  smoothed  out  quan¬ 
tities  like  currents  or  density  functlcsxs,  the  electrons  were  traced  In 
phase  speuse  after  emission  which  V€ub  simulated  by  sosie  stochastic  pro¬ 
cesses  with  appropriate  distributions.  Fluctuations  In  the  electron  beam 
would  then  be  eveLIuated  at  several  planes  beyond  the  potential  minimum 
using  statistical  methods.  Appreciable  noise  reduction  was  found  at  low 
frequencies  and  a  sudden  dip  In  the  noise  reduction  factor  was  observed 
Just  belov  the  plasma  frequency  corresponding  to  the  electron  concentra¬ 
tion  at  the  potentlcQ.  minimum.  Besides  these,  no  coupling  was  observed 
between  the  noise  current  and  the  velocity  fluctuations. 

The  same  approach  was  again  einployed  by  Lambert  and  Dayem^^  to 
analyze  the  transi>ort  problem  In  an  open-circuit  diode  which  has  dimension 
and  emission  characteristics  Identical  to  Tien  and  Moshman's  diode. 

Slmllnr  to  other  open-circuit  diode  noise  anedyses,  the  dip  In  noise 
reduction  factor  no  longer  exists  and  some  correlatlcms  between  the  two 
Independent  noise  sources  were  observed. 

1.3.5  Reactive  Damping  of  Noise  In  an  Electron  Beam  with  Finite 
Diameter.  Berghammer^^  demonstrated  the  possibility  of  the  existence 
of  a  new  noise  reduction  mechanism  In  a  lew  average  velocity  electron 
beam  of  finite  diameter  using  the  square  velocity  dlstrlbutlcm.  It  Is 
believed  that  In  addition  to  the  two  space -chsurge  waves  described  by 
Ramo  and  two  decaying  waves  propagating  In  opposite  directions 

might  exist.  This  Is  similar  to  the  danqplng  mechanism  proposed  by 
Landau  In  plasma  jdiyalce. 


•9- 


In  sunmiary,  the  literature  eontaliu  numerous  nolae  transport  and 
reduction  analyses  In  electron  beans.  Hcwever,  each  of  the  above  anal¬ 
yses,  except  the  nevly  developed  reactive  dashing  hypothesis.  Is  United 
to  one -dimensional  models.  It  appears  that  a  thorou^  study  of  the 
nultl-dlmenslonal  effects  on  the  noise  characteristics  of  an  electron 
beam  has  not  yet  been  carried  out.  An  Investigation  of  this  type  would 
be  desirable  because  the  extremely  low-noise  figures  achieved  experimen¬ 
tally  are  still  unaccounted,  for  by  existing  noise  theories. 

1.1>-  Noise  Transport  In  a  Two-Dimensional  Space-Char ge-Llmlted  Diode 

During  the  past  three  years  noise  figures  lower  than  3  db  have 
been  achieved  In  longltudlnsLl  beam  amplifiers  at  S-band  frequencies. 
Because  of  the  discrepancy  between  experimental  results  and  the  existing 
one -dimensional  theories.  In  which  the  Fry-Langnulr  potential  distribution 
w€us  adopted,  it  was  suggested  to  the  writer  by  Professor  J.  E.  Rowe  that 
the  multl-dlmenslonal  effect  in  the  vicinity  of  the  potentlea  minimum 
may  better  explain  the  transport  and  reduction  phenomena  In  a  space- 
charge-limited  hollow-beam  diode. 

1.4.1  Two-Dlmens lonal  Model.  In  the  following  chapters,  the 
noise  transport  characteristics  of  a  theoretlced.  multl-dlmenslonal 
high-frequency  diode  will  be  studied  by  two  different  numerlced  methods. 
Since  the  noise  characteristics  of  thin  annular  beams  are  found  to  be 
superior  to  the  ccxiventlonal  pencil  beams,  a  mathematlccd.  model  Is  chosen 
to  simulate  a  hollow  beam  with  a  small  thickness -to-radlus  ratio.  If  all 
the  veu-latlons  In  the  circumferential  direction  are  neglected,  each  snail 
section  of  this  actual  electron  beam  may  be  expanded  and  developed  Into 
a  thin  ribbon  beam  of  Infinite  width  aa  shown  In  lig.  1.1,  resulting  In 
a  reduction  by  one  dimension.  Noise  analyses  are  to  be  carried  out  for 
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the  two-dimensional  electron  heam  model  thus  obtained,  flowing  between 
two  parallel  plane  electrodes. 

1.4.2  Methods  of  Investigation.  Among  the  methods  of  emalysis 
adapted  by  previous  workers  only  the  density  function  method  and  the 
Monte  Cetrlo  techniques  are  valid  In  the  cathode-potential  minimum  region 
where  reflected  flows  exist.  Both  methods  will  be  employed  in  the 
following  chapters  to  solve  the  problem  of  noise  propagation  for  a  two- 
dimensional  discrete-beam  diode  with  identlceuL  emission  characteristics 
to  the  one -dimensional  diode  selected  by  Tien  and  Moshman  and  Lambert  and 
Dayem. 

In  the  density  function  method,  the  analysis  of  noise  or  signal 
propagation  along  an  electron  beam  depends  on  the  simultaneous  solution 
of  a  set  of  linearized  multl-coupled  a-c  transport  equations  corresponding 
to  various  velocity  classes  In  the  discrete  beam.  Numerical  solutions 
may  be  obtained  with  the  aid  of  a  high-speed  electronic  computer.  Because 
of  the  linear  natxure  of  the  working  equations,  the  well  established  linear 
transformation  theories  may  be  applied  to  establish  the  boundary  condi¬ 
tions  for  the  transmitted  flows  at  the  potential  minimum,  making  It 
possible  to  evaluate  the  effects  due  to  the  electrons  returning  to  the 
cathode.  Coupling  between  streams  not  experiencing  the  same  d-c  electric 
field  Is  accomplished  by  the  a-c  electric  field  associated  with  the 
convection  current  In  the  diode.  It  Is  believed  that  correlations  between 
velocity  and  ciirrent  fluctuations  will  be  developed  through  the  a-c 
coupling  electric  fields. 

The  Monte  Carlo  method  of  noise  transport  calculation  Is,  on  the 
other  hand,  extremely  simple  in  mathematical  formulation.  Random  numbers 
are  generated  corresponding  to  the  stochastic  emission  process  and  each 
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electron  emitted  is  traced  In  phase  space  throuc^out  "by  n\merieal  Inte¬ 
grations.  Both  linearization  and  singularity  extractions  are  avoided  and 
hence  the  true  nature  of  electrons  Inside  a  beam  Is  preserved.  Coupling 
between  streams  Is  achieved  by  the  open-circuit  assumptlcxi.  Potential 
distribution  Inside  a  two-dimensional  diode  with  certain  boundary  condi¬ 
tions  are  evaluated  by  taking  the  tine  average  of  the  self-consistent 
potentleil  profile  under  space-charge-llmlted  conditions  and  the  instan¬ 
taneous  electric  fields  are  computed  by  the  Image  method  sunmlng  over  all 
electrons  Inside  the  diode. 

Noise  parameters  at  some  planes  where  only  transmitted  flows  exist 
are  evaluated  and  the  results  obtained  by  these  two  methods  of  analysis 
are  compared  with  previous  ones  derived  using  the  one -dimensional  theory. 


CHAPTER  II.  DENSm-PtWCTION  TOBMJIATIOW 

2.1  Introduction 

In  the  following  sections  >  the  Boltzaaxu  transport  equation  Is 
formulated  In  terns  of  the  density  function,  the  electric  fields  and 
the  operating  frequency  for  a  space-charge-linlted  thermionic  diode  In 
two-dimensional  space.  The  density  function  Is  defined  as  the  density 
of  electrons  in.  a  slx-dlmenslonal  phase  space  in  which  electric  charges 
and  masses  are  considered  to  have  a  continuous  probabilistic  distribu¬ 
tion.  Instead  of  stvidylng  the  microscopic  movements  of  Individual 
electrons,  the  density- function  method  of  noise  transport  studies  deeds 
with  the  average  macroscopic  behavior  of  the  density  functions  which  Is 
In  tiu:n  governed  by  the  Boltzmann  tremsport  equation.  Discrete  velocity 
classes  and  a  discrete  beam  model  cure  Introduced  to  reduce  the  nonlinear 
transport  equation  to  a  set  of.  llneeu*  smedl-slgned  dlfferentled  equations 
mutually  coupled  through  the  space-charge  electric  fields.  Linear  trans¬ 
form  techniques  are  thus  applicable  In  solving  the  boundary-value 
problems  with  uncorrelated  (independent)  inputs. 

The  method  of  extraction  of  singularities  Is  described  In  detail 
and  the  source  of  the  singularities  Is  also  Investigated.  Special 
attention  Is  directed  to  the  disappearance  of  Watkins's^  low-frequency 
limit  In  the  transport  equation  for  a  two-dlmensloncd  electron  beam. 

Details  of  the  density-function  formulation  are  Included  and  the 
expressions  for  the  noise  parameters  In  terms  of  the  density  functions 
are  developed.  Haus's  noise  theory  for  small-signal  longltudlzial-beaffl 
devices  Is  eaqiloyed  to  evaluate  the  minimum  achievable  noise  figure 
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for  an  anpllfler  In  vhlch  Important  quantities  depend  upon  tvo  space 
dl»ensl(»is. 

2.2  A  Discrete  Beam  Model  Tuo-Dlmenslcoal  Diode 

In  the  following  noise  transport  analyses,  the  actual  two- 
dimensional  electron  beam  Is  replaced  by  a  mathematlccdly  discrete  beam 
model  to  facilitate  solution  by  digital  compter  techniques.  In  two- 
dimensional  models,  eQl  vcurlatlons  in  one  of  the  transverse  directions 
are  Ignored.  For  a  z-dlrected  electron  beam  of  finite  thickness  In  the 
y-dlrectlon,  no  variation  of  any  kind  Is  allowed  frcn  x  *  -os  to  x  >  oo. 

If  the  ribbon  beam  Illustrated  In  Fig.  2.1  Is  divided  Into  smsdl  sections 
of  thickness  the  total  current  per  unit  length  In  the  x-dlrectlon  may 
be  replaced  by  a  current  sheet  of  equsd  magnitude  at  the  center  of  the 
corresponding  Increment  dy.  The  resulting  model  thus  consists  of  layers 
of  ctirrent  sheets  of  Infinitesimal  thickness  located  sty«0,  y^t  dy, 
y  ■  ±  2  ...  etc.,  separated  by  a  small  Increments^,  distance  dy.  The 

mathematics^,  electron  beam  model  may  be  made  em  arbitrarily  close  approxl- 
matlon  to  the  reeil  ribbon  electron  besun  flowing  between  two  electrodes 
of  a  diode  lnnersed  In  a  strong  socls^  magnetic  field  by  reducing  the 
magnitude  of  each  current  sheet  suid  Increeislng  the  number  of  current 
sheets  per  unit  length  In  the  y-dlrectlon. 

The  discrete  beam  method  used  In  the  suialysls  gives  good  approxl- 
metUons  sub  long  as  dy  la  smsdl  coopsured  to  the  wavelength  of  the  highest 
frequency  of  Interest.  In  the  density  function  method  of  Investigation, 
only  the  z-dlrected  current  density  is  of  Importsmce,  the  width  of  the 
electron  beam  In  the  x-dlrectlon  Is  tsdwn  to  be  unity  at  sdl  times.  Due 
to  geometric  symnetry,  sheet  beams  of  the  same  distance  from  the  center 
of  the  ribbon  beam  should  be  Identicsd.  in  d-c  characteristics. 
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2.3  AsaumptioPB 

Certain  asaunptlonB  are  Bade  to  alapllfy  the  problem  vhieh  la  to 
be  aolved  numerically  on  a  digital  ccnputer. 

1.  The  aneOyaia  la  two  dlmenalonal>  yet,  no  tranaverae  motion  of 
electrona  la  eJJ.Otfed.  An  infinite  magnetic  field  parallel  to  the 
direction  of  the  electron  beam  ia  aaauBted. 

2.  There  are  no  direct  colUalona  between  electrona;  thia,  together 
with  the  prevloua  aaavmqittlon,  limlta  the  coupling  mechanlam  to  apace> 
cheurge  forcea  alone. 

3.  NonrelativlBtic  aaaua^tlon  prevaila. 

4.  The  initial  velocity  diatrlbution  at  the  cathode  la  taken  to 
be  the  Rayleigh  diatrlbution. 

3.  Aa  shown  in  Fig.  2.2  the  beam  model  is  divided  into  many  velocity 
claeaea;  the  number  of  velocity  claaaea  is  large  enouedi  to  inaure  a  amall 
velocity  spread  within  each  velocity  claas  compeured  to  its  reapectlve 
mean  velocity. 

6.  The  linearized  small-aignal  Boltzmann  transport  equation  ia  uaed 
in  the  anstlyaia. 

7.  The  fluctuation  at  the  cathode  plane  in  any  one  velocity  claaa 
is  assumed  to  be  uncorrelated  with  the  fluctuation  in  any  other  velocity 
class. 

8.  Those  electrons  returning  to  the  cathode  before  reaching  the 
potentica  minimum  produce  no  effect  on  cathode  emission.  There  is  no 
secondary  emission  at  either  of  the  two  electrodes. 

9.  There  is  no  accumulation  of  electrons  at  the  potential  minimum. 
The  probability  that  an  electron  ahould  possess  Just  enough  energy  to 
reach  the  potential  mlnlBiaii  and  stay  there  la  zero. 


FIG.  2.2  RETARDIl*}  FIEID  DIODE  VfITH  DISCRETE  VELOCITI  CLASS.  VELOCOT 
VS.  DISTAHCE  FOR  EACH  LAYER  OP  THE  ELBCTROH  BEAM. 
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10.  The  d-c  potential  dlstrlhutim  la  obtained  from  the  time 
average  of  the  potential  distributions  from  the  Monte  Carlo  calculatloas. 

2.5.1  Transport  Equations.  The  transport  equations  written  In 
terms  of  the  density  functions^  the  electric  fields  and  the  operating 
frequency  eure  obtained  from  the  combination  of  LlouvlUe's  theorem,  the 
force  lavs  and  the  small-signal  eissumptlon.  The  density  function 
^y}Z,y,z,t)  for  a  two-dimensional  model  Is  defined  such  that  [F(y,z,y, 
z,t)  dy  dz  dy  dz]  represents  the  n\mber  of  electrons  found  between  y  and 
y  +  dy,  z  and  z  -t-  dz  with  y-directed  velocities  between  f  and  f  and 
z-dlrected  velocity  components  between  z  and  z  +  di  at  time  t.  According 
to  LlouvlUe's  theorem,  the  density  function  P(y,z,5^,i,t)  in  phase  space 
Is  Invariant  with  respect  to  tlmej  l.e.,  the  total  derivative  of  F(y,z, 
y,z,t)  with  respect  to  time  Is  IdentlceLLly  zero. 

dF(y,z,y,z,t)  _  3F(y,z,y,z,t)  dF(y,z,y,z,t)  .  8F(y,z,y,&,t)  „ 

dt  dt  dy  dy  ^ 
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For  a  linear  small-signal  perturbation  analysis,  the  density 
function  as  veil  as  the  electric  field  Is  assumed  to  be  separable  into  a 
d-c  component  and  a  sinusoidal  time-varying  component  In  such  a  manner 
that  the  magnitude  of  the  r-f  component  Is  very  small  compared  to  that  of 
the  d-c  part. 

F(y,z,z,t)  =  F^(y,z,z)  +  F^(y,z,z)  e’^'*^  ,  (2.5) 

vhere 

T^(y,z,z)  ^  F^(y,z,z)  F^^(y,z,z) 

and 

F^Jy,z,z)  «  1  . 


E,(y,z,t)  =  E^^(y,z)  +  E,^(y,z) 


iz' 


(2.4) 


\rtiere  <b  is  the  radian  frequency  of  Interest.  The  subscripts  "o"  and  "i" 
denote  the  d-c  and  the  a-c  parts  respectively. 

The  nonlinear  partial  differentlsd.  equation,  Eq.  2.1,  may  now  be 
reduced  to  two  linear  transport  equations  which  govern  the  behavior  of 
the  d-c  and  the  a-c  components  of  the  density  fvmctlons  by  neglecting  the 
cross  product  term  Involving  two  a-c  quantities. 


^o(y,z,z)  SF  (y,z,z) 

— S - — si —  -  0 


(2.5) 


amd 


SF  (y,z,z)  QF  (y,z,z) 

JaiF^(y,z,z)  +  z - ^ - hE^^(y,z) 


hi 


§2 


=  0  ,  (2.6) 


Equation  2.6  Is  the  basic  transport  equation  for  the  multivelocity  elec¬ 
tron  beam  In  which  the  third  term  describes  the  effect  due  to  the  potential 
gradient  while  the  last  term  specifies  the  amount  of  coupling  between 
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varlous  velocliy  classes  and  adjacent  electron  streams  throu^  the  r-f 

electric  field  E  .  For  a  two-dlmensloxial  model,  E  may  be  expressed  In 
iz  'xz 

terms  of  the  density  functions  (see  Appendix  A.l) . 

00  00 

E^2(y,z)  f  f  P’  i'Pi(y',z,z')e'^’|y'y’l  dy'  dz'  , 

^  .00  .00 

(2.7) 

where  the  y' -Integration  is  actually  taken  across  the  thickness  of  the 
electron  beam. 

E  (y,z)  Is  the  a-c  electric  field  at  (y,z)  due  to  the  existence 

XZ 

of  an  a-c  density  function  at  (y',z)  which  belongs  to  the  u'  ■  z' 
velocity  class. 

P'  = 

For  a  drifting  beam  without  coupling  between  density  functions  of  dif¬ 
ferent  classes,  a  solution  for  Eq..  2.6  may  be  obtained  by  setting  E 

iz 

to  zero. 

Fi(y^z,t)  =  Fj(y)  ,  (2.8) 

where 

P  “  V  . 
z 

This  is  designated  as  the  solution  of  the  homogeneous  transport  equation 
in  which  signals  travel  along  the  stream  with  only  a  change  in  phase 
angle. 

2.3.2  Excess  Energy  Parameter.  In  noise  transport  studies,  a  new 
parameter  Vf,  the  excess  energy  parameter,  is  often  introduced  to  replace 
the  velocity  parameter  z  for  Identification  purposes. 

W  ■  z®  -  2h[V(y,z)  -  V'(y)]  , 


(2.9) 


-21- 


wliere  V(y,z)  (see  Appendix  A.Il)  is  the  d-c  potential  at  (y,z)  with 
respect  to  the  cathode,  V'(y)  Is  the  potential  minimum  at  plane  y,  and 
V'(y)  Is  generally  negative.  W  represents  the  amount  of  excess  energy 
of  an  electron  over  Its  potential  energy  In  the  electrostatic  field;  It 
Is,  In  the  steady  state,  determined  solely  by  the  InltlsLL  velocity  of  the 
associated  electron  and  Is  Invariant  throughout.  The  Introduction  of  W, 
a  scalar  pEurameter,  to  replace  a  vector  variable  z  however  raises  some 
complications  In  the  region  where  reflected  flows  exist.  Because  the 
amount  of  excess  energy  possessed  by  an  electron  of  a  reflected  flow 
class  Is  Independent  of  Its  direction  of  propagation,  a  plus  or  minus 
sign  must  be  attached  to  the  differential  of  the  parameter  W  Indicating 
whether  the  group  of  electrons  Is  traveling  In  the  forward  or  reversed 
direction,  l.e.,  towasrd  the  anode  or  cathode. 

dW  =  ±  2z  dz  .  (2.10) 


2.3.3  Dimensionless  Parameters .  The  following  dimensionless 
parameters  are  brought  together  with  Eqs.  2.7  through  2.10  to  reduce  the 
a-c  transport  equation  into  dimensionless  form. 


n(W)  S 

^C 

^  y 

6  ^  pa®/*  z 

g  ^  _<S_ 
paV* 


Potential  Parameter, 


Real  Space  Coordinate 
Parameters , 


Frequency  Parameter, 


A 

m 


P(y,z,W) 


Density  Function  Parameter, 


Excess  Energy  Pcurameter, 


w  «  aW 


(2.11) 
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vhere  a  *  in/2]iT  . 


/ell 


m 


mass  of  an  electron > 

average  current  density  of  the  electron  heam, 
cathode  temperature  In  degrees  K, 

Boltzmann's  constant, 

■  dielectric  constant  of  free  space. 

The  transport  equation  may  now  be  written  as  (see  Appendix  A. Ill): 


c 

k 


r  ^  1  p  r 

“L  dw  J  lj8>/;r  J  J 


•  exp  r  y  --Si^C’l^l  dw  dC'l  .  (2.12) 

*-  Vw(^',5)+tj(SM)-I  Vw(C’,6)+n(&',S)  ^ 


The  upper  signs  are  for  the  forward  going  electrcxis  of  the  density  func¬ 
tion  parameter  Immediately  following,  while  the  lower  signs  are  for  the 
reflected  flows. 


2.k  Velocity  Claeses 

The  concept  of  velocity  classes  Is  Introduced  to  split  the  trans¬ 
port  equation  Into  a  set  of  llneetr  small-signal  multl-coupled  dlfferentlciL 
equations  In  the  multivelocity  region.  As  shown  In  Fig.  2.3,  all 
electrons  emitted  with  cm  excess  energy  parameter  between  - Av^/2 
andw^  ••■Av^/2  are  considered  to  be  In  the  velocity  class  characterized 
by  v^.  If  the  rms  fluctuatloi  within  each  velocity  class  Is  wwan  cod- 
peured  to  the  mean  velocity,  a  delta  function  can  be  placed  at  the  center 
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of  each  of  the  velocity  classes  to  represent  their  respective  ms 
fluctiiatlons.  This  means  of  Incorporating  the  fluctuations  Into  the 
mean  velocity  of  every  velocity  class  may  be  accomplished  by  choosing  a 
large  number  of  classes  to  nmke  certain  of  a  negligibly  small  deviation 
from  the  mean  In  each  class.  Since  no  correlatlcms  between  Initial 
fluctuation'’  from  different  velocity  classes  are  postulated,  the  sum  of 
the  mean  square  current  fluctxiatlons  at  the  cathode  plane  Is  taten  to  be 
the  full  shot  noise  Input. 

One  of  the  most  controversial  problems  brought  forth  by  the  dis¬ 
crete  velocity  class  model  Is  the  existence  of  the  zero  velocity  class 
whose  characteristic  excess  energy  parameter  U  Is  zero.  In  a  space- 
charge  -limited  diode  with  the  formation  of  a  potential  minimum,  only 
those  electrons  emitted  with  W  >  0  have  enough  energy  to  go  pewt  the 
potential  minimum.  W  ■  0  merely  implies  that  the  group  of  electrons 
possesses  Just  enough  energy  to  reach  the  potential  minimum  and  stay 
there.  In  an  actual  diode,  the  probability  measure  of  l  electrons 
having  W  »  0  Is  zero.  The  accumulation  of  electrons  Is  therefore  not  a 
true  representation  of  what  really  happens.  Velocity  classes  should  be 
rearranged  so  that  electrons  sure  either  In  the  lowest  W  >  0  class  or  the 

lowest  V  <  0  class.  Similar  surgument  nmy  also  be  applied  to  the  W  ■  -n 

'c 

velocity  class  In  which  electrons  sure  emitted  with  zero  velocity  from 
the  cathode. 

2.$  Extraction  of  Singularities  (see  Appendix  A.IV) 

In  density  function  ansULysls,  the  step  function  discontinuities 
in  the  d-c  psurt  of  the  density  function  F^(w)  would  any  function 
Involving  the  term  &F^(W)/SW  nonauisdytlc  at  one  or  more  vsU.ue8  of  W. 

The  presence  of  this  Ill-behaved  fact(»r  often  Introduces  some  undesirable 
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complications  In  solving  the  transport  equation  hy  certain  nvsaerlcsd 
techniques .  Fortunately,  a  method  vaa  derived  hy  Sle^uan^^  hy  'which  the 
slngvilarlty  factor  In  the  a-c  density  functions  may  he  extracted  leaving 
only  analytic  functions  In  the  Boltzmann  equation. 

The  a-c  density  fvinctlon  may  he  separated  Into  two  parts  which 
will  he  called  the  homogeneotis  part  and  the  response  part  respec'tl'vely . 

The  homogeneous  or  kinematic  part  consists  of  a  Dirac  delta  fvmctlon 
corresponding  to  the  current  fluctmtlons  originating  from  'the  'thermionic 
emission  process .  This  part  undergoes  no  change  In  magnitude  and  propa- 
ga'tes  along  -the  electron  beam  under  d-c  effects  alone.  The  response  part 
of  -the  a-c  densl-ty  function  is  -the  result  of  in-tercoupling  -throu^  the  a-c 
electric  fields.  If  the  singularities  of  "the  response  parts  are  postu- 
la-ted  to  he  of  "the  same  nature  as  a  new  variable  p^(^,|,w)  may  he 

defined  for  a  unit  a-c  current  Input  at  -the  veloci-ty  class  characterized 
hy  an  excess  energy  parame'ter  w^. 


r  (5,i,v)  I  -je  I 

"  26(w-Wj^)e  L" 


"ST 


,  (2.13) 


where  p^(?»l,w)  does  not  con-tain  any  singularity, 

Iq(5)  is  'the  average  current  densi-ty  of  "the  beam  layer,  cen-tered 
at  e 


®i  =  8 


JL 


"±"  signs  are  for  -the  forward  and  reflected  flow  respectively. 
6(w-Wj^)  is  defined  as  the  Dirac  del'ta  function  wl-th  -the  following 


proper-ties; 


6(w-w^)  ^  0 


for 


w  ^  w^. 


CO 


for  w 
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Wj^+€ 


J  6(w-Wj^) 


dw  ^  1  , 


Wj_-e 


w^+e 


J'  8(w-w^)  P(w)  dv  =  F(wj^)  , 


w^-e 


From  Eq.  A. 24 


SPq(w) 

SvT 


=  -  2 


^5) 


r  -]  p  5P  (w) 

|^8(w^)  -  S(v^)J  J  -§^dw 


=  -  2 


'o(0 


{/ 


+  n)  -  s(w^  +  n^,) 


dv 


0 

-  f  [^B(w_  +  \)  -  S(w_  +  \)  dv_|  , 


dw  =  -  dv  • 

+  -  ' 

•vrtiere  S(w^)  =  a  unit  step  function, 

S(w_|_)  =  1  for  v^  =  0,  z  >  0 

S  ( )  =  0  for  ®  • 

Equation  2.12  may  now  be  expressed  in  terms  of  the  response  part 
p^(5,6,w)  of  the  a-c  density  functions.  In  the  preceding  section,  an 
argument  was  presented  showing  that  the  population  of  electrons  charac¬ 
terized  by  the  energy  parameters  w  =  and  w  =  0  are  zero.  Ihe  current 
fluctTiattons  brought  forth  by  the  singularities  In  at  these 

two  VEiLues  of  v  are  actuauLly  included  in  their  respective  adjacent  vel¬ 
ocity  classes.  After  some  simplification,  the  transport  equations  may 
be  written  for  an  lxq;nat  at  velocity  class  w^^  as 


'/w+n(5,t) 
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± 


>v) 


± 


If 


•  00 


-’>c 


8(w-w^) 

Tw^JTFTT) 


exp 


-je. 

e  ^  dw  dC 


3^  y  y  exp  [  -  glS-e’l  1 


J4 


•00  0 


MV ,t)+y\U' A)  -I  Vw(C’,6)+ti(CSS) 


exp 


[-w(cso] 


•dwdC '  -  J  y'  exp  ^ - -  -  -l  ^~^-[ 


J4.  TiT 


>w) 


-Ir 


^w(c',j)+Ti(5’,e)  J  >RF7rWFT) 


•exp  |^-w(C’p6)j  dw  dC'l  .  (2.11^) 

where  Is  the  potential  parameter  at  the  cathode. 

2.6  Boundary  Conditions  for  the  A-c  Density  Functions 

The  following  conditions  are  attached  to  the  a-c  transport 
equations : 

1.  Full  thermionic  emission  shot  noise  Is  assvnaed  at  the  cathode 
of  the  diode,  and  the  response  part  of  the  density  function  for  all  the 
forward  going  electrons  Is  Identically  zero  at  the  cathode  plane  for  all 
velocity  cleuBses. 

2.  C(»itlnulty  of  the  a-c  density  functlcms  Is  assumed  at  the 
potential  minima  as  well  as  the  points  where  the  w  >  0  velocity  classes 
reverse  their  direction  of  propagation.  The  a-c  density  function 
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assoclated  with  the  forward  going  electrons  Is  equal  to  that  Msociated 
with  the  baclGward  going  electrons  of  the  sane  velocity  class  at  the 
point  of  reflection. 

3.  The  Initial  conditions  specified  cure  Insufficient  If  numerlccd. 
solutions  are  to  be  evaluated  by  the  step-by-step  forward  Integration 
nethod  stcurtlng  from  the  cathode  plane,  because  of  the  difficulties 
encountered  In  an  attempt  to  deteimlne  uniquely  the  response  part  of  the 
a-c  density  functions  associated  with  the  backward  going  electrcms  of 
the  reflected  flows  at  the  cathode.  Another  boundcury  Is  selected  as  the 
starting  point  Instead.  As  a  natter  of  convenience,  the  integrations 
will  be  carried  out  beginning  at  the  potential  minima  and  the  points  of 
reflection.  A  method  by  which  the  boundary  conditions  of  the  a-c 
density  functions  are  determined  at  the  newly  selected  boundary  will  be 
fully  described  in  Section  2.8. 

2.7  Working  Equations 

A  mathematically  discrete  beam  model  which  can  be  made  to  approxi¬ 
mate  the  d-c  density  function  arbitrarily  well  in  the  velocity  domain  Is 
introduced  so  that  the  Integrations  In  the  a-c  transport  equation  may 
be  replaced  by  summations  during  numerical  computations.  The  discrete 
beam  model,  as  shown  in  Fig.  2.4,  consists  of  a  number  of  velocity  classes 
each  characterized  by  an  excess  energy  peurameter  w^  and  a  d-c  density 
function  parameter  where  Pq(w)  is  assumed  to  be  a  constant  for 

all  electrons  Inside  any  velocity  class.  In  Section  2.9  the  a-c  current 
and  the  kinetic  potentlcd  are  shown  to  be  proportional  to  the  first  and 
third  velocity  moments  of  the  d-c  density  function,  respectively.  There¬ 
fore,  a  satisfactory  discrete  beam  model  may  be  obtained  provided  Its 
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first  and  third  velocity  monents  are  equal  to  the  sane  moments  of  the 
half-Mexvelllan  distribution  for  each  velocity  class. 

From  Eq.  2.l4,  working  equatl(»is  for  electrons  going  In  both 
forward  and  reversed  directions  are  derived  with  a  unit  Input  to  the 
homogeneous  peurt  of  the  a-c  density  function  at  layer  ^  and  velocity 
class  v^. 


Similarly,  with  a  unit  Input  to  the  response  part  of  the  a-c  density 
function  at  layer  and  velocity  class  w^^,  the  working  equations  are; 


as 


(2.15c) 


where 
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and  the  superscripts  or  cure  Indications  of  the  propagation  direc¬ 
tions,  1*  defined  as  the  response  part  of  the  a-c 
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denslty  function  at  (Cj>0  telonging  to  velocll^  class  due  to  a  \inlt 
Input  to  the  homogeneous  part  of  the  a-c  density  function  at  layer  ^ 
and  velocity  class  v^,  ^ J ^ ^ ^  response  part  of  the  a-c 

densl-ty  function  of  velocity  class  Wj^  at  (5j>0  due  to  a  unit  Input  to 
the  response  part  of  the  a-c  density  function  at  layer  and  velocity 
class  w^,  ^ 


vhere  0  =  the  phase  shift  per  vudt  time  of  the  homogeneous  psurt  of 

the  a-c  density  function, 

X(ti,w)  =  lifw<0,  w+Ti  =  0.02, 

=  0  otherwise , 

Y(ti»w)  =  lifv>0orw<0  and  w  +  tj  =  0.02. 

Only  fifteen  discrete  velocity  classes  are  employed  hecatise  of  the 
limitation  in  computer  time  eind  expense .  For  a  Taeam  of  m  layers  and  n 
velocity  classes,  a  [2  x  n  x  m]  x  [2  x  n  x  m]  sq.'uare  matrix  will  have  to 
be  solved  for  the  initial  conditions  at  the  new  boundary.  (The  factor  of 
two  comes  from  the  complex  nature  of  the  density  functions.)  The  amount 
of  computing  time  recLulred  goes  up  as  the  square  of  the  number  of  velocity 
classes . 

2.8  Linear  Treuisf ormati ons 

The  linear  characteristics  of  the  working  equations  and  the 
simultaneous  emission  assumption  allow  the  application  of  the  superposi¬ 
tion  theorem  at  the  cathode  to  determine  the  boundary  conditions  at  the 
newly  chosen  botmdary.  For  the  m-layer  n-beam  model,  the  set  of  coupled 
transport  equations  defines  a  unique  linear  transformation  between  two 
sets  of  2n  X  m  linearly  Independent  functions,  i.e.,  the  homogeneous  and 
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the  response  parts  of  the  a-c  density  functions  associated  with  the 
forward  going  electrons  at  the  cathode  plane  and  those  at  the  new  bound¬ 
ary.  Due  to  the  second  boundary  condition,  the  a-c  density  functions 
associated  with  flows  of  opposite  directions  are  nutually  dependent  In 
natxire.  Now  If  n  x  m  sets  of  linearly  Independent  Inputs  consisting  of 
"ones"  and  "zeros"  for  the  response  parte  are  Injected  at  the  x>otentlal 
minima  and  the  reflection  points,  the  boundary  conditions  of  the  a-c 
density  functions  at  the  cathode  plane  must  be  satisfied  by  a  linear  com¬ 
bination  of  the  n  X  m  sets  of  transforms  thus  obtained  plus  the  transform 
of  a  full  shot  noise  input  (homogeneous  part)  at  the  potential  minima. 

In  Appendix  A. VI  the  weighting  factor  R(Cju>Wj^)  for  the  homogeneous  part 
of  the  a-c  density  function  at  layer  m  and  velocity  class  1  is  evaluated 
using  the  shot  noise  formula  by  Schottky  and  the  uncorrelated  input 
assumption  at  the  cathode  plane  for  different  velocity  classes. 

For  each  input  in  the  homogeneous  part  of  the  a-c  density  func¬ 
tion  at  layer  m  and  velocity  class  1,  a  weighting  factor  G(  5^,w^) 

for  the  response  part  of  the  a-c  density  function  P^(  at 

the  potential  minimum  may  be  obtained  from  the  following  equation.  At 
the  cathode  plane  where  the  algebraic  sum  of  the  response  parts  in  each 

velocity  class  is  identically  zero, 

2  a 

n«-2  f=-7 

The  weighting  factor  G  needed  to  satisfy  the  boundary  conditions  at  the 
cathode  then  specifies  the  initial  conditions  at  the  new  boundary 
selected  for  the  transmitted  flow  Induced  by  the  unit  noise  current 
input  at 


2.9  Congputatlon  of  Noise  Parameters  euad  the  Minimum  Noise  Figure 

Noise  parameters  expressed  In  terns  of  the  self-pover  spectral 
density  (SPSD)  euid  the  cross-power  spectral  density  (CPSD)  of  the 
z-dlrected  a-c  current  and  voltage  fluctuations  according  to  the  small- 
signal  unlveloclty  theory  of  Haus  and  Robinson  are  employed  In  the  noise 
transport  analysis  In  a  two-dimensional  multlveloclty  discrete  beam  model 
provided  that  the  equivalent  current  fluctuations  and  the  kinetic 
potential  are  properly  defined. 

2.9.1  Current  Fluctuations.  The  noise  current  density  at  any 
point  (5>6)  per  imlt  bandwidth  associated  with  a  unit  current  Input  at 
the  w^  velocity  class  Is  defined  such  that 

00 

=  -  I  y'  P^(C,l,w;C„,w^)  dw  ,  (2.17) 

0 

or  In  terms  of  dimensionless  parameters. 


(2.18) 


k=i  n=-2  i*-7 


where  H  and  R  Indicate  the  homogeneous  Input  and  response  Input  respec¬ 
tively.  For  full  shot  noise  Input  at  the  cathode,  the  varleince  of  the 
noise  current  at  any  point  downstream  Is,  by  the  superposition 

theorem. 


-55- 


a  e 

<i®(Cj,S)>  -  Y,  ^  Aw^(AC)® 

Om-Z  l«-7 

‘[^iH^  ■*■  }  *  (2*19) 


where  "*"  denotes  complex  conjugates. 

2.9«2  Kinetic  Potentleil.  In  noise  transport  analyses,  the 
velocity  fluctuation  In  an  electron  beam  Is  often  replaced  by  a  variable 
which  has  the  dimension  of  voltage  and  Is  designated  as  the  kinetic 
potential.  From  an  energy  standpoint,  may  be  regarded  as  the  kinetic 
energy  density  expressed  In  electron-volts  and  it  Is  therefore  continuous 
through  a  velocity  Jump,  l.e.,  a  sudden  acceleration  or  deceleration  of 
the  electron  beam.  In  a  single-velocity  perturbation  analysis,  the 
kinetic  potential  Is  defined  to  be 

Vi  -  -  ,  (2.20) 

where  u^  Is  the  d-c  or  avereige  velocity  of  the  electrons  in  the  beam,  and 
Vj^  is  the  a-c  velocity  component. 

In  a  multivelocity  region,  the  expression  for  the  kinetic  potentleQ. 
Is  not  nearly  as  unique.  It  can  be  shown  that  a  direct  substitution  of 
the  density  function  Into  Eq.  2.20  does  not  yield  a  quantity  with  the 
invariant  characteristic  over  potential  gaps.  In  Appendix  A.V,  Is 


shown  to  be 
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V^(y,z)  -  J  53  rP^(y,z,£)  +  F^(y,z,i)j  dz 

0 

00 

where  1^(0  and  l(0  are  the  d-c  current  density  and  the  total  current 
density  at  layer  ^  respectively.  After  some  algebraic  manipulations, 


£3p 


j(y,z,£)d£j- 


(2.21) 


00 

■  i,hiJs)A5  I 


(2.22) 


where 


for  half -Maxwellian  velocity  distributions. 

It  can  be  seen  that  Eq.  2.22  will  not  reduce  to  Eq.  2.20,  the 
univelocity  definition  of  kinetic  potential,  unless  the  velocity  spread 
involved  is  negligibly  small.  Similarly,  a  direct  substitution  of  the 
density  functions  into  Eq.  2.20  will  not  give  a  satisfactory  equivalent 
kinetic  potential  for  a  multivelocity  electron  beam,  since  the  quantity 
thus  obtained  will  not  be  invariant  over  a  velocity  gap. 

In  terms  of  dimensionless  parameters,  the  kinetic  potential  at 
(^j,6)  due  to  a  unit  homogeneous  input  at  ^  and  the  w^  velocity  class 
is 


-klT  I  r'  r— I  -V, 

"  “e“  iTfy  L  L 

°  kal 


Aw 
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and  for  unit  response  Input, 


8 


i— 7 


The  variance  of  the  kinetic  potential  at  Inside  a  diode  with  full 

shot  noise  input  at  the  cathode  is  then 
a  a 


<Vj(tj,S)>-  ^  ^  {RCE„,v^)iWj 


m=-2  l=-7 


(2.24) 


and  the  covariance  between  the  current  and  voltage  fluctuations  is 
defined  as 

2  B 

. .  ■  E  E  ,w^)  Aw^  AC 

m=-a  i=-7 


(2.25) 
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Slightly  different  expresslone  for  the  fluctuations,  their 
i^spectlve  variances  and  the  coveurlance  will  be  obtained  If  all  layers 
are  taken  Into  account  simultaneously.  Comparisons  may  be  made  between 
the  noise  parameters  of  the  Individual  layers  and  those  of  the  whole  beam. 
Noise  current  densities  at  $  are 


2  a 


js-2‘  k»i 


+  -|^e  (2.26a) 


and 


2  8  2  8 


“  5°  Z  Z  .  Z  Z 


ja-2  k»i  n=-2  i»-T 


*V  ~i 

‘  (2*26b) 


Similarly,  the  voltage  fluctuations  are  written  as 
-KT  ^  ® 


J.:-2  k»i 


n,  -J04  1 

+  (w^  -  l)e  Av^  (2. 


27a) 


and 


ZZl  Z  E 

J»-2  k»i  n»-2  I ■-7 


•V 

-  l)e  ^AWj^  .  (2.27b) 


The  respective  variances  are  respectively 
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2  8 


JOm-Z  im-7 


•[v^(t;£,,.^)  +  V^(Si£«,.i)] 


and 


2  a 


<  i!(t)  >  ■  i  i  Aw^(5A  0‘ 


ma-2  l=-7 


L  J 


{iy5;E„,Vi)  ■nys;C„,''i)]}  ' 


and  the  covariance  between  the  two  variables  is 

2  8 

<  li(S)V^(8)  >  -  Y  I  Aw^(5A  C) 


niB-a  i=-7 


Since  the  kinetic  power  density  is  a  function  of  the  cathode  temperature 
alone,  <  ij^(O)  Vj^(O)  >  should  be  equal  to  <  ij^(5j>0)  Vj^(^j,0)  >  at  the 


cathode. 
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2.9»3  Noise  Parameters.  Normalized  noise  parameters 
the  self -power  spectral  density  of  the  current  fluctuatlcxis , 
the  self -power  spectral  density  of  the  voltage  fluctuations, 
the  real  part  of  the  cross-power  spectral  density  of  the  current  cmd 
the  voltage  fluctuations,  A(Cy5)>  imaginary  part  of  the  cross-power 
spectral  density  and  S(5j>6)  for  each  Individiaal  layer  of  the  electron 
beam  are  defined  according  to  Haus's^®  notations. 


-  <  l2(Cj,S)>/2e|l(Cj)|AC  , 

=  t<  v2(5j,6)>/2(kT^)2]e|l(Cj)lAC  , 

n(5j,0  -  Re[<  l^(Cj,t)  V^(Cj,6)>/2Krg]  , 

and 

-  A®(Cy6)  .  (2.51) 


where  S  and  H  are  invariant  quantities  in  one-dimensional  single-velocity 
flow  perturbation  theory.  As  will  be  pointed  out  In  the  Monte  Carlo 
method  of  noise  transport  calculations,  the  noise  parameters  are  the 
Fourier  transforms  of  the  auto-correlation  functions  and  the  cross- 
correlation  functions  of  the  a-c  current  and  the  equivalent  a-c  voltage. 

The  minimum  noise  figure  of  a  longltudlnal-beam  microwave  amplifier 
with  one  of  the  m  layers  of  the  electron  beam  Is  again,  according  to 
Haus's  theory: 
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-"(tj'Sln)]  (2-») 

for  a  fixed  frequency,  where 

T  is  the  ambient  temperature  and 

is  the  distance  between  the  input  plane  of  the  interacting 
circuit  and  the  cathode. 

Considering  all  layers  together,  theoretically,  it  is  possible  to 
match  the  a-c  impedance  of  the  electron  beam  (sdl  layers  and  all  velocity 
classes)  simultaneously  provided  the  lengths  of  the  transformation  region 
and  the  drift  regions  are  unlimited. 

Noise  parameters  similar  to  those  for  one  layer  beams  are  written 
in  terms  of  the  total  current  fluctuations  and  the  kinetic  potential. 

t(i)  -  <  lf(|)>/2e|l^|(5AO  , 

®(6)  -  [<Vf(5)>/2(kT^)2]  e|l^l(5AO  , 

n(|)  ^  Re[<ij^(6)V^(6)>J /2kT^  , 

A(|)  -  I^[<ij,(l)Vj(|)>]/2kTj, 

and 

3^(6)  -  ♦(5)t(l)  -  A^(6)  . 

The  minimum  noise  figure  is  given  by  the  expression 

^  *  T  -  ><«>] 

for  a  fixed  operating  frequency. 


(2.35) 


(2.34) 


The  mlnlnum  noise  figure  viU  cease  to  be  a  function  of  distance 
as  sow  as  the  average  velocity  of  the  electrons  becomes  very  large 
comi>ared  to  the  equivalent  velocity  f luctuatlcsis .  From  then  on,  the 
small-signal  single -velocity  theory  will  take  over  and  no  further  reduc¬ 
tion  in  electron  beam  noise  content  may  be  achieved. 

2.10  Choice  of  Velocity  Classes  and  Integration  Steps 

A  compromise  must  be  made  between  two  contreuiictory  factors 
governing  the  nianber  of  velocity  classes  selected  and  the  sizes  of  the 
numerical  Integration  steps.  The  available  computing  time  limits  the 
maximum  number  of  discrete  velocity  classes  in  the  theoretical  beam 
model;  computing  time  varies  as  the  square  of  the  number  of  classes  and 
Inversely  with  the  size  of  integration  steps.  On  the  other  hand,  the 
small-signal  assumption  poses  a  lower  limit  on  the  number  of  velocity 
classes  employed. 

2.10.1  Choice  of  Velocity  Classes.  For  the  two-dimensional  diode, 
1^  velocity  classes  are  employed  in  the  noise  transport  calculations. 

In  Fig.  2.k  and  Table  2.1,  the  first  and  third  moments  of  each 
velocity  class  in  the  discrete  beam  model  described  in  Section  2.7  are 
identically  set  to  the  corresponding  elements  of  the  actual  beam  with  a 
half -Maxwellian  velocity  distribution  at  the  cathode. 

The  slowest  velocity  class  in  the  transmitted  flow  is  character¬ 
ized  by  B  0.04961  and  the  slowest  initial  velocity  class  has  a  w  of 


approximately  -1.^ 
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Table  2.1 

The  First  and  Third  Mcnents  of  the  Velocity  Classes 
In  the  Discrete  Bean  Model 


Velocity 

Class 


K 

1 

0.100 

2 

0.150 

3 

0.180 

4 

0.220 

5 

0.250 

6 

0.400 

7 

0.700 

8 

2.000 

-1 

0.150 

-2 

0.350 

-3 

0.550 

-4 

0.250 

-5 

0.100 

-6 

0.073 

-7  Layer 
II  &  IV 

0.278 

Leyer  III 

0.373 

Mean  Excess 
Energy  Parameter 

_ _ 

0.0496105 

0.174360 

0.338663 

0.538010 

0.7723636 

1.093326 

1.629715 

0.693112 

-0.0757257 

-0.3299202 

-0.6800619 

-0.9777244 

-1.150256 

-1.317179 

-1.396850 

-1.465326 


First  Velocity 

-V 

k 

Moment  e 

0.09516 

0.12604 

0.12829 

0.12846 

0.11548 

0.13404 

0.13719 

0.13554 

0.1618 

0.4868 

0.6909 

0.6646 

0.5159 

0.255 

1.1238 

1.6147 


Third  Velocity 
Moment  v^^e  £ai^ 

0.0047209 

0.0219763 

0.0434471 

0.0691128 

0.0891925 

0.1465494 

0.223580 

0.3644858 

-0.0122524 

-0.16060515 

-0.4698547 

-0.649795 

-0.365366 

-0.358931 

-1.56978 

-2.36606 


2.10.2  Integration  Steps  for  the  A-c  Density  Functions. 

2.10.2a  Computation  Method.  The  set  of  linear  small- 
signal  differential  equations  Is  solved  by  the  step-by-step  Integration 
method.  Variations  of  the  a-c  density  functions  vlth  respect  to  distance 
are  governed  by  the  following  fonmila: 
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Pi(C,53^,,w) 

where 


3Pi{C»6g,W) 

Ti 


A6, 


dp 


9 


(2.35) 


i 


S+l 


and 


dPj^(5,lg,w) 

dp 


dp^(C.6g^l,W)  _  ^Pi(C>6g.w) 

_ 


This  Is  actually  the  Taylor's  series  approximation  taken  to  second -order 

terms.  In  order  to  avoid  large  computation  errors,  Integration  steps 

AS  are  chosen  to  be  extremely  small  with  their  magnitudes  being  deter- 
s 

mined  by  the  increment  of  transit  angles  at  the  frequency  of  interest. 

2.10.2b  Integration  Steps  in  Errors  introduced  by  the 
discrete  step-by-step  integration  method  may  be  minimized  if  the  transit 
angles  of  the  slow  electrons  are  limited  to  only  a  small  fraction  of  a 
cycle  for  any  step  taken.  The  transit  angle  is  defined  as 


9  ^  g  /  ■  . 

J  Vw  +  ti(s) 

and 


(2.36) 


A9 


gAS 

+  T}(S) 


A| 


(2.37) 


For  the  transmitted  flows,  w  ■  0.02.  If  A9  is  chosen  to  be  less  than 
0.20  radian. 
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+0.002J  .  (2.38) 

For  the  reflected  flows,  w  -1-  ij  is  chosen  to  be  larger  than  0.02  at  all 
times  so  that  the  problem  of  infinite  phase  shift  at  the  points  of 
reflection  may  be  by-passed.  It  can  be  shown  mathematically  ( see 
Appendix  A.VIl)  that  the  phase  shift  would  not  really  increase  without 
bound  when  (w+tj)  approaches  zero  as  a  limit.  The  singularities  encoxm- 
tered  here  are  believed  to  be  the  result  of  the  discrete  beam  analysis. 
Physically,  not  all  the  electrons  in  a  velocity  class  change  their 
direction  of  propagation  simultaneously  and  those  reflected  electrons 
actually  spend  no  time  (or  undergo  no  phase  change)  at  the  point  of 
reflection. 


2.11  High-Frequency  and  Low-Frequency  Limits 

In  the  two-dimensional  diode,  the  noise-transport  analysis  is  valid 
as  long  as  both  the  distance  between  the  discrete  beam  layers  and  the 
spacing  between  the  electrodes  are  small  compared  to  a  free  space  wave¬ 
length  corresponding  to  the  highest  frequency  of  interest.  On  the  other 
hand,  at  low  frequencies,  the  limitation  described  by  Watkins^®,  Slegman 
and  Hsieh  for  one-dimensional  models  with  an  open-circuit  assumption 
fails  to  appear  in  the  present  formulation.  According  to  Eq.  k.9,  the 
z-directed  a-c  electric  field  in  a  two-dimensional  diode  no  longer  varies 
inversely  as  frequency  and  it  stays  in  bo\jnd  when  the  operating  frequency 
drops  to  zero. 


Ejy,z) 


2Jo)€„ 


OD  CQ 

J  J'  p’z'  Fj^(y’,z,z')« 


•00  •CP 


dy'dz'  ,  (a.9) 


where 


Therefore, 


•00  >00 

(2.39) 

which  Is  Independent  of  the  radian  freqxiency  (d  for  y  >  y ' .  The  a-c 
electric  field  evaluated  this  way  depends  on  the  a-c  space  charge  In  a 
way  quite  similar  to  the  dependence  of  the  d-c  electric  field  on  the  d-c 
space  charge. 


CHAFFER  111.  MOHTE  CARLO  FOBMJUTION 


$.1  Introduction 

The  Mc»ite  Carlo  method  Is  a  numerical  method  of  approximation 
with  which  Investigations  In  stochastic  models  of  physical  or  mathe¬ 
matical  processes  may  he  carried  out  hy  direct  simulation.  The  method 
hy  itself  Is  far  from  new,  nor  Is  there  any  novelty  In  the  theory  of 
stochastic  processes.  Yet,  the  appllcatloa  of  the  Mmte  Csurlo  techniques 
did  not  gain  popularity  until  high-speed  electronic  computers  with  a 
large  memory  capacity  became  available. 

For  some  equations  arising  In  a  nonprobablllstlc  context,  neither 
a  closed  form  nor  a  numerical  solution  Is  easily  obtainable  althou^  thert 
may  exist  a  stochastic  process  with  a  proper  distribution  which  satisfies 
the  equation.  Also  It  may  be  more  efficient  to  construct  such  a  process 
and  compute  statistically  than  to  attempt  using  the  more  conventional 
mathematical  tools.  In  the  noise  analyses  where  the  physical  process 
involved  Is  stochastic  In  nature,  the  Monte  Carlo  method  Is  generally 
found  adaptable.  Because  of  the  low  average  electron  velocity  near  the 
potential  minimum  of  a  thermionic  diode  operated  under  space-charge- 
llmlted  conditions,  the  velocity  spread  causes  nonlinearity  In  the 
transport  equation  governing  the  macroscopic  propagation  phenosiena.  It 
Is  rather  doubtful  that  the  llnesu*  smstU-slgnal  assumption  made  In  most 
other  methods  of  analysis  (e.g. ,  density  function)  for  a  practical  model 
Is  valid.  The  Monte  Carlo  techniques  employed  here  In  solving  the 
multivelocity  transport  problem  would  not  only  by-peuss  the  difficulties 
which  arise  from  nonlinearity,  but  would  cOao  preserve  the  psurtlcle 
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nature  of  the  electrons  and  have  the  slngulcurlty  factors  ellolnated 
from  the  mathematical  formulatlan. 

Simple  as  It  may  sound,  the  Monte  Carlo  calculations  are  usually 
extremely  lengthy  If  a  hl£^  degree  of  resolution  Is  demanded.  Sufficient 
samples  must  be  taken  so  as  to  make  the  statistical  approach  feasible. 

One  of  the  prereq\il8ltes  for  the  application  of  the  method  Is  that  edl 
physical  or  mathematical  stochastic  processes  satisfy  the  ergodlc  theorem 
In  which  the  time  average  of  a  statloneury  random  process  Is  equal  with  a 
probability  one  to  the  statistical  average.  In  this  way,  the  stochastic 
features  are  dominated  by  a  statlstlceil  approach  with  time  as  the  Inde¬ 
pendent  parameter.  Data  obtained  over  a  long  period  of  time  are 
Interpreted  the  same  as  a  large  sample  drawn  at  an  arbltreory  moment. 

The  rather  extensive  arithmetic  Involved  In  the  average  Is  easily  CGu:Tled 
out  by  modern  machine  methods. 

All  In  all,  the  Monte  Carlo  method  Is,  with  the  aid  of  high-speed 
automatic  computing  machinery,  a  means  of  returning  to  the  stochastics  In 
phase  space  to  obtain  the  thermodynamic  functions  required  rather  than 
averaging  over  all  the  phase  space  to  deed  with  the  thermodynamic  quan¬ 
tities  directly.  During  the  process  of  calctilatlons ,  all  the  work  Is 
carried  out  In  the  time  domain;  It  Is  therefore  quite  difficult  to 
observe  the  reed  mechanism  In  noise  reduction  which  Is  a  function  of 
frequency.  In  spite  of  this  disadvantage,  this  method  of  calculation  Is 
still  ccxisldered  to  be  the  most  accurate  numerical  method  for  noise  trans¬ 
port  Investigations  In  the  low-potential  region  of  a  multivelocity 
electron  beam.  Due  to  the  similarities  between  a  two-dlmensloncd  Inmersed 
flow  model  and  Its  one-dlmenslonal  counterpeurb.  It  Is  convenient  to 
describe  the  noise  transport  problem  In  one-dlmenslonal  formulatlm  first 


and  add  on  some  necessary  modiflcatlois  reqvdred  by  the  nultl-dlaenslonal 
model. 

3.2  One -Dlmens lonal  Analysis 

In  a  two-dimensional  noise  transport  problem,  the  Immersed-flow 
assumption  is  necessarily  made;  otherwise  the  model  diode  would  generalize 
into  a  three-dimensional  model  automatically  once  transverse  velocities 
of  any  kind  begin  to  Interact  with  the  axial  magnetic  field.  Mixing  of 
electrons  which  originate  at  various  parts  of  the  cathode  will  Increase 
the  complexity  of  the  problem  to  such  a  degree  that  It  Is  unlikely  to  be 
solved  In  a  reasonable  amount  of  time  even  by  the  highest  speed  conqputer 
available.  In  the  present  analysis  an  Infinite  magnetic-field  etssumption 
is  made  and  the  two-dimensional  model  may  be  approximated  by  a  set  of 
tightly  coupled  electron  sheet  beams  with  slightly  different  d-c  charac¬ 
teristics.  Each  sheet  beam  is  regarded  as  if  it  were  a  one -dlmens lonal 
electron  beam  in  the  Monte  Carlo  treatment. 

Some  assumptions  are  made  for  the  theoretical  electron  beam  model 
employed  to  simplify  the  analysis. 

1.  Transverse  motion  of  electrons  is  not  allowed,  a  direct  conse¬ 
quence  of  the  infinite  axial  magnetic  field  assumption. 

2.  All  random  processes  involved  are  ergodlc  stationery  stochastic 
processes  with  equal  time  averages  and  ensemble  averages. 

3.  Electrons  are  emitted  normal  to  the  cathode  surface  with  a 
Rayleigsh  initial  velocity  distribution. 

4.  A  potential  minimum  exists  Just  outside  of  the  cathode  surface. 

A  full  space-charge-limited  condition  prevails. 

3.  No  effect  is  produced  by  the  returning  electrcnis  on  the  emission 
characteristics  of  the  cathode. 
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6.  Direct  collisions  between  electrons  eore  prohibited.  The 
coupling  mechanism  between  electrons  comes  from  the  space-charge  fields. 

7.  There  Is  no  secondary  emission. 

8.  Accumulation  of  electrons  at  the  potential  minimum  Is  disallowed. 
Electrons  emitted  either  will  eventually  reach  the  anode  or  return  to  the 
cathode . 

9.  Relativistic  effects  are  not  considered. 


5.3  The  Open -Circuit  and  the  Short-Circuit  Assumption 

As  will  be  shown  later,  coupling  between  convection  current  density 
and  velocity  fluctuations  is  achieved  by  the  relationship  between  the  a-c 
electric  fields  and  the  a-c  current  through  Ampere’s  law,  and 


V  X  H 


J  = 


1 

c 


(3.1) 


where  €^(8E/dt)  Is  the  displacement  current  density,  1^  is  the  convec¬ 
tion  current  density,  V  is  the  a-c  voltage  across  the  diode,  Z  is  the 
r-f  impedance  of  the  diode  as  seen  by  the  electron  beam,  and  A  Is  the  area 
of  the  cathode.  The  r-f  impedance  Z  of  a  high-freqviency  diode  Is  gener¬ 
ally  considered  to  be  either  an  open  circuit  or  a  short  circuit  in  noise 
transport  problems  because  it  is  otherwise  difficult  to  determine  a 
unique  finite  equivalent  impedance  independent  of  frequency  variations. 

In  the  short-circuit  diode  theory,  the  argument  Is  based  on  the  fact  that 
the  large  inter-electrode  capacitance  represents  virtually  a  short  circuit 
at  high  frequencies.  Both  the  a-c  electric  fields  and  the  a-c  potential 
across  the  diode  approach  zero  as  a  limit  in  a  short-circuit  diode,  yet 


(3.2) 
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Is  assumed  to  hold  true  alvays.  It  Is  easy  to  see  that  no  correlatim 
between  the  current  fluctuations  and  velocity  fluctuations  can  be  obtained 
In  a  short-circuit  diode  because  of  the  absence  of  a-c  electric  fields. 

In  the  open-circuit  assumption,  coupling  between  the  electrcni  beam  and 
the  diode  circuit  Is  assumed  to  be  negligible  and  the  impedance  of  an 
electron  beam  Is  taken  to  be  an  open  circuit  at  high  frequencies  where 
the  average  transit  angle  for  the  electrons  becomes  appreciable. 


a-c 

Zk 


and 


5t 


(3.5) 


(3.1^) 


Since 


^  =  --E 
dt  m 


the  correlations  between  the  current  and  the  velocity  fluctuations  are 
a  direct  result  of  the  a-c  space-charge  effect.  In  the  Monte  Carlo 
analysis  presented  In  this  chapter  the  open-circuit  assumption  slmllau:  to 
the  one  employed  by  Dayem  and  Lambert  is  adopted  because  it  seems  to  be 
a  better  approximation  to  the  real  physical  phenomenon  than  the  short- 
circuit  assumption  made  by  Tien  and  Moshman.  In  Section  5.5;  further 
discussion  on  the  subject  of  the  a-c  impedeuice  of  two-dimensional  diodes 
will  be  presented. 


^,k  Statistical  Numbers 

3.4.1  Generation  and  Transformation  of  Random  Numbers.  There  sure 
countless  niuiber  of  ways  by  which  random  numbers  can  be  generated.  In 
the  modem  digital  computers,  rstndom  numbers  are  Internally  generated 
using  built-in  routines.  To  name  a  few  there  are  the  mid-square  method, 
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the  modulus  method,  the  Rand  method,  etc.  Because  the  end  result  of  the 
Monte  Carlo  calculation  depends  heavily  on  the  quality  of  the  random 
numbers,  it  Is  necesseory  to  make  sure  the  subroutines  are  checked  statis¬ 
tically  prior  to  application. 

By  means  of  a  random  number  transformation,  a  set  of  positive 
random  numbers  ■  1>2,3...)  may  be  obtained  having  a  desirable 

distribution  function  f(v),  where  ^v)  denotes  the  probability  of  finding 
a  number  between  v  and  v4dv  and  f(v)  Is  called  the  distribution  function 
of  V.  If  an  integral  distribution  function  P(v)  Is  defined  such  that 

V 

F(v)  »  J  f(v)  dv  ,  (5.5) 

•OO 

F(v)  must  fall  between  zero  and  one  and  it  must  be  a  single-valued  memo- 
tone  increasing  function  of  v.  In  other  words,  there  exists  a  one-to-one 
correspondence  between  v  and  F(v).  For  a  set  of  equally  probable  random 
numbers  R^(l  =  l,2,3-.0  generated  so  that 

0  <  <  1  (5.6) 

a  transformation  may  be  accomplished  by  giving  this  set  of  random  numbers 
a  one-to-one  correspondence  with  F(vj^), 

\  •  (5.7) 

The  set  of  variables  v^  corresponding  to  F(v^)  thus  obtained  would  then 
have  the  characteristics  of  random  numbers  governed  by  a  distribution 
function  f(v). 

3.4.2  Velocity  Distribution  of  Thermionic  Emission.  With  the 
Infinite  axial  magnetic  field  assumption,  all  transverse  velocities  may 
be  neglected  and  the  Initial  velocity  of  electrons  emitted  normal  to  the 
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surface  of  the  cathode  follows  approxlnately  a  Rayleigh  distribution 
provided  that  the  temperature  of  the  emitting  surface  Is  reasmably  hlc^ 
compared  to  the  equivalent  temperature  of  Its  work  function.  In  Fig. 

3.1  the  probability  density  function  that  an  electron  be  emitted  with 
velocity  between  v  and  v+dv  is  plotted  against  Its  initial  velocity  v. 
This  function  Is 


f(v)  dv  -  ^  exp  ^  ^  dv  . 


(3.8) 


c  '  c 

The  Integrated  distribution  function  is  shewn  In  Fig.  3-2  as, 

V 


F(v) 


/«v) 


dv 


1  -  exp 


[-1^]  • 


(3.9) 


If  R  Is  a  random  number  between  zero  and  one,  It  is  easy  to  see  that  1-R 
is  also  a  random  number  In  the  same  range.  As  long  as  R  is  uniformly 
distributed,  a  set  of  random  numbers  R^(i  ■  1,2,...)  may  be  transformed 
in  the  following  manner. 


1-R. 


or 


R(v^) 


1  -  exp 


(-4) 


(3.10) 


(3.11a) 


;ZED  miTIAL  VELCKTETY  DISTRIBUTIC®  AT  THE  CATHCHJE  (RAYIEIGH  Dll 
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(3.11b) 


where  is  the  initial  velocity  of  the  1th  electron  and  the  velocity 
distribution  is  governed  by  a  probability  density  function  f(v).  In 
emission  processes,  a  random  number  must  be  generated  corresponding  to 
the  initial  velocity  of  every  electron  emitted. 

3. 4-. 5  liiml ssion  Probability.  The  Instantaneous  rate  of  emlssl(»i 
during  a  very  small  time  interval  At  follows  Poisson's  distribution  as 
shown  in  Fig.  3.3,  because  of  the  postulate  that  the  emission  process  of 
thermionic  cathodes  is  stationary.  From  a  unit  area  on  the  surface  of 
the  cathode,  if  the  average  number  of  electrons  emitted  during  a  time 
Interval  At  is  n,  the  probability  that  s  electrons  will  be  emitted  during 
a  time  interval  between  t  and  t+At  is 

f(8)  =  •  (5.12) 

Since  s  assumes  only  discrete  values,  the  integration  in  Eq.  3*3  must  be 
replaced  by  a  series  summation,  such  that 


F(s)  « 


y  f(s)  • 


s 

The  plot  of  F(s)  versus  s  in  integer  values  is  shown  in  Fig.  3.4. 
a  numerical  integration  is  carried  out, 


(3.13) 


After 


F(s)  -  I  f(0)  +  f(l)  +  f(2)  +  ...  +  f(s-l)  +  i  f(s)  .  (3.14) 


It  is  convenient  to  designate  the  number  of  electrons  emitted  per  time 


Interval  At  to  be  s  if 
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3.5  PROBABILITr  DEaUSiry  of  EIECTROK  EMISSICW  IER  UKET  time  UJTEBVAL.  (n  =  7.5509  EIECTROiS) 
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where  R  Is  a  raxidom  number. 

^  8 

F(s  -  1/2)  and  F(s  +  I/2)  may  be  obtained  by  means  of  graphical 
method  from  the  F(8)  versus  s  curve. 

3. *<■.4  Emission  Time  of  Electrons .  Because  of  the  discrete  nature 
of  the  digital  computer  techniques,  the  exact  emission  time  of  all  elec¬ 
trons  during  a  time  Interval  At  may  be  simulated  by  a  random  variable  with 
a  uniform  distribution  function.  However,  for  computation  purposes  It 
would  be  much  more  convenient  If  all  the  electrons  are  considered  to  be 
emitted  at  the  end  of  each  time  Interval  at  a  small  distance  from  the 
surface  of  the  cathode,  where  1  »  1,2,...  corresponding  to  the  1th  number 
of  electrons  emitted  during  the  immediately  preceding  time  Interval.  The 
small  displacement  z^  Is  actually  the  distance  traveled  by  the  1th 
electron  during  the  fraction  of  At  between  the  time  of  Its  emission  and 
the  end  of  the  time  interval. 


1  ^ 

1  +Tf 


cathode 

dz 


(3.16) 


where  v^^  is  the  initial  velocity  of  the  ith  electron  at  the  cathode 
plane,  R^  Is  the  random  number  corresponding  to  the  emission  time  of  the 
ith  electron,  Is  the  electric  field  as  seen  by  electrons  at  the 

cathode  surface.  No  random  number  transformation  Is  needed  for  the  set 
of  random  numbers  corresponding  to  the  emission  time  of  electrons. 

The  velocity  of  the  1th  electron  at  Is  then 


v!  ■  V.  -  h  E 


'cathode 


dE 


'cathode 

dz 


R,At 


(3.17) 
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3.5  Electric  Flelde  In  a  TVo-Dlmenalaial  Diode 

Inside  the  diode,  electric  fields  experienced  hy  each  electron  due 
to  the  presence  of  the  space  charge  are  to  be  evaluated  at  the  beginning 
of  every  unit  time  Interval  and  the  d-c  electric  fields  are  assumed  to 
be  unchanged  throughout  At.  Because  of  the  divergence  conditions 
encountered  In  an  attempt  to  evaluate  the  electric  field  Inside  a  tvo- 
dlmenslonal  diode  through  tem-vlse  differentiation  of  the  Fourier  series 
representation  of  the  potential  function,  the  Image  method  Is  Instead 
employed  to  solve  Poisson's  equation.  Details  of  both  the  Fourier  series 
method  and  the  Image  method  are  presented  In  Appendix  B.I. 

As  a  result  of  the  infinite  axial  magnetic  field  assumption,  the 
only  current  flowing  perpendicular  to  the  two-dimensional  ribbon  beam  Is 
the  displacement  current  Induced  by  the  variations  of  the  surface  space- 
charge  density  of  the  z -directed  sheet  beams.  According  to  the  equation 
of  continuity  which  requires  the  net  current  flowing  out  of  any  volume  In 
space  to  be  zero,  the  existence  of  a  net  current  In  the  transverse  direc¬ 
tion  must  be  accompanied  by  a  ncxizero  algebraic  sum  for  the  longitudinal 
currents  from  the  same  volume.  In  other  words,  the  Impedance  of  a  truly 
two-dimensional  diode  will  remain  finite. 

In  the  present  analysis,  since  It  Is  difficult  to  determine  how 
rapidly  the  a-c  electric  fields  will  decay  In  the  transverse  dlrectlcm, 
some  additional  assumption  on  their  boundary  condltlcxis  are  made. 

Consider  a  ribbon  beam  whose  thickness  Is  very  sms^JL  compared  to  a 
free-space  wavelength  corresponding  to  the  highest  frequency  of  Interest, 
where  the  Induced  a-c  electric  fields  may  be  assumed  uniform  across  the 
electron  beam.  Outside  of  the  beam  boundcurles,  where  the  large 
Inter-electrode  capacitance  of  the  diode  behaves  like  a  short  circuit 
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at  hl£^  frequencies,  no  z -directed  a-c  electric  fields  or  dlsplaceaent 
current  Is  allowed.  Now  If  the  electron  heam  Is  divided  axially  Into 
otany  sections,  the  transverse  displacement  current  flowing  traa  any 
section  may  be  made  arbitrarily  small  coopared  to  the  respective 
z-dlrected  component  by  reducing  the  ratio  between  the  axial  dimension 
of  the  section  and  the  thickness  of  the  electron  besa.  Under  such  condi¬ 
tions,  a  close  approximation  of  the  axial  a-c  electric  field  nay  be 
obtained  by  means  of  the  open-circuit  assumption.  It  follows  from  the 
added  boundary  conditions  that  one-dlmenslcxial  a-c  electric  fields  are 
actually  employed  In  the  Monte  Carlo  noise  transport  analysis  In  a  diode 
with  a  two-dimensional  d-c  potential  distribution.  A  comparison  of  the 
results  from  the  present  analysis  with  those  from  the  one-dlmenslonal 
diodes  will  therefore  reveal  the  effects  aa  the  noise  transport  charac¬ 
teristics  by  the  d-c  space  charge. 

3.5.1  Grouping  of  Electrons.  The  effects  of  space  charge  on  d-c 
electric  fields  and  potential  distribution  may  be  obtained  by  equations 
derived  in  Appendix  B.I.  Unfortunately,  the  large  amount  of  time  needed 
to  compute  the  d-c  electric  fields  alone  makes  the  Monte  Ceurlo  analysis 
highly  Impractical  even  for  the  highest  speed  computer  available.  One 
way  of  reducing  the  length  of  the  problem  is  to  divide  the  diode  Into 
some  compartments,  each  containing  a  few  electrons  as  Illustrated  In 
Appendix  B.II.  These  small  groups  of  electrons  are  In  turn  replaced  by 
x-dlrected  line  charges  equal  to  the  sum  of  all  the  electric  chsurges  In 
the  respective  compeurtments  and  are  located  at  their  respective  effective 
centroids.  After  the  potential  distribution  and  the  electric  field  are 
computed  for  all  the  boundary  lines  between  cosqpartments,  the  d-c 
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Bpece-charge  field  axiyvhere  Inside  the  diode  may  be  evaluated  by  linear 
Interpolation  methods. 

Further  slnqpllflcatlon  of  the  coe^tatlon  procedure  nay  be  achieved 
by  fixing  the  centroids  at  their  corresponding  coeipartments  regardless  of 
the  actual  space-charge  configurations.  Techniques  of  linear  transform 
are  applicable  once  the  centroids  are  considered  fixed  In  space,  because 
the  resulting  electric  field  Is  a  linear  coadDlnatlon  of  all  the  electric 
fields  due  to  the  presence  of  space  charge  at  the  various  centroids. 

3.5.2  Potential  Distributions .  Due  to  the  absence  of  space  charge 
outside  of  a  tvo-dlnenslonal  ribbon  bean,  the  potential  distribution 
Inside  a  tvo-dlnenslonal  space-charge-llmlted  diode  will  deviate  consid¬ 
erably  from  lansnulr's  curve  for  one-dlnenslonal  diodes  In  which  a 
uniform  space-charge  distribution  In  the  transverse  dlrectlcm  Is  assumed, 
unless  an  appropriate  applied  potential  distribution  Is  introduced  to 
replace  the  potentleJ.  depression  effect  of  the  missing  electrons.  In  a 
synthetic  model,  the  desirable  potential  distribution  may  be  realized  by 
properly  placing  a  fixed  fictitious  space  charge  between  the  electrodes. 

An  applied  potential  distribution  corresponding  to  Child's  lav  Is  assumed 
in  the  two-dimensional  diode  employed  in  the  present  analysis  to  account 
for  the  effect  of  the  missing  space  charge.  As  a  result  of  the  Image 
method,  the  potential  at  a  point  between  compartments  (i,4)  and  (l,J+l) 

Is 

vCi,j)  -  Y  I ''h-Ji-.”)  +  ,  (3.18) 

mal  n>nl 

where  V(l,J;m,n)  Is  the  potential  at  a  point  due  to  the  electric 

charge  at  the  centroid  of  the  (m,n)  compartment,  the  last  tens  Is  the 
applied  potential  at  z^  according  to  Child's  law,  a  Is  the  distance 
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between  the  cathode  and  the  anode  and  is  the  applied  potential  at 

the  anode. 
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^  Is  the  charge  per  unit  length  In  the  (m,n)  CQanpea*tinent.  For  Zj  +  z^ 
<  0.2a 


Pm  n  r  r  1  5Z  .Z  'I 

V(l,Jjm,n)  m  ^  4in  ^  J.  ^  f  •  (5-20) 
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The  time  average  of  the  discrete  potential  functions  thus  obtained 
will  be  put  into  series  form  for  each  layer  of  the  electron  beam  so  that 
a  continuous  potential  distribution  will  be  available  for  the  noise  trans¬ 
port  analysis  In  Chapter  II  where  the  density  function  method  Is  employed. 

3.^.3  D-c  Electric  Fields .  The  z-dlrected  d-c  electric  field 
experienced  by  an  electron  within  the  (l,j)  compartment  Is  obtained  by 
term-by-term  differentiation  of  the  potentisa  function.  Electric  fields 
at  points  between  compartments  are  first  evaluated,  then,  linear  inter¬ 
polation  in-betwem  will  be  made  for  each  Individual  electron  Inside  the 
compartments. 

The  z-dlrected  electric  field  at  a  point  between  (i,j)  and  (l,J+l) 


may  be  written  as 


.6k- 


-  2  v^/® 

E(i,j)  -  ^  ^  ,  (3.21) 

BbI  nal 

where  E(l)J;m,n)  Is  the  electric  field  at  dvie  to  the  presence  of 

a  charge  ^  at  the  centroid  of  compartment  (m,n) ,  and  -0<V^^^^/3a) 
(zVa)^/®  Is  the  d-c  applied  electric  field  corresponding  to  the  last 

J 

term  In  Eq.  3*18« 

Again  from  Appendix  B.I 


E( 


(Zj.2 


^ _  1  -  (z  +Z  )  [  _ i _ 

-  (2a)2  J  J  n  L  (zj+Zq)®  +  (Vi-yJ 


J  n' 


2  -  (2a)2  ]  2a  [ 


(2a+z.+z„)t2a  -  (z  -z^)] 


x-j 

(z +z  )2  -  (2a)2  J  ■‘='®  L  (2a+z,-z„)[2a  -  (z,+z. 


J  n' 


j«„)i  Jr 


or.  If  z,+z  «  0.2a,  the  equation  above  may  be  reduced  to 
0  n 


E(l,Ji»,n)  -  a€°  { 


Z,+  Z 

J  n 


(5.22) 


(5.23) 

For  an  electron  in  compartment  (i,j)  such  that  z  ■  Zj,  where  (y^^iZj)  are 
the  coordinates  of  the  centroid  of  compartment  (l,j),  the  d-c  electric 
field  may  be  evaluated  by  linear  Interpolation. 


z-z. 


E(y,z)  .  E(i,j)  +  [E(1,J+1)  -  E(i,j)]  - - i 


^J+1)  ■ 


or,  if  z  • 


(3.24) 
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E(y,z)  -  E(i,j\+  [E(i,J-l)  -  E(1,J)] -  .  (3.95) 

’“(J-I) 

3»5.^  Initial  Potential  PlBtribntloo  In  the  Diode.  A  one- 
dlnenslonal  potential  distribution  as  Illustrated  In  Fig.  3.5  Is  assuned 
In  the  tvo-dlmenslonal  diode  for  the  first  tvo>hundred  tine  intervals  so 
that  the  diode  nay  be  properly  filled  up  with  electrons  before  fomulas 
listed  In  the  previous  sections  are  employed  to  compute  the  electric 
fields.  Approximation  of  Langmuir's  one-dlmenslonal  space-charge-llnlted 
diode  potential  distribution  curve  Is  made  for  various  sections  of  the 
diode.  In  dimensionless  form 

1.  0  <  6  <  1.9. 

Tj  -  0.25(1.9-1)^  +  0.01^5(1.9-1)*  +0.0045(1.9-1)® 

and 

-  -  0.50(1.9-6)  +  0.18(1.9-6)®  +  0.027(1.9-6)®  . 

2.  1.9  <  6  <  6. 

T,  =  [0.424(6-1.9)]^-'^^^  , 

-  =  -  1.735  [0.424(6-1.9)]°*^^^  X  0.424  . 

3.  6  <  6  <  12. 

n  =  [0.4405(6-1.9)]^'^^®  , 

.  M  .  .  1.618  X  0.4405  [0.4405(6-1.9)]°'®^®  . 


(I0“*  METER) 
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12  <  e  <  if2. 

n  -  [0.5055(1-1.9)]^'*^^  , 

-  .  -  1.483  X  0.5035  [0.5035(6-1.9)]°**^^  . 

5.  42  <  6  <  200. 

n  -  [0.5925(1-1.9)]^'**°®  , 

-  |a  =  -  1.408  X  0.5925  [0.5925(6-1.9)]°'**°®  , 

where  q  and  6  are  the  potential  and  distance  parameters  defined  In 
Chapter  II,  -  (dtj/dj)  Is  proportional  to  the  d-c  electric  field. 

The  potential  minimum  plane  for  a  one-dlmenslonal  diode  with  the 
same  emission  characteristic  as  the  model  employed  here  Is  located  at 
approximately  «  1.9.  a  two-dimensional  model  the  potential  minimum 
of  the  center  layer  of  the  electron  beam  Is  expected  to  be  approximately 
the  same  distance  away  from  the  cathode  surface. 

3.5.5  A-c  Electric  Fields .  Since  the  diode  Is  assumed  to  be  an 
open  circuit  at  high  frequencies,  the  total  r-f  current  flowing  across 
any  plane  parallel  to  the  surface  of  the  cathode  Is  zero.  Or  In  terms  of 
current  densities, 

SE(y. ,z  )  ^ 

«o - -*5  L  +Jd.c  '  °  '  <5.26) 

1-1 

where  [pv](yj^,Zj)  Is  the  convection  current  density  carried  by  the  elec¬ 
trons  In  compartment  (l,j),  Is  the  average  current  density  of  the 
electrraa  beam,  dE/8t  Is  the  time  rate  of  change  of  the  electric  field 
experienced  by  all  electrons  In  canqpartments  (l,j),  where  1  >  1,2,... 5. 


-68- 


^  -  (3-27) 

I  J 

where  i  Is  the  number  of  electrons  In  compartment  ds/dt  Is  computed 

at  the  beginning  of  each  time  Interval  At  and  Is  assumed  to  be  a  constant 
throughout . 

The  relation  given  between  the  a-c  electric  fields  and  the  current 
fluctuations  In  this  two-dimensional  model  suggests  that  the  a-c  electric 
fields  drop  off  abruptly  beyond  the  boundaries  of  the  electron  beam  In 
the  y-dlrectlon.  In  a  real  diode,  the  field  strength  Is  expected  to  decay 
exponentially  Instead. 

3.6  Numerical  Integration 

Electrons  are  traced  In  phase  space  by  means  of  numerical  Integra¬ 
tion  techniques.  The  incremental  changes  In  velocity  and  distance  are 
governed  by  the  electric  field  throu^  the  use  of  Newton's  force  laws 
to  the  first-order  approximation  In  both  time  and  space. 


Interval  At 
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E.  Is  the  total  d-c  electric  field  at  the  point  where  the  electron 

<1*C 

Is  located.  E^_g  ■  ^applied  *  ^space  charge* 
dE/dt  Is  the  rate  of  chemge  of  the  a>c  electric  field  obtained  by  the 
open-circuit  diode  assumption.  It  Is  assumed  to  be  constant 
over  a  small  time  Interval  At,  l.e.,  fAt  Is  much  smaller  than 
unity  If  f  Is  the  highest  frequency  of  interest. 
bE/hz  Is  the  slope  of  the  d-c  electric  field  at  the  Initial  position 
of  the  electron  at  the  beginning  of  the  time  IntervcLI. 

Az  Is  the  Incremental  change  of  velocity  during  At  for  an  electron 
traveling  with  a  velocity  z  and  located  at  z  at  the  beginning 
of  the  time  Interval. 

Az  Is  the  incremental  change  of  distance  occurring  during  At. 
Numerical  Integrations  are  carried  out  for  every  electron  between 
the  two  electrodes  for  every  time  Interval  At  until  they  cross  the  anode 
plane  or  return  to  the  cathode.  Positions  and  velocities  of  electrons  are 
stored  temporarily  and  then  updated  at  the  end  of  each  Iteration. 

The  second-order  term  l/6  (d^E/St®)(At)®  is  assumed  to  be  negli¬ 
gible  compared  to  the  first-order  effect  of  the  a-c  electric  field  in  the 
Monte  Carlo  analysis. 

3.7  Computation  of  Noise  Parameters 

Noise  parameters  are  evaluated  at  several  planes  beyond  the  poten¬ 
tial  minimum  for  a  set  of  discrete  frequencies.  According  to  the  noise 
theory  derived  by  Haus  and  Roblnscai,  the  minimum  noise  figure  achievable 
by  a  longltudlnal-beam  amplifier  depends  solely  on  these  parameters,  which 
are  the  Fourier  transforms  of  the  correlation  functions  of  the  ccmvectlon 
current  density  and  the  voltage  fluctuations. 
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3.7.1  Cvtrrent  and  Velocity  Fluctuations.  The  current  and  the 

velocity  fluctuations  are  computed  for  planes  located  at  z  ■  z^^,  z^,... 

z  for  each  tine  Interval  between  t.  and  t.4At.  The  current  fluctuations 
n  1  i 

at  (y  )  during  the  ith  tine  interval  nay  be  expressed  as 


where  the  subscript  "1"  denotes  an  r-f  quantity,  and  "o"  denotes  a  d-c 

quantity,  number  of  electrons  crossing  z^^  at  layer  y^^^ 

between  t^  and  t.+At,  and  N  (y  ,z  )  is  the  average  number  of  electrons 
1  1  o  n  n 

crossing  plane  z  at  the  y  th  layer  for  an  Interval  At.  Similarly,  for 
n  m 

the  electron  beam  as  a  whole,  the  current  fluctuation  at  z^  during  the 
1th  time  Interval  is 

5 

mal 

where 

5 

W  -  5  I  • 

mol 

As  to  velocity  fluctuations,  the  problem  is  much  more  complicated 
if  the  retarding  field  region  of  the  diode  Is  considered.  Because  of  the 
reflected  flows,  it  will  be  necessary  to  assign  each  electron  a  name  tag 
and  follow  through  to  see  whether  it  will  eventually  pass  the  potential 
minimum  or  not.  However,  it  is  less  difficult  to  determine  velocity 
fluctuations  at  planes  beyond  the  potential  minimum.  Chu's  kinetic 
potential  may  be  Introduced  in  place  of  velocity  fluctuations.  By 


deflnitlcm 


-TI¬ 


MS’ 

♦ 


"o'V'n) 


I  -  ^V'n>]  ' 

i 

(3.33) 


where  I  is  the  number  of  electrons  crossing  plane  at  layer  y^^  between 

t.  and  t.+At,  v,(y  ,z  ,t.  )  Is  the  velocity  of  the  1th  electron  at  z  as 

shown  in  Appendix  B.III,  and  is  the  mean  square  velocity  of 

electrons  at  the  z  plane  of  the  layer  y  . 

n  m 


v®(y  ,z  ) 
o'^m'  n' 


and  for  the  whole  beam. 


an  V,  ; - --i -  £  Y  ' 

(JOS*) 

where 


v^(  z  ) 
o  n 


Y.  I  ' 


III  i<V'n'\c>  • 


tj^  m«i 


The  noise  power  carried  by  the  two-dimensional  electron  beam  as  a  whole 
is  expected  to  be  somewhat  higher  than  that  carried  by  Individual  layers, 
because  of  the  non-Maxwellian  velocity  distribution  in  the  multi-layer 
model. 

3.7«2  Time  Auto-Correlation  and  Time  Cross-Correia tioo  Puncti<»e . 
The  time  auto-correlation  functlcxi  is  an  indication  of  the  amount  of  mean- 
square  fluctvtatlons  as  well  as  the  time  dependence  characteristics  of  a 
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sample  functloi  of  a  random  process.  On  the  other  hand,  the  tine 
cross-correlation  functlcm  of  the  sanq^le  fvmctlons  of  two  random  pro¬ 
cesses  describes  the  dependency  of  the  two  random  processes. 

From  the  definitions  given  In  Appendix  B.IV.l,  the  time  auto¬ 
correlation  and  cross-correlation  functions  of  the  kinetic  potential  aM 
the  noise  current  density  fluctuations  are  given  by  the  following 
equations: 

(p-s^At 

*  P^  Z  '  ^5-55) 

t-0 


(p-^At 


i|f(y  ,z  ,8)  »  - 

n'  p-8 


^5.36) 


t«0 


and 


(  p-^At 

9(y^,z^,8)  -  ^  ^  '  ^5.37) 

t-0 


where  s  ie  an  Integer  less  than  or  equal  to  q  such  that  q  <  p,  and  p  Is 
the  maximum  number  of  computer  time  intervals  (total  number  of  Iterations 
carried  out) .  It  Is  postulated  that  negligible  correlation  can  be  found 
between  events  separated  by  a  time  interval  larger  than  qAt.  The  magni¬ 
tude  of  p  is  limited  mainly  by  the  available  time  and  money  for  computer 
operations. 

The  choice  of  q  depends  on  two  conflicting  considerations.  First 
of  all,  the  frequency  resolution  of  the  statistical  noise  analysis 
increases  with  q;  no  new  information  in  frequency  domain  may  be  extrsicted 
at  discrete  values  of  frequency  less  thanAf  apart,  where  Af  -  l/qAt.  On 
the  contrary,  the  confidence  Interval  of  the  statistlcsd  approximation 
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deterlorates  with  any  Increase  of  the  ratio  of  q/p.  It  is  therefore 
necessary  to  sane  kind  of  cosqproiBlse  between  the  two  llaltlaag  feetors 
keeping  In  mind  that  q  should  be  of  the  order  of  one  hundred  so  that  At 
Is  much  smaller  than  l/dt,  the  upper  frequency  limit  of  the  spectral 
analysis. 

3.7.3  Power  Spectral  Density  and  the  Noise  Parameters.  The 
application  of  the  power  spectral  densities  and  the  correlation  functions 
of  random  processes  are  usually  referred  to  as  the  generalised  harmonic 
analysis.  For  a  stationary  random  process  such  as  a  signal  or  noise,  the 
power  spectral  density  may  be  considered  m  the  power  distribution  In  the 
frequency  domain.  The  normalized  self-power  spectral  density  and  the 
normalized  cross-power  spectral  density  may  be  obtained  by  expanding  the 
auto-correlatlon  and  the  cross-correlation  functions  respectively  in 
Fourier  series  as  shown  In  Appendix  3.17.2. 

el(y  )At  r 

+  2  ^  <p(y^,z^,s)  cos  cos  (nO)j  ,  (3*38) 

s«i 


'm'  n'  ’  eHy 


0) 


+  2  cos  +  t(y„>an;<l)  cos  (xn)j  ,  (5*39) 


(5^0) 


-  J 


At 

kT. 


{  Y  -  ®2(  (x) 


s«i 


®,(y  >z  >q.)  -  s  (y  »z  >q.) 

1  am'  n' 


Bln 


(«n)}  , 


where  fl  ■  2qfAt,  a  frequency  parameter, 
f  ■  the  freqviency  variable. 

If  f  -  0,  l/2qAt,  l/qAt,  3/2  (l/qAt),  5/2  (l/qAt),  0  -  0,1,2, ...q, 
the  last  term  in  Eq.  3«^0  vanishea  identically.  *(yjjj>*Q>^)  is  the  power 
spectral  density  of  the  kinetic  potential  at  is  the 

power  spectral  density  of  the  convection  current  density  fluctuation  at 
(yjjj,z^),  e(y^,z^,n)  represents  the  cross-power  spectral  density  of  the 
kinetic  potential  and  the  convection  current  density  fluctuations  at 


(p-^At 

t-0 

( p-OAt 

0_(y  ,z  ,s)  -  y  i,(y  ,z  ,t46At)  V,(y  ,z  ,t)  . 

a'^m'  n'  '  p-s  ^  n'  I'-'m'  n' 

t«0 


Similar  expressions  may  be  obtained  for  the  power  spectral  densities  of 
the  entire  electron  beam  if  the  parameter  y^  is  dropped  frcn  Eqs.  3.35 
throu^  3.40. 


-75- 


The  noise  peoraiaeters  defined  by  Haus  and  Robinson  are  enqplayed  in 
the  noise  transport  analysis. 

n(yju,2n>o)  -  Re 

A(yj^,Zn'^^  “  ^  (3.^2) 

and 

For  one-dlnenslonal  longitudinal  beam  single  velocity  amplifiers, 
both  S(y^,z^,fl)  and  Il(y^,z^,fl)  are  independent  of  z^,  i.e.,  they  are 
invariant  quantities.  S(fl)  represents  the  total  a-c  power  carried  by  the 
fluctuations  and  11(0)  specifies  the  amount  of  correlation  between  velocity 
and  current  fluctuations.  In  general,  ll(y^,0,n)  is  assumed  to  be  zero  at 
the  cathode  plane  for  a  diode  with  uncorrelated  velocity  and  current 
fluctuations  at  the  input. 

3.7.^  Frequency  Aliasing  and  Lag  Windows . 

3.7.4a  Cutoff  Frequency  f^.  In  a  generalized  harmonic 
analysis  for  a  random  time  process  with  samples  drawn  at  eqvially  spaced 
intervals  in  time,  part  of  the  information  contained  by  the  real  con¬ 
tinuous  process  is  lost  because  of  the  assumption  of  the  invariance  of 
fluctuations  and  their  correlation  functions  within  small  unit  time 
Intervals  At.  If  the  correlation  function  R(t)  Is  an  absolutely 
Integrable  function  of  t,  the  Fovirier  treuisform  of  R(t)  exists  and  is 
called  the  power  spectral  density  Sp(f), 
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00 

Sj^(f)  .  J  R(t) 


-JOftT  . 
e  dT  . 


(5.44) 


If  the  average  value  of  R(t)  in  a  time  interval  between  t  -  (At/2)  and 
t  +  (At/2)  is  designated  as  R('r),  then 


*  At 


s  /  »(t) 

.  At 

t  -  -TT 


(3.45) 


S=(f)  is  defined  as  the  estimated  spectrum  of  R(t),  such  that 
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Froo  Fig.  }.6  It  can  b«  shown  that  sin  (]CfAt)/]cfAt  Is  approfxiaately  equal 
to  one  for  small  arguiKnts.  The  estimated  spectrum  falls  off  very  rapidly 
as  fAt  approaches  unity.  At  f^  ■  estimated  spectrum 

drops  off  to  zero,  and  f^  is  called  the  cutoff  frequency  of  the  spectrum 
analysis.  Therefore  the  true  spectrum  may  be  recovered  cxily  at  frequen¬ 
cies  very  low  compared  to  the  cutoff  frequency  f^. 

3.7»‘^b  Aliasing  Effects .  As  seen  in  Fig.  3*7>  aliasing 
between  spectra  centered  at  different  frequencies  arises  from  the  discrete 
nature  of  the  parameter  t  of  the  correlation  functions.  If  s  ■  Jo),  the 
Fourier  transform  of  an  eqvially  spaced  Dirac  comb  (a  set  of  delta  func¬ 
tions  equally  spaced  in  the  time  domain)  may  be  expressed  in  the 
following  manner. 

CO  CO 

s(s)  =  y' 

•CO 

The  right-hand  side  of  Eq.  is  called  a  two-sided  Laplace  transform, 

a  generalized  case  of  the  Fourier  transform.  If  the  transforms  of  both 
R(t)  and  8(T-nAt)  exist,  then 
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R(t)  6(T-nAt)e"®^  dr 


(3.‘^7) 
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By  the  coDvolutlon  theorem, 

s(8)  •  ^  f  Sj,(s-u)  a^(u)  du  .  (5.50) 

By  the  convolution  theorem, 

where  the  first  line  Integral  Is  taken  along  the  Jtn-axls  and  a  semicircle 
of  infinite  radius  on  the  right  half  plane  to  account  for  all  the  poles 
in  the  right  half  plane.  Then,  a  second  integration  is  performed  along 
the  Ja>-axls  and  a  semicircle  of  infinite  radius  on  the  left  half  of  the 
S  plane  (Fig.  3*8) •  For  e"^^  ■  1 


u  = 


1 

At 


n  =  0,  ±1,  ±2,  ... 


(5.52) 


Therefore 


S(s) 
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Residues  of 


Sj^(  8-u) 


-uAt 

1-e 


n>-oa 


(5.53) 


The  resulting  spectrum  in  Fig.  3*9  consists  of  a  fundamental  mode 
centered  at  zero  frequency  and  repeats  itself  every  f  >  n/At.  The 
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spectrvm  observed  anywhere  will  therefore  Include  a  certain  amount  of 
undesirable  effects  from  all  the  modes  other  than  the  fundaaiental.  It 
Is  recalled  from  the  previous  section  that  the  fundamental  mode  drops 
off  rapidly  to  zero  as  f  approaches  n/At.  If  a  good  approximation  to 
the  true  spectrum  Is  desired,  the  hl^st  frequency  of  Interest  should 
be  taken  to  be  less  than  l/lOAt  where  all  aliasing  modes  are  negligible 
compared  to  the  fundamental. 

3.7.^c  Lag  Windows.  During  the  computation  of  the  noise 
parameters,  all  samples  taken  farther  than  a  time  Interval  qAt  apart 
are  assumed  to  be  completely  uncorrelated.  Power  spectral  densities 
under  such  as  assumption  may  be  written  as 


S(8) 


00  00 
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D^(T)6(T-nAt)e 
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(5.5!^) 


where  =  1  for  1t|  <  qAt, 

s  0,  otherwise, 

D^(t)  Is  called  the  lag  window  (Pig.  3»10)* 

By  the  convolution  theorem 

00 

S(s)  -  J  dv  ,  (5*55) 

•00 

where  Q^Cf)  =  2qAt(3ln  2»fqAt)  approaches  a  delta  function  aa  qAt  -*00 
and  the  introduction  of  the  finite  time  correlation  lag  window  constitutes 
an  averaging  process  with  as  the  weighting  factor  oa  the  frequency 

spectrum.  Unfortunately,  this  functlcm  has  positive  and  negative  peaks 
of  appreciable  magnitudes  ccmpared  to  the  main  peak  at  f  ■  v  and  these 
undesirable  characteristics  make  It  difficult  to  relate  S(s)  to  S(s). 


-8lt- 


FIG.  5.10  LAG  WIMDOWS  Dq(t)  =  1,  t  <  (^t;  Dq(t)  =  0,  t  >  c^t; 
D^(t)  =  1/2  (1  +  cos  rtr/^t). 
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In  this  analysis  a  new  lag  vindcw  D]^(t)  ■  l/2  [1  *  cos  XT/qAt]  Is 
selected  which  has  the  Fourier  transform 

9i(«  -  I  l«o  (f  *  sfe)  *  «<,  (' -  ife)  ]  ♦  I  %<')  • 

As  Illustrated  In  Fig.  the  side  loops  drop  off  much  more  rapidly 

than  those  of  the 


V 


s 

SMOQNIM  nVMlOadS 


CHAPTER  IV.  RESUMS  AHD  IMTERPRB'IATICWS  OF  THE  DENSin 


FUNCTION  CAICUIATIONS 


4.1  Introduction 

In  this  chapter,  the  results  of  the  noise  transport  analysis  of  a 
discrete  "beam  two-dimensional  space -charge -limited  diode  hy  means  of  the 
density  function  method  are  presented.  Correlations  between  the  trans¬ 
port  phenomena  and  the  coupling  mechanism  are  sou^t.  A  comparison  Is 
made  between  the  a-c  coupling  term  employed  In  the  two-dimensional 
density  function  calculations  and  the  one-dlmensloneLL  a-c  coupling  term 
resulting  from  the  open-circuit  assumption. 

Calculations  are  performed  for  three  discrete  values  of  the  fre¬ 
quency  parameter  a  =  0.25,  0.5  and  1.0,  corresponding  to  2  Gc,  4  Gc  and 
8  Gc  respectively.  Equivalent  noise  parameters  Y,  4>,  II,  S  and  for 
the  transmitted  flows  are  evaluated  at  ei^t  planes  located  at 
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■vrtiere  is  1he  average  distance  between  the  cathode  and  the  potential 
minima  of  beam  layers  Nos.  2  and  4,  ^  >  6  x  10*^  m.  Because  of  the 
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multlveloclty  nature  of  the  electron  130801,  especleLLly  In  the  vicinity  of 
the  potential  mlnlmmi,  the  meaning  of  the  equivalent  noise  parameters 
must  1)e  Interpreted  vlth  great  care.  According  to  the  noise  theory  ly 
Haus,  the  noise  parameters  vere  originally  defined  for  single-velocity 
one -dimensional  electron  beams .  The  noise  parameters  presented  In  this 
chapter  are  actually  the  equivalent  noise  parameters  evaluated  after  the 
electron  beam  Is  suddenly  accelerated  to  a  very  hl^  velocity.  In  the 
retarding  field  region  between  the  cathode  and  the  potential  minimm 
where  reflected  flows  exist,  the  real  i^iysical  meemlng  of  the  noise 
parameters  corresponding  to  the  transmitted  flow  is  not  clearly  estab¬ 
lished  and  neither  are  these  quantities  physically  measurable .  However, 
they  may  still  serve  to  give  some  indication  of  the  noise  transport 
characteristics  of  the  electron  beam. 

It  takes  approximately  82  minutes  of  computing  time  on  an  IBl  709O 
computer  to  complete  the  calculations  for  each  selected  value  of  the 
frequency  parameter.  No  attempt  has  been  made  to  refine  the  calculation 
for  better  accurewy. 

4.2  Noise  Parameters 

Normalized  noise  parameters  for  the  entire  beam  are  plotted  against 
the  normalized  distance  from  the  cathode  for  the  three  chosen  frequencies. 
At  the  cathode  plane, 

S(0)  =  f(0)  =  $(0)  =  1.0 
n(o)  =  a(0)  =  0 

according  to  the  initial  conditions  and  the  minimum  noise  figure  Is  6.96 
db  corresponding  to  full  thermionic -emission  shot  noise. 
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Llttle  variation  In  noise  transport  characteristics  Is  observed 
from  layer  to  layer  of  the  electron  'beam.  As  will  be  pointed  out  In 
Section  k.'i,  the  coupling  a-c  electric  field  decays  too  rapidly  In  the 
transverse  direction  to  bring  forth  any  significant  two -dimensional 
effects . 

4.2.1  Self -Power  Spectral  Density  fSPSP).  The  normalized  self- 
power  spectral  density  '¥  of  the  convection  current  density  fluctuations 
and  the  normalized  self -power  spectral  density  $  of  the  kinetic  potentlel 
are  shown  as  functions  of  normalized  distance  from  the  cathode  In  Figs. 

4.1  and  4.2  respectively.  Spatial  variations  are  observed  for  both  Y 
and  <I>  for  all  frequencies.  In  the  low  potential  region,  the  spatieJ. 
variations  of  Uf  and  ®  appear  to  be  In  quadrature,  In  agreement  with  the 
transmission-line  analog  of  electron  beams.  In  addition,  the  wavelengths 
of  the  spatial  variations  are  inversely  proportional  to  frequency.  For 
a  =  1.0,  equivalent  to  8.25  Gc,  a  full  wavelength  Is  measured  to  be 
approximately  24  x  10"®  meter.  The  velocity  v  may  be  evaluated  frcm 

V  =  f\ 

=  8.25  X  10"®  X  24  X  10"® 

=  1.9752  X  10®  m/sec. 

This  turns  out  to  be  the  average  velocity  of  electrons  in  the  vicinity 
of  the  potential  minimum  of  the  diode.  Large  amplitude  variations  of  f 
are  observed;  in  fact,  the  convection  current  density  fluctuations  tend 
to  dominate  the  noise  characteristics  of  the  electron  beam  as  evidenced 
by  the  following  discussion.  !!3ie  basic  noise  reduction  mechanism,  namely, 
the  fluctuating  potential  mlnlmm  effect.  Is  not  founa.  The  spatial 
variation  of  the  fluctuations  seems  to  have  come  frcm  some  beatlngHrave 


effects  Instead 
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4.2.2  Holse  Parameters  S  and  H.  In  Fig.  4.3^  noznallzed  noise 
parameters  S  az]d  IT  are  shown  as  frinctlons  of  normalized  dlstaiuse  from 
the  cathode  with  frequency  as  a  parameter.  Because  of  the  dominating 
role  played  "by  the  convection  current  density  fluctuations,  "both  S  and  n 
seem  to  have  a  similar  spatial  variation  as  As  a  result,  the  sum  S+H 
Is  far  from  Invariant  with  respect  to  distance,  contrary  to  the  one- 
dimensional  linear  multivelocity  electron  heam  noise  theory  vdilch  predicts 
an  Invariance  of  S+n  for  electron  beams  with  a  half 4taxwellian  velocity 
distribution.  It  m^lst  be  pointed  out  that  the  velocity  distribution  In 

a  two-dimensional  diode  is  no  longer  heiLf Maxwellian  when  the  space-charge 
potential,  depression  effects  are  taken  Into  account.  Beyond  the  multi¬ 
velocity  region  little  change  in  II  takes  place  and  the  kinetic  power 
remains  slightly  negative. 

Hie  ratio  of  II  and  S  for  the  three  frequencies  sure  shown  In  Pig. 

4.4  as  a  function  of  normalized  distance.  Unlike  the  results  from  the 
Monte  Ceurlo  calculations  presented  In  the  next  chapter,  the  absolute 
value  of  n/S  remains  less  than  0.2  most  of  the  time. 

4.2.3  Theoretical  Minimum  Noise  Figure  In  &  Two-Dimensional 

Diode.  The  equivalent  minimum  noise  figures  computed  In  terms  of 

the  noise  parameters  S  and  II  are  shown  in  Fig,  4,5  as  a  function  of 
normalized  distance.  little  or  no  reduction  from  thermionic  emission 
shot  noise  is  observed  across  the  two-dimensional  diode.  At  2  Grc 

(a  =  0.25)  the  small  Increane  of  noise  figure  near  the  anode  Is  probably 
due  to  small  errors  accumulated  In  the  process  of  numerical  Integration. 
The  second  derivative  term  of  the  Taylor's  series  employed  In  the  numer¬ 
ical.  step-by-step  forward  integration  process  Is  foimd  to  be  quite 
Influential  beyond  z/z^  «  10.5  for  low  frequencies.  Ihe  minimum  xx>lse 
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flgiire  versus  oormallzed  distance  curves  also  take  on  a  similar  shape  as 
the  convection  current  density  fluctuation  curves. 

In  Fig.  k,6,  the  minimum  noise  figure  Is  presented  as  a  function 
of  frequency  at  three  z-planes,  located  at  z/z^^  =  1.05,  ■  10.5  and 

z/Zjjj  *  30.5,  the  anode  plane.  At  Ihe  anode  of  the  two-dimensional  diode ; 
a  general  deterioration  of  noise  performance  from  I.5  dh  to  2.5  dh  la 
observed . 

h,3  The  Two-Dimensional  Coupling  Mechanism 

After  the  density  calculation  results  are  reviewed,  an  attempt  Is 
made  to  correlate  the  two-dimensional  electron  heaia  noise  transport 
characteristics  with  the  coupling  mechanism.  From  Eq.  A. 9,  the 
z-dlrected  a-c  electric  field  at  (y,z)  Induced  by  a  current  sheet  at 
(y',z)  Is 

00  00 

=  2jfe"  f  f  P'z’  dy'  dz'  ,  (A.9) 

^  "00 

where  p '  =  joj/ z  | . 

In  the  one -dimensional  noise  transport  analysis,  the  coupling 
z-directed  a-c  electric  field  based  on  the  open-circuit  assumption  may 
be  expressed  as  follows; 

00 

E^^(OHE.D)  =  J  F^(z,z)  z  dz  .  (4.1) 

^  —00 

A  comparison  of  the  two  expressions  shows  that  the  two-dljnenslonal 
Induced  a-c  electric  field  depends  not  only  on  the  convection  current 
density  of  the  original  disturbance,  but  It  also  varies  Inversely  as  the 
velocity  of  the  movlxig  electrons.  Fluctuations  eussoclated  with  slow 
electroxis  are  given  more  wel^t  In  the  a-c  coupling,  equlvulent  to  an 


;CAL  MZNIMAl  NOISE  FIGURE  VS 
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effectlve  alteration  of  the  d-c  dis-brlhixtlon  function  In  the  velocity 
domain.  A  net  reduction  of  the  kinetic  power  carried  hy  the  electron 
heam  will  therefore  result,  causing  deterioration  In  the  noise  character¬ 
istics  of  the  electron  heam. 

The  optimum  ratio  of  the  two-dimensional  a-c  electric  field  to 
the  one -dimensional  a-c  electric  field  Is,  for  y  ^  y' , 
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A  conclusion  may  Tae  reached  based  on  Eq.  k.3  that  the  coupling  between 
beam  layers  is  weak.  The  rapid  decay  of  the  a-c  coupling  fields  in  the 
transverse  direction  prevents  significant  interactions  between  beam 
layers.  In  the  density  function  calcxilations  of  noise  trsuisport  in  a 
two-dimensional  diode,  the  lack  of  reduction  In  shot  noise  may  be 
explained  on  the  basis  of  the  coupling  mechanism  employed. 


CHAPTER  V.  RESUITS  AND  INTERPRETATIONS  OF  THE 


MONTE  CARLO  CAICULATIONS 


$.1  Introduction 

In  the  following  sections  the  results  pertaining  to  the  noise 
transport  characteristics  of  a  two-dlmensloneLl  discrete -layer  space-charge- 
limited  diode  obtained  by  means  of  the  Monte  Carlo  technique  are  presented. 
Special  attention  Is  devoted  to  the  Interpretation  and  the  evaluation  of 
the  data.  An  attempt  Is  made  to  correlate  the  noise  transport  etnd 
reduction  phenotnena  observed  with  the  physical  properties  of  the  electron 
beam  under  various  conditions.  Previous  results  from  the  analyses  enqploy- 
ing  one-dimensional  diodes  of  similar  d-o  configurations  are  often  brought 
forth  for  the  purpose  of  comparison  so  that  the  multi -dimensional  space- 
charge  effects  may  be  revealed. 

Before  the  results  are  reviewed.  It  Is  necessary  to  point  out  that 
the  accuracy  of  the  Monte  Ceurlo  method  depends  heavily  <»  the  equivalent 
size  of  the  sample  selected.  In  the  present  calculation,  sample  collection 
is  carried  out  for  a  period  of  two  thousand  time  intervals.  At,  after  a 
quasl-statlonary  condltimi  Is  established  Inside  the  diode.  The  toted 
length  of  the  collection  period  together  with  the  maximum  edlowable  time 
lag  q  would  serve  to  determine  the  quality  of  the  estimates  made.  If 
Gaussian  error  Is  assumed,  for  q  «  12^,  the  freqviency  resolution  Is  4  Oc 
and  the  eqxilvalent  size  of  the  sample  corresponds  to  a  X-square 
distribution  of  32  degrees  of  freedom.  The  confidence  Interval  of  the 
estimated  spectra  Is,  according  to  Fig.  t  32  percent  at  80  percrat 

probability.  In  other  words,  the  estimated  spectra  fcdl  within  ±  3^ 

-99- 


•101 


percent  of  tlie  true  spectra  80  percent  of  the  time  disregarding  errors 
introduced  hy  the  nonstatistical  approxlnatlons . 

From  the  sane  figure  it  is  shown  that  the  confidence  Interval  nay 
be  reduced  to  ±  20  percent  with  8l  degrees  of  freedom  in  the  X-square 
distribution,  corresponding  to  five  thousand  time  Intervals.  The  slow 
improvement  of  the  quality  of  the  estimates  and  the  expense  Involved  in 
the  lengthy  computations  limit  the  accuracy  obtainable  by  this  statistical 
method  of  noise  transport  analysis. 

All  the  noise  parameters  d,  II,  A  and  S  in  the  frequency  domain 
are  presented  in  normalized  lag  window  average  forms  as  described  in  Eqs. 
3.39  through  3.44.  Spectral  estimations  are  made  at  discrete  frequencies 
from  0  Grc  to  40  Gc  at  2  Gc  intervals.  Because  of  the  4  Gc  frequency 
resolution,  no  additional  information  may  be  obtained  if  the  spectral 
densities  are  evaluated  at  frequencies  less  than  2  Gc  apart. 

Following  the  one -dimensional  diode  model  employed  by  Tien  and 
Moshman,  the  temperature  ratio  T^/T  is  taken  to  be  3 *96,  equivalent  to  a 
theoretical  minimum  noise  figure  of  6.96  db  at  the  cathode  plane.  Limita¬ 
tion  in  computing  time  prevents  numerical  evaluation  of  the  noise 
parameters  at  the  cathode.  However,  since  the  random  numbers  employed 
in  the  mathematical  model  are  statistically  checked  beforehand,  there  is 
no  reason  to  doubt  the  behavior  of  the  noise  parameters  at  the  plane  of 
emission.  From  the  initial  conditions  at  the  cathode,  the  fluctuations 
in  kinetic  potentiea  and  current  are  uncorrelated  at  emission  which  leads 


to 


t(n,o)  .  ♦(0,0)  -  s(o,o)  -  1 
Il(Q,0)  .  A(n,o)  -  0  . 


(5.1) 
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Noise  paraioeters  axe  coaputed  for  four  z-planes  located  at 

-  1*5 

Zg/z^  -  3.5 

•  30.5  >  (5.2) 

where  z  Is  the  average  distance  between  the  cathode  and  the  potential 
n 

minima  of  layers.  Nos.  2  and  4. 

z  ■  6  X  10“®  meter  . 

m 

Due  to  the  variations  in  distance  between  the  cathode  and  the  potential 
minima,  spectra  at  the  first  z-plane  *  I.5)  for  different  beam 

layers  should  be  viewed  separately  to  avoid  unnecessary  confusion.  At 
the  first  z-plane,  the  effective  z/z^^^  changes  from  layer  to  layer. 

As  pointed  out  in  Chapter  II,  the  noise  parameters  are  originally 
defined  for  small-signal  one -dimensional  electron  beams  with  no  velocity 
spread.  At  the  first  two  z-planes  where  the  multivelocity  nature  of  the 
electron  beam  prevails,  the  results  presented  are  actually  the  noise 
parameters  evaluated  after  the  electron  beam  is  given  a  sudden  acceleration 
to  reduce  the  velocity  spread  to  the  meEUi  velocity  ratio.  As  was  shown 
previously,  both  the  current  fluctuations  and  the  kinetic  potential  in 
an  electron  beam  are  continuous  over  a  sudden  change  in  d-c  potential. 

Due  to  the  rapid  change  in  potential  gradient  in  the  vicinity  of 
the  cathode,  it  is  necesseury  to  reduce  the  basic  time  Interval  for 
numerical  Integration  steps  to  4  x  10-^®  seconds  for  the  electrons  inside 
the  first  compartnient  of  each  of  the  beam  layers.  For  the  two  thousand 
time  intervals  the  average  ccmvection  current  density  measured  at  the 
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anode  plane  agrees  with  the  conputed  current  density  corresponding  to  the 
average  potential  minima  of  the  layers  to  within  five  percent.  Because  of 
the  half -Maxwellian  velocity  dlstrlhutlcai  of  the  electrons  any  fluctuation 
of  the  potential  minima  about  their  respective  averages  tends  to  Increase 
the  number  of  electrons  crossing  slightly  beyond  the  value  corresponding 
to  the  average  potential  minima.  The  five  percent  discrepancy  Indicated 
is  not  considered  to  be  very  significant  due  to  the  short  time  average 
employed  in  the  estimation  of  the  average  potential  profile. 

3.2  The  Two-Dimensional  Potential.  Distribution 

Both  the  space -charge  distribution  and  the  potential  profile  of 
the  two-dimensional  space-charge-limited  diode  are  shown  compared  with 
their  respective  one -dimensional  predictions  in  Figs.  5*2  and  5.3.  A 
great  deal  of  simllaority  is  found  between  the  curves  for  the  center  layer 
of  the  two-dimensional  model  and  their  one -dimensional  counterpeurts . 
However,  an  appreciable  amount  of  variation  In  both  space-charge  density 
and  d-c  potential  are  observed  in  the  transverse  direction,  which  causes 
a  50  percent  difference  In  current  density  between  the  center  layer  and 
the  edge  layers. 

As  a  direct  result  of  the  two-dimensional  space-charge  effect, 
certain  discrepancies  in  d-c  potential  are  found  between  nei^boring  layers 
at  the  same  distance  from  the  cathode.  In  the  region  beyond  the  potential 
minima,  each  layer  has  its  own  low  velocity  cutoff;  in  other  words, 
additional  velocity  spread  Is  effectively  induced  and  the  nonlinear  region 
of  the  electron  beam  is  greatly  extended.  At  the  four  z-planes  the 
average  d-c  potentials  are  -0.12  volts,  +0.12  volts,  l.T  volts  and  10 
volts,  respectively.  The  existence  of  the  artlflcleUL  multivelocity  region 
is  shown  vital  to  noise  reduction  later  in  this  chapiter. 
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FIG.  5.3 
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5.3  Correlatloo  Pimctlons 

Auto-correlatloa  fvtnctlois  for  both  the  kinetic  potential  and  the 
current  fluctuatl(»is  In  the  time  domain  are  displayed  In  Figs.  3*^  and  5«3* 
Appreciable  deviations  from  the  Ideal  delta  function  are  observed.  A 
certain  amount  of  correlatlcm  has  developed  betveen  fluctuations  sepeurated 
by  short  time  intervals.  Prom  these  two  figures,  negligible  amounts  of 
correlation  between  fluctuations  are  detected  farther  than  30  time  Inter¬ 
vals  apart;  therefore,  the  123  At  maximum  time  lag  assumption  seems  to  be 
Justified. 

1:1  Noise  Parameters 

Noise  peirameters  are  generally  presented  In  two  different  forms. 

They  are  shown  first  as  a  function  of  frequency  with  the  distance  from 
the  cathode  as  a  parameter.  In  addition,  they  are  plotted  against  distance 
for  a  selected  number  of  discrete  frequencies.  Comparisons  between  noise 
parameters  of  individual  beam  leiyers  and  those  of  the  whole  electron  beam 
are  also  made  at  the  anode. 

3.4.1  Self-Power  Spectral  Density  ( SPSP) .  The  normalized  self¬ 
power  spectral  density  9  of  the  kinetic  potential  is  shown  as  a  functioi 
of  frequency  and  then  as  a  function  of  normalized  distance  In  Figs.  3.6 
and  3«7>  respectively.  At  the  first  z-plane,  ♦  is  approximately  unity  as 
expected.  In  agreement  with  the  theory  that  the  gating  action  of  the 
fluctuating  potential  minimum  does  not  alter  the  velocity  distributlcxi 
of  the  electrons.  Formation  of  the  peak  and  the  valley  of  4  Gc  and  8  Oc 
seems  to  have  originated  In  the  low  potential  reglcm  Just  beyond  the 
potential  minima  and  the  causes  of  their  existence  are  not  yet  known. 
Probably  these  points  are  outside  of  the  confidence  Interval  by  chance  sus 
indicated  at  the  beginning  of  this  chapter. 
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The  SPSD's  of  kinetic  potential  are  fovmd  to  be  gronring  funetiooi 
with  respect  to  distwce  and  d-c  potential  which  is  siinilar  to  the  results 
obtained  from  the  one -dimensional  multivelocity  bean  analysis. 

In  a  slmilEur  manner  the  normalized  self- power  spectral  density  T 

of  the  current  fluctuations  is  shown  for  the  whole  beam,  the  edge  layer 

and  the  center  layer  as  a  function  of  frequency  in  Figs.  5*8/  5.9>  5*10 

and  ^.11,  respectively.  Also  Included  are  the  respective  shot  noise 

reduction  factors  calculated  by  Tien  and  Moshman  and  by  Dayem  and 

Lambert.  Since  the  SPSD  of  the  kinetic  potential  is  essentially  unity  in 

the  vicinity  of  the  potential  minima,  the  normalized  SPSD  of  the  current 

fluctuations  taken  at  z/z  >1.5  may  be  regarded  equivalent  to  P^.  Results 

m 

from  the  two-dimensional  diode  are  found  to  be  quite  different  from  those 
obtained  for  either  the  open-circuit  or  the  short-circuit  one-dlnensicual 
diode  except  at  the  extremely  high-frequency  limit  where  there  is  little 
reduction  in  shot  effect.  No  sharp  rise  and  fedl  similar  to  Tien's  dip 
is  observed  in  the  self -power  density  spectrum  of  current  fluctuations. 

At  the  anode  plane  there  is  a  substantial  amount  of  reduction  in  current 
fluctuations  below  the  full  shot  noise  level  up  to  40  Gc  which  is  at  least 
ten  times  the  electron  plasma  frequency  corresponding  to  the  electron 
concentration  at  the  potential  minimum  of  the  center  layer. 

A  comparison  is  made  for  the  SPSD  of  the  current  fluctuations  for 
the  edge  layer,  the  center  layer  and  the  whole  beam  at  the  anode  plane. 

At  low  frequencies,  a  lesser  amount  of  reduction  in  shot  noise  current 
takes  place  at  the  Indlvidueil  layers  compared  to  the  r-f  current  reduc¬ 
tion  in  the  whole  beam.  This  phenomenon  may  be  readily  explained  by 
recalling  the  a-c  coupling  mechanism  described  in  Chapter  III  in  which 
the  a-c  electric  field  is  considered  to  be  a  function  of  the  tota^.  a-e 
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The  nomcLLlzed  self -paver  speetrcd  density  of  the  current  fluctu¬ 
ations  in  individual  bean  layers  and  the  vhole  bean  are  presented  as  a 
function  of  nomalized  distance  in  Pigs.  ^.12  and  5«13>  Very  little 
reduction  in  noise  cvirrent  is  found  beyond  the  potenticd  nininun  at  lov 
frequencies.  On  the  other  hand  the  gating  action  of  the  potential  ninisiun 
does  not  effectively  cut  down  the  high  frequency  component  of  the  noise 
current  which  is  in  turn  rapidly  damped  in  the  low  potential  multivelocity 
region.  Beyond  the  10. 3  z^  plane,  where  the  average  potential,  reaches 
1.7  volts,  V  settles  down  and  becomes  Independent  of  distance  as  predicted 
by  the  single  velocity  theory. 

Shollar  to  the  results  obtained  by  Dayem  and  Lambert,  the  reduction 
in  current  in  the  low  potential  region  is  accompanied  by  a  rise  in  kinetic 
potential,  contrary  to  the  reduction  mechanism  offered  by  the  gating  motion 
of  the  potential  minimum. 

$.4.2  Cross -Power  Spectral  Density  (CPSP).  The  reed  and  the 
imaginary  part  of  the  cross-power  spectral  density  between  the  kinetic 
potential  and  the  current  fluctuations  cure  presented  in  Figs.  3.14  throxxgh 
3.16,  respectively.  According  to  the  postvilated  input  conditions,  both 
n  eind  A  are  equal  to  zero  at  the  cathode  plane.  In  the  retarding- field 
region,  heuxlly  any  correlations  are  developed  between  the  two  types  of 
fluctuaticms  as  evident  by  the  small  values  of  n  and  A  at  the  first 
z-plane.  Negative  values  for  n  are  observed  at  the  low  frequency  end  of 
the  spectrum  In  the  vicinity  of  the  potential  minimum,  which,  according 
to  Appendix  C,  Indicates  that  the  magnitude  of  the  slow  wave  is  leurger 
than  that  of  the  feist  wave.  Some  kind  of  damping  or  gradual  reflection 
must  have  happened  to  the  fast  space -charge  wave  in  the  low  potentied 
region,  leaving  an  excessive  amount  of  slew  wave  with  negative  power. 


DISTASCE  FROM  THE  CAXHCSS 


FIG.  NaRMALI2:&D  NOISE  PABAMBTER  A  VS.  FREQUENCI  (THE  WHOIE  BEAM) 
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It  Is  observed  that  the  real  part  of  H,  the  ooraallzed  eross^power 
spectral  density,  varies  between  0.32  and  0.66  at  the  anode  plane.  At 
low  frequencies  H  Is  asymptotic  to  0.4  which  Is  rather  high  compared  to 
the  results  by  Watkins,  Slegman  and  Hsleh^^  or  Dayem  and  Lambert^^.  The 
rapid  development  of  the  large  amount  of  correlation  In  the  region  beyond 
the  potential  minimum  la  likely  to  be  a  direct  result  of  the  d-c  velocity 
spread  Injected  by  the  two-dimensional  space-charge  effects.  The  CPSD's 
of  the  center  and  edge  layers  at  the  anode  plane  are  found  to  be  slightly 
less  than  the  previously  mentioned  figures,  therefore  less  correlation 
between  the  kinetic  potential  and  the  current  fluctuations  Is  expected. 

3.4.3  Noise  Parameter  S.  (kily  a  small  reduction  In  S  takes  place 
at  the  center  beam  layer  across  the  two-dimensional  diode.  In  Pig.  5 ‘I?; 

S  of  the  center  layer  Is  shown  to  be  relatively  Independent  of  frequency. 
On  the  other  hand,  similarities  are  found  between  the  noise  parameter  S 
versus  frequency  curves  Illustrated  in  Figs.  5 *18  through  5*20  for  the 
edge  layer  and  the  whole  electron  beam  and  the  self-power  spectral  density 
\|r  versus  frequency  curves  in  Pigs.  5*8  and  5-9 •  At  zero  frequency,  S  is 
reduced  to  0.43  at  the  anode  plane  compeured  to  0.23  obtained  by  Dayem  and 
Lambert.  At  very  high  frequencies,  S  approaches  unity  m  expected. 

In  the  small-signal  single-velocity  theory,  both  S  and  11  are 

invariants  under  any  d-c  acceleration.  However,  these  two  noise 

parameters  are  found  to  be  noiconservatlve  even  beyond  the  z  ■  10.3  z 

m 

plane  of  the  two-dlmenslonsd  beam  model.  Ihe  noise  parameters  sure  showti 
In  Pig.  3*21  fts  a  function  of  distance.  In  the  multivelocity  region, 
nearly  Identical  growth  rates  of  S  and  H  sure  observed  at  some  ehoemi 
frequencies  aM  the  cause  Is  tmced  back  to  the  growth  phenosMacn  cX  • 
or  the  Increase  In  aagnltude  of  kinetic  potmtlal  with  dlstanee.  A 
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perelstent  antount  of  addltloial  velocity  spread  Induced  by  the  space-charge 
potentlcd  depression  may  change  the  noise  paraaeters  S  and  n  by  the  sane 
amount  upon  acceleration,  leaving  S-II  practically  Invariant. 

A  similar  display  of  the  noise  parameters  3  and  H  may  be  seen  In 
Figs.  5.22  and  5*23  for  the  center  layer  and  the  edge  layer  of  the  ribbon 
beam.  More  gentle  elopes  for  the  S  and  n  curves  are  observed.  As  It  was 
previously  pointed  out,  the  edge  Icyer  Is  actually  ferther  away  from  Its 
potential  mlnljmim  compared  to  the  center  layer  at  the  same  normalized 
z-plane.  Consequently,  less  variations  In  magnitude  are  found  In  the 
noise  parameter  vers\is  distance  curves  for  the  edge  layers. 

Noise  Parameters  S-n,  S+H  and  n/S.  Prom  the  smeLll-slgaal 
noise  theory  by  Haus,  the  theoretical  minimum  noise  figure  of  a  longi¬ 
tudinal  beam-type  amplifier  is  directly  proportional  to  the  difference  of 
the  two  normalized  noise  parameters  S  and  n 

■W  •  •  (5.5) 

The  same  equation  may  also  be  written  as 


separating  the  bracket  into  two  independent  terms  S  and  1  -  ll/S.  There 
has  always  been  a  controversy  as  to  how  noise  reduction  is  accomplished, 
whether  S  or  the  factor  (l  -  H/s)  is  responsible  for  the  nmconvervative 
noise  transport  phenomena  in  electrcm  beams. 

Curves  for  S-lI,  S+Il  and  II/S  at  the  anode  plane  are  presented  in 
Figs.  5.24  and  5<25,  respectively.  In  a  linearized  one-dimensionsd 
emalysis  both  S-lI,  the  amount  of  negative  energy  carried  by  the  slow 
space-charge  wave,  and  S+H,  the  amount  of  positive  energy  carried  by 
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the  fast  space-charge  wave,  should  stay  ccaistant  throughout.  For  the 
tvo-dlnensioiml  diode  S-Il  is  found  to  be  under  0.4  in  the  frequency  range 
of  Interest  and  up  to  3  db  increase  in  S+n  is  observed  between  10  and 
40  Gc.  Accompanying  the  growth  of  II,  the  Increase  in  fast  wave  amplitude 
is  attributed  to  the  multi-dimensional  space-charge  effects. 

Unprecedentedly  large  valvies  for  the  ratio  Il/S  for  the  two- 
dimensional  model  are  observed,  fluctviating  between  0.6  and  0.9  up  to 
40  Gc.  The  n/S  curves  for  the  individual  layers  are  found  to  be  somewhat 
below  that  of  the  whole  beam,  indicating  that  the  a-c  coupling  effects 
through  the  open-circuit  assumption  is  the  Important  mechanism  in  bringing 
forth  the  correlations  between  the  velocity  fluctuations  and  current  fluc¬ 
tuations. 

The  same  ratio  is  again  shown  as  a  function  of  distance  in  Figs. 
5.26  through  5«28.  For  the  region  beyond  the  third  z-plane,  very  little 
change  in  Il/S  is  observed  for  the  individual  layers.  On  the  other  hand, 
a  substantial  increase  in  Il/S  for  the  whole  beam  tends  to  compensate  for 
the  effects  introduced  by  S  on  the  right-hand  side  of  Eq.  5*4  making  it 
constant  with  distance  at  frequencies  above  4  Gc.  The  decay  in  S(l  -  Il/s) 
for  low  frequencies  will  be  discussed  in  the  foUcwing  section. 

3.3  Theoretical  Minimum  Noise  Figure  in  a  Two-Dimensional  Diode 

Theoretical  minimum  noise  figures  F  computed  in  terms  of  the  noise 
parameters  S  and  II  are  presented  in  Figs.  5.29  through  5«36  for  the  two- 
dimensional  diode  operating  under  conditions  described  in  Chapter  III. 
Approximately  3  db  noise  reduction  is  found  at  the  anode  plane  for  the 
edge  and  the  center  layer  of  the  electron  beam  up  to  20  Gc.  At  low  fre¬ 
quencies,  a  small  rise  in  F  is  observed  beyond  the  multivelocity  region 
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PIG.  5.29  THEORETICAL  MDEEMUM  NOISE  FIGURE  VS.  FREQUEHCY  (CEREER  lAXER  OF  USE  ELECTRON  BEAM). 
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FIG.  ^.31  THEORETICAL  HDIIMUM  NOISE  FIGURE  VS.  FREQUENCY  (THE  WHOIE  BEAM) 
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posslbly  caused  by  the  Increase  in  S  previously  mentioned.  As  for  the 
two-dimensional  beam  as  a  whole,  less  than  3  db  minimum  noise  figure  Is 
predicted  for  frequencies  below  20  Gc  and  It  stays  under  4  db  up  to  to  Gc. 
As  expected,  the  noise  figure  of  the  two-dimensional  diode  tends  to  rise 
with  frequency,  reaching  full  shot  noise  level  at  approximately  70  Qc. 

Results  previously  obtained  for  (xie-dlmenslonal  noise  transport 
analyses  are  shown  ccnapared  with  the  minimum  theoretical  noise  figure  of 
the  two-dimensional  diode.  For  the  two-dimensional  model  both  the  spectral 
estimates  and  their  lag  window  averages  are  presented.  It  Is  obvious  that 
the  point  representing  a  noise  figure  of  less  than  zero  db  Is  outside  of 
the  32  percent  confidence  Interval.  Aside  from  a  few  points,  the  smoothed 
curve  seems  to  be  a  good  representation  of  the  noise  characteristics  In 
the  frequency  domain.  A  theoretical  minlmiam  noise  figure  of  less  than 
2  db  Is  predicted  for  a  small  frequency  range,  and  from  the  two-dimensional 
analysis  the  prediction  appears  to  be  In  quantitative  agreement  with  some 
experimental  results  recently  reported.  It  is  necessary  to  point  out  that 
the  results  for  the  one -dimensional  analyses  by  Vivian^®,  Dayem  and 
Lambert^®  and  Siegman,  Watkins  and  Hsleh^^  are  evaluated  at  a  small 
distance  beyond  the  potential  minimum  where  the  noise  parameters  are 
assumed  to  become  Invariant  with  distance,  while  the  estimates  for  the 
two-dimensional  diode  are  obtained  at  the  anode  plane.  In  the  multi¬ 
velocity  region  created  by  the  space-charge  effects,  substantial  noise 
reduction  is  observed  for  the  two-dimensional  model  at  low  frequencies. 

It  is  clear  from  the  foregoing  conparlson  that  the  multi -dimen¬ 
sional  effects  may  be  substantlcd.  and  should  not  be  disregarded  In  noise 
transport  analyses  and  the  minimum  noise  figure  of  longitudinal  beam-type 


devices  nay  be  greatly  reduced  when  an  additional  dlnenslon  la  taken  Into 
account . 

^.6  Space-Charge  Effects  at  the  Anode 

In  the  two-dlmensl(»ial  theoretical  model  employed  In  the  Monte 
Carlo  calculation  of  noise  transport  analysis,  the  space-charge  effect 
Is  approaching  zero  at  the  anode.  As  a  result,  the  slopes  of  the  noise 
parfjneters  S  and  H  are  expected  to  level  off  at  the  anode  plane  where 
the  single-velocity  theory  begins  to  take  over.  Unfortunately,  this 
phenomenon  is  not  confirmed  because  of  the  lack  of  information  concerning 
the  transport  characteristics  of  the  elef'tron  beam  between  the  anode  and 
the  third  z-plane.  In  Fig.  5.25,  no  attempt  Is  made  to  fit  the  S  and  IT 
curves  with  a  prescribed  zero  gradient  at  the  anode  plane. 

5.7  Errors  in  Electric  Field  Evaluation 

The  nature  of  errors  introduced  In  the  d-c  electric  field  evalua¬ 
tion  is  actually  two-fold.  First  of  all,  the  method  described  in  Appendix 
B  involves  certain  averaging  processes  when  the  compartments  are  Installed 
Second,  the  replacement  of  the  infinite  series  with  an  Integral  brings 
in  a  further  approximation  in  the  electric  field  computation.  However, 
sampling  of  the  crossing  velocity  at  the  z-planes  shows  good  agreement 
between  the  d-c  potential  and  the  velocity  of  the  slowest  electrons.  This 
together  with  the  average  current  and  potential  minimum  characteristics, 
leads  to  the  belief  that  the  method  for  the  d-c  field  computation  is 
satisfactory. 


CHAPTER  VI .  SXMMARY  AND  CONCLUSIONS 


6.1  Svamnary  and  Conclusions 

Obe  noise  transport  phenonena  In  a  hollow-team  thermionic  diode 
operated  under  space -charge -limited  conditions  have  teen  Investigated  hy 
means  of  two  different  nunerlcal  methods .  A  two-dimensional  mathematical 
discrete  team  model  weis  constructed  to  simulate  the  actual  diode  with  a 
large  team  radius  to  team  thickness  ratio. 

A  small-signal  linearized  Boltzmann  equation  expressed  In  terns 
of  the  two-dimensional  density  function  and  the  electric  lields  has  teen 
derived  for  an  electron  team  with  an  infinite  axial  magnetic  field. 
Because  of  the  discontinuity  of  the  d-c  density  function  in  the  velocity 
domain,  Slegman's  slngulari'fy  extraction  method  weie  employed  to  take 
away  the  singularity  factor  from  the  transport  equation;  details  of  the 
extraction  techniques  are  presented.  To  facilitate  a  solution  for  the 
transport  equation  hy  digits^,  computer  techniques,  velocity  classes  are 
Introduced  as  an  approximation  to  the  distribution  function.  Singular¬ 
ities  brought  forth  by  the  velocity  class  approximations  In  the  boundary 
conditions  were  Investigated  and  their  i^yslcal  meanings  are  discussed  In 
detail.  Linear  transform  techzilques  were  employed  to  solve  the  set  of 
small-signal  multi -coupled  linear  first-order  differential  equations 
governing  the  behavior  of  the  density  functions  In  phase  space.  Haus's 
small-signal  noise  theory  was  employed  to  evaluate  the  noise  transport 
chareuiterlstlcs  of  the  two-dlmenslonsd  electron  beam. 

The  two-dimensional  team  noise  transport  phencmena  was  clLso 
treated  by  means  of  the  Monte  Carlo  techx>ique .  Contrary  to  the  densily 
function  calculations,  the  stochastic  emission  process  at  the  thermi«aic 
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cathode  was  directly  slsmlated  and  the  electrcos  onltted  were  traced  in 
phase  space  by  nunericail  integration  techniques  until  they  left  the 
diode.  Random  number  transform  techniques  were  presented  and  the  merits 
of  the  Monte  Carlo  method  of  noise  transport  analysis  is  discussed.  In 
Appendix  B,  the  method  employed  to  eval\iate  the  d-c  and  a-c  electric 
fields  is  described  and  the  implication  of  the  open-circuit  assunq)tlon  is 
elaborated.  Following  the  thin  beam  asswption,  the  axial  a-c  electric 
field  is  talcen  to  be  uniform  in  the  transverse  direction  eu:ross  the 
electron  beam,  equivalent  to  a  one -dimensional  assumption  for  the  a-c 
electric  field.  As  a  result,  only  the  d-c  multi -dimensional  effects  are 
Included  in  the  Monte  Carlo  calctilations .  A  brief  description  of  the 
generalized  harmonic  analysis  and  an  Introduction  to  sampling  techniques 
are  presented.  Haus's  noise  parameters  are  estimated  statistically  by 
means  of  sampling  techniques. 

Results  frcmi  the  density  function  ceJ.culations  carried  out  at  three 
discrete  frequencies  are  presented.  Little  or  no  reduction  from  the 
themdoiic  field  emission  shot  noise  is  observed.  There  is  no  evidence  of 
any  m\iltl -dimensional  effects  postulated  and  even  the  well  known  noise 
reduction  mechanism  due  to  the  gating  action  of  the  potential  mlnlmvan 
appears  to  be  missing.  A  close  examination  of  the  two-dimensional 
coupling  mechanism  revesds  that  the  induced  a-c  electric  fields  in  the 
two-dimensional  analysis  decay  too  rapidly  in  the  transverse  direction  to 
bring  forth  significant  multi -dimensional  effects .  As  a  coxisequence  of  the 
velocity  dependent  characteristics  of  the  coupliiig  field,  heavy  oiqihasls 
is  placed  on  the  fluctuations  associated  with  the  slow  electrons,  causing 
a  net  reduction  of  the  kinetic  power  carried  by  the  electron  bean  and 
some  deterlation  in  noise  performance. 
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Frcm  the  results  of  the  Monte  Carlo  calculations,  substantial 
noise  reduction  at  a  vide  band  of  frequencies  up  to  20  Gc  is  observed. 

A  minimum  noise  figure  of  less  than  2  db  is  predicted  for  a  annll 
frequency  range  and  the  prediction  appesirs  to  be  in  quantitative  agree¬ 
ment  with  some  experimental  results  recently  reported.  In  addition  to 
the  noise  current  smoothing  mechanism  euiileved  by  the  fluctuation  of  the 
potential  minlmvim,  additional  reduction  of  noise  at  frequencies  several 
times  the  electron  plasma  oscillation  frequency  at  the  potential,  mini¬ 
mum  is  observed  in  the  extended  multivelocity  region  created  by  the 
multi -dimensional  space-charge  potential  depression  effects.  Unusually 
large  values  of  the  ratio  n/S  for  the  two-dimensional  beam  model  are 
realized  in  the  nonlinear  multivelocity  region  beyond  the  potential, 
minlmm,  resulting  in  a  net  reduction  of  the  minimum  achievable  noise 
figure.  Neither  S-II  nor  S+II  is  found  to  be  invariant  with  respect  to 
distance  from  the  cathode  as  predicted  by  the  one -dimensional  linearized 
single -velocity  and  multlvelocii^  noise  theories.  IQie  importance  of 
the  multi -dimensional  d-c  space-charge  effects  are  clearly  brou^t  out 
in  the  Monte  Carlo  analysis. 

It  is  evident  from  the  discrepeuicy  between  the  results  from  the 
two  numerical  methods  of  noise  transport  anelysis  in  two-dimensional 
multivelocity  electron  beams  that  the  coupling  mechanism  plays  em 
influential  role  in  the  transport  characteristics .  Fvirther  improvement 
of  the  electron  baam  noise  transport  theory  would  depend  on  a  more 
thorough  understanding  of  the  actual  coupling  mechanism  in  a  multi¬ 
dimensional  realistic  beam  model. 
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6.2  Suggestions  for  Fvirther  Study 

It  vould  "bs  Interesting  to  replace  the  tvo -dimensional  a-c 
coupling  tern  In  -the  density  function  calculations  by  Its  one -dimensional 
counterpart  and  see  If  the  results  thus  obtalxied  vould  be  similar  to  the 
results  from  the  Monte  Carlo  calculations.  Bie  validity  of  the  linear¬ 
ized  analysis  employed  In  the  highly  nonlineeLT  region  of  the  electron 
beam  will  be  revealed. 

Coupling  mechanisms  in  multi -dimensional  electron  beans  should 
be  further  Investigated.  It  might  even  be  feasible  to  construct  experi¬ 
mental  measurement  to  help  determine  the  coupling  a-c  electric  fields 
induced  by  the  fluctuation  of  convection  current  densities . 


APPENDIX  A.  A-C  AND  D-C  ELECTRIC  FIELDS  IN  A  TWO-DIMENSIONAL  DIODE 


AND  SOME  DETAIIfi  IN  THE  DEHSITI  PUNCTICaJ  FORMULATION 

A. I  A-c  Electric  Fields 

In  a  two-dimensional  diode,  the  a-c  electric  fields  Induced  by 
the  moving  space  charge  may  he  evaluated  using  Maxwell's  equatlcms.  At 
high  frequencies,  because  the  electromagnetic  field  strength  drops  off 
rapidly  with  distance  from  the  origin  of  the  disturbance,  effects  due  to 
the  Image  currents  are  neglected  In  the  following  calculations.  Cc«i- 
slder  l^(y',z)  In  Fig.  A.l  to  be  the  a-c  current  density  of  a  z-dlrected 
sheet  beam  Infinite  In  the  x-dlrectlon. 


00 

Ij^(ySz)  =  -e  y  4  F^{y',z,z)dz  .  (A.l) 

-00 

If  the  convection  current  l^(y',z)  Is  assumed  to  have  an  e”*^^  ^  varia¬ 
tion,  according  to  the  Laplace  law,  the  temgentlal  magnetic  fields 
H^i(y,z)  and  may  be  expressed  as 


H  (y,z)  = 

XI 


Ae' 


P'(y-y') 


for 


y  <  y' 


(A.2) 


H  (y,z)  =  Be"^  ^  for  y  >  y'  ,  (A. 3) 

X2 


wti^ere  b  |<n/u'|  and  3*  /  0.  Continuity  of  magnetic  fields  tangential 
to) a  sheet  beam  of  Inflnltesinua  thickness  requires  that 


i^(ySz)dy'  -  * 
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(A.4) 


FIG.  A.l  CAICUIATIOR  OF  EIKTRIC  FIEUDS  ASSOCEIfG  A 


SHEET  OF  CURHEIfT  AT  y  ■  y' 
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Since  the  z-directed  electric  fields  at  y'  associated  with  either  or 
H  must  Ise  identical. 

X2  * 

A  -  -  B  (A.5) 

8U3d 

ii(y'>z)  dy*  »  2A  =  -  2B  .  (A.6) 

From  Maxwell's  equations 


’  ■■57  ' 

the  z-directed  electric  field  at  (y,z)  induced  "by  a  current  sheet  at 
(y',z)  is  therefore. 


li(y’»2)  -S'lv.v'l 

E  (y,z)  ,  =  -  -rk - e  ^  ^  '  dy' 

iz'*'^  'y'  J2a3€  ^ 


(A.8) 


and  the  total  electric  field  at  (y,z)  is 


“23^;  /  /  P'z’Fi(y',z,z’)  dy'  dz* 

-CO  —00 

This  analysis  assxnnes  P'  =  constant,  i.e.,  no  acceleration. 


(A.9) 


A. II  Potential  Dlstrlhutlon  of  &  Two-Dimensional  Space  -Charge  -Limited 
Diode 

Techniques  similar  to  Langmuir's  employed  to  reduce  Poisson's 
equation  for  a  space -charge -limited  two-dlmensionsd  diode  to  dimension¬ 
less  form.  Unfortunately,  no  general  solution  can  be  obtained  corres¬ 
ponding  to  the  one  obtained  in  the  one -dimensional  case.  Some  arbltrazy 
boundary  conditions  have  to  be  applied  to  the  tcq?  and  bottom  layers  of 
the  electron  beam  to  facilitate  obtaining  a  special  solution. 

For  any  sheet  beam  located  at  y  with  current  flowing  in  the 
z-dlrection  if  the  region  between  the  cathode  and  the  potential  minima  is 


called  Region  I  and  the  region  heyond  is  Region  II,  the  space-charge 
densities  In  respective  regions  are 


Pi  = 


o 

m 


V  F^(v  )dv 
0  0^  o'  0 


oo 

/ 


V  F  (v  )dv  i 
o  0 '  ^  ' 


where  p 
v^(y,z) 


v  F  (v  )dv 
0  o'  o  o 


€C 

I 

V' 

0 


V  F  (v  )dv 

r\  rs^  r\' 


0  o'  0 


(A.IO) 


=  p(y>z)  space -charge  density  at  (y,z), 

=  V®  +  2hV(y,z)j  V  Is  the  velocity  of  electrons  at  (y,z) 
with  v^  as  Initial  velocity, 
a  Initial  velocity  of  electrons, 

=  *7 -2hV' (y),  initial  velocity  of  the  slowest  electrons  in 

the  layer  located  at  y  >dilch  possesses  harely  enough 

energy  to  reach  the  potential  minimiaa, 

=  number  of  electrons  emitted  with  their  Initial  velocity 

between  v  and  v  +dv  per  unit  «u:ea  per  second. 

0  0  0 


The  initial  velocities  of  the  electrons  emitted  perpendicular  to 
a  thermionic  cathode  surface  follow  a  Raylei^  distribution  according  to 
assumption  7  in  Chapter  II. 


F  (v  )dv 
o'  O  0 


-ov! 


-  2a  I-  V  e 
den.  o 


dv_ 


(A.ll) 


^den  emission  current  denslly  at  the  cathode.  For  an  Itmersed 

flow  diode,  mixing  of  electrons  originating  from  two.  different  current 
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sheets  Is  prohibited.  It,  In  addition,  the  effect  of  colUalona  between 
electrons  Is  neglected,  Eqs.  A. 10  and  A. 11  may  be  combined  to  fozm 


/ai(V-V')a 


II 


or 


■geov]  r 

.  “  J  1 

J  ^  J 

r\  r\ 

-v2 


dvl 

“J 


(A.12) 


Poisson's  equation  may  then  be  expressed  in  dimensionless  variables. 


1 

2 


1  ±  erf 


f 


(A.15) 


where 


A.  Ill 


2eaV* 
m  * 

the  "±"  signs  are  for  Regions  I  and  II  respectively. 

A-c  Transport  Equation  In  DlnenslonlesB  Parametric  Foim 
After  the  Introduction  of  the  excess  energy  parameter  W,  Eq. 


2.6 


is  reduced  by  one  term 
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±  -jw  +  £b[V(y,z)  -  V'(y)]|  *1^ 


oz 


SF,(y>z»w)  ^,(yyz»w)>,  f 

- ^  "  2hE^(y,z)  ^  2bEj_(y,z)  |W  +  2h[v(y,z) 


-  V(y)3} 


1/2  aF^(y,z,w) 

iw 


=  0  ,  (A.1^^) 


where  W 
But 


=  £2  _  2h[V(y,z)  -  V'(y)]  . 


^  _  .Oh 

dz  dz 


=  2hE^(y,z)  . 


The  transport  equation  Is  therefore. 


(A.  15) 


JcuFj^(y>z,w)  ±  -^w  +  ai[v(y,z)  -  v'(y)3j- 


1/2  ^^iy>z,v) 


dz 


r  'I  ^Jy>z,v) 

2hE^(y,z)  jw  +  2h[V(y,z)  -  V'(y)]J.  — S-gjj -  =  0  .  (A.l6) 


This  equation  combines  with  Eq.  2.7  to  form  an  a-c  treuisport  equation 
expressed  In  terms  of  the  density  functions  alone. 
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SF  (y,z,W) 

§r— 


-  Ja]g^(y,«»W) 
s/w  +  2h[V(y,z)  -  V'(y)] 


±  r-=^ 


w  w 

// 


P'F^(ySz,W')e 


-P*|y-y'l 


dy* 


(M 


(A.17) 


Each  tern  In  Eq.  A.17  may  be  put  Into  dlmenBlonless  font,  yielding, 


^i(y>z,W) 

§z 


■""ai 


aPj(c*i,w) 

z 


(A.18) 


J<i^’j^(y,z,w) 

>/w  +  2h[V(y,z)  -  V'(y)J 


g  •  a  -^P^(C,6,w) 


w  C>6) 

o'*’  a 


Iq  J8P^(C>6>w) 

®  ■^w+n(5»l) 


(A. 19) 


and 


G 
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he 

2Ja)€p 


CO  60 


dy»  dW 


Pi(CS6»w)  exp  • 


“g|K' 


N/w(c',t)+(£',e) 


Vw(S',0  +  T|(S',6) 


3Po(C»t»w) 

dw 


(A.20) 


The  transport  equation  In  dimensionless  form  is  therefore. 


^5  '/w  +q(5,6) 


/  /  “»[ 
-60  -60 


-...-gu-s'i  1 

'/w(CM)+n(CM  J 


P^(SM,w) 

'/w(5'>i)+n(CS0 


5w 


I  .  (A.21) 
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A.IV  Singularities  In  the  A-c  Density  PunctlonB 

In  a  diode  In  vhlch  electrcDS  hava  a  half -Maxwellian  velocl'^ 
dlstrihution,  the  d-c  component  of  the  density  function  may  he  expressed 
as 

2aljy)  fr  /  “OW4.  -QW* 

F^(y,z,W)  - - ^ j  S(c«^  +  Tic)  -  S(W_^)J  e  +  S(W^)  e 

■  -ciV  ■] 

e  ■  I  ,  (A.22) 

>diere  ^  W  for  the  forward  going  electrons,  z  >  0, 

W  ^  W  for  the  hackward  going  electrons,  z  <  0, 

T]^  =  the  dimensionless  potential  parameter  at  the  cathode, 

S(oiW^+Ti^)  =  a  unit  step  function,  and 

S(oW^+Tlc)  =  1  for 

0  for  aW.  <  n  . 

The  "±"  superscripts  of  the  square  brackets  correspond  to  the  forward 
going  and  the  "backward  going  electrons  respectively,  and  these  two 
"brackets  vanish  in  the  region  "beyond  the  potential  mlnlmvBi.  Due  to  the 
existence  of  a  step  function,  or  step  functions  In  the  d-c  density  func¬ 
tions,  one  or  more  singularities  are  expected  in  any  f\mction  Involving 
the  derivative  of  F^(y,z,W)  depending  on  the  discreteness  of  the  velocity 
spectrum  employed.  If  the  d-c  density  function  is  reduced  to  dimension¬ 
less  forms, 

=  -  2  |s(w_^)e  =  S(w_^)  e 

r  i"  1  Ms) 

+  |^S(w_  +  ti^)  -  S(w_)J  e  "  1“^-”  t  (A.23) 
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Ihe  derivative  of  P^CC^i^v)  vith  respect  to  may  "be  obtained  in  Dirac 
delta  fuziction  form. 


5w~ 


21.(0  fr  1 

- - - - 1  8(w^)  -  S(w^)  e 

o  ^  L  -J 


6(\  +  n^j)  -  B(w^)  -  S(w^  +  Ti^)  +  S(w^) 


-w. 


8(w 


\)  -  sC'f.)  -  S(w_  + 


\)  +  S(wJ 


-w 


}- 


(A.2U) 


vhere  8(w)  is  a  Dirac  delta  function. 


A.V  Kinetic  Potential 

If  the  definitions  of  the  a-c  and  d-c  density  functions  8a«  sub¬ 
stituted  into  the  expression  of  Chu's  kinetic  potential  the  invariant 
characteristic  over  a  d-c  velocity  Jisnp  gap  is  found  missing  once  velocity 
spread  is  considered.  However,  a  generalized  form  of  this  equivalent 
voltage  fluctuation  may  be  derived  if  the  problem  is  approached  from  the 
energy  standpoint. 

In  a  two-dimensional  miiltlvelocl'ty^  electron  beam,  the  kinetic 
potential  (or  equivalent  a-c  voltage)  of  a  sheet  beam  at  y  may  be 
defined  as  the  a-c  energy  expressed  in  electron  volts. 


v^(y,z) 


00 

o'"' '  0 


+ 


-I6l- 


where  I^(y)  “  -  |  Fo(y>*»W)<iW 

i(y)  *=  igCy)  +  i^Cy)  » 


■  I  J  “  I  /  ^i(y»*»w)dw  , 


=  current  density  of  the  beam  layer  at  y. 

The  subscripts  "o”  and  "l"  are  for  d-c  and  a-c  quantities 


respectively. 


Vi(y,0  S  ^  .-_3^^(y)  J  tw  +  2tiV(y,z)l  P^(y,z,w)dW 
°  ^  0 

00 

-  Ii(y)  J  [W  +  2nV(y,z)]F^(y,z,W)dwJ  (A.26) 


and 


Vj^(y,z) 


4hl  (y)  _ 
o  0 


CO 

[  /»->, 


z,W)dW 


I  J  F^(y,i,W)dW  J  Wjiy,z,\l)m 


I<,(y) 


] 


4hl 


J  (W^g)  F^(y,z,W)dW  ,  (A. 27) 


where  W 


w 

y'  •  WF^(y,ZFW)dw 
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A. VI  Evaluation  of  the  Weighting  Factor 

Evaluation  of  the  weighting  factor  R( ^,w)  for  the  honogeneoua 
peurt  of  the  a-c  density  functions  for  an  electron  beam  with  a  continuous 
velocity  spectrum  may  be  carried  out  In  the  following  fashion.  If  the 
weighting  factor  Is  made  up  of  two  factors  A(0  and  e”'^,  each  a 

function  of  a  single  variable >  by  means  of  the  shot  noise  formula  and 
Eq.  2.18  for  unit  a-c  current  Input  per  unit  bandwidth,  the  following 
relationship  may  be  derived. 

2ell^(0lAC  =  <if(C,0)> 


or 


R(?,w) 


2etl^(C)l 

~a1  ® 


(A.28) 


(A. 29) 


The  value  of  the  weighting  factor  R(^,w)  thus  obtained  also  satisfies 
the  boundeoy  conditions  of  the  kinetic  potential  at  the  cathode.  From 
Eq.  2.2k 


<  V2(C,0)  > 


(w2-2w+1)  ^  R{  C,w)dw 
■*■0 


( 


•) 

eyt)  J 


00 


J 


(w*-2w+l) 


ge|lp(S)| 


e  dw 


2e|lJC)ld; 


(A.30) 
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aad 

<  V^(C,0)  1*(S,0)  >  -  0  .  (A.31) 

Similarly,  the  kinetic  power  theorem  for  the  electron  beam  as  a  i^le  may 
also  be  verified. 


<  V(0)  1*(0)  >  -  0  .  (A.32) 

A .VII  Points  of  Reflections 

In  the  vicinity  of  the  reflection  points  where  w+t)  approaches 
zero  as  a  limit,  the  first  tern  on  the  right-hand  side  of  Eq.  2.13  in 
Chapter  II  tends  to  become  infinity.  It  is  necessary  to  show  that  the 
incremental  change  of  Pj(5»l»w)  stays  in  bound  at  the  reflection  points 
and  the  transport  equation  by  itself  is  valid  anywhere  in  the  diode. 

If  digital  techniques  are  not  employed  in  the  numerical  calculations,  the 
assumption  of  w-tq  =  0.002  is  absolutely  lumecessary. 

On  the  right-hand  side  of  Eq.  2.13;  all  terms  except  the  first  one 
drop  out  in  the  neighborhood  of  a  reflection  point  for  a  velocity  class  w^ 
leaving 


or 


^  -Jgd6 


(A.55) 


(A.ik) 


where 


11m  +  ti(5,0  -*0. 

5  -"6, 
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As  an  exaople,  the  potential  distribution  cvurve  for  a  one -dimensional 
beam  is  approximately 


»l(6)  =  J  +  0.045  I*  +  0.0045  I®  . 


(A.55) 


Equations  A. 54  and  A. 35  combine  to  give 


dp^(5,|,w^)  ^ 

+  '■■■■■'  — — — 1^  SB  I  Ml  11  . . I  ■  II  I  I  n 


p,(5>s>w.)  rr7T~7.  TTTTs  ' 

^  i  (iMp  +  0.045{|*.|J) 


(a.36) 


■j 


Pi, 


I 


"ML 


2  i  Sp 


(A.37) 


11^  ±  /n  ~ 

I  -»6j.  Pi,r 


=  -  J26  in 


{ 


-  it 


^r 


(A.38) 


j^?e_  •'*  [’  (nT^  *  t )] 


or 
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liffl 


®i,r+i 


(A.ltO) 


A6  -»0 

that  Is,  the  increneatal  change  of  Is  zero  vhen  t)+v^  -»0  for  the 


velocity  class. 


APFENDK  B.  POTENTIAL  DISTRIBUTICai  IN  A  TtfO-DlMEaiSIOKAL 


DIODE  AND  GOSNERAL  HARHCRIC  ANALISIS 


B.I  Pie  D-c  Electric  Field  and  Potential  Distribution 

Within  a  tvo -dimensional  discrete-beam  diode  like  the  model 
en^loyed  In  the  density  function  method  of  noise  transport  Investiga¬ 
tion,  It  Is  possible  to  solve  Poisson's  equation  for  the  Instantaneous 
space -charge  potential  distribution  If  the  solutions  of  a  set  of  Laplace 
equations  In  Dlrlchlet  form  corresponding  to  the  speu^e-charge  distribu¬ 
tion  In  the  diode  are  superposed  on  each  other.  Unlike  the  point  charge 
assumption  In  a  real  three-dimensional  model,  electron  are  taken  to  be 
infinitely  long  lailform  line  cheurges  extended  In  the  x-dlrectlon.  For  a 
line  charge  located  at  Poisson's  equation  may  be  written  as 

**'  **  p 

~  f  f  BCy-y^)  ^  *  (b*i) 

0  "«o  ^ 


where  V(y,z) 
a 


&(z-z^) 


BCy-y^) 


f  6(y-yjj)dy 


the  d-c  potentlsLL  at  (y,z), 

the  distance  between  the  cathode  and  the  anode, 
the  charge  per  unit  length  in  the  x-directlon. 


00 

at 

z 

= 

^0' 

0 

for 

z 

00 

at 

y 

= 

0 

for 

y 

^0' 

1, 


00 

/ 


6(z-z^)dz 


1, 
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vlth  the  following  boundary  coadltlona 

V(0,y)  -  0  , 

V(a,y)  =  0  , 

V(z,±«.)  =  0  , 


3v(z,y^) 


Sv(z,y^) 


y<y. 


y>y. 


(B.aa) 

(B.2b) 

(B.2c) 

(B.2d) 


If  the  Laplace  equation  Is  to  be  solved  by  the  separation  of  variables 
method,  a  general  solution  of  the  potential  function  V(z,y)  Is  found 
such  that 

y/  njt(y-y  )  \ 

sin  exp  ^  ,  (B.J) 

n 

where  the  sign  is  for  the  region  y  >  y^, 

"+"  sign  Is  for  the  region  y  <  y^,  from  symmetry. 

The  Fourier  coefficients  may  be  evaluated  by  means  of  the  boundary 
condition  of  Eq,  B.ad. 


2A^  y  ^  Sln2  SS  dz  =  ^  J  6(z-z^)  sin  Sp  dz  (B.k) 


or 


p,  ,  nnz 
A  =  ^  .  A.  Bln  — a 

n  €  nn  a 
o 


(B.5^ 


and  therefore 
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Unfortunately,  the  Fourier  series  representation  of  potential  dlsti'lbutlm 
function  V  thus  obtained  Is  not  differentiable  at  some  boundaries  irtiere 
texm-by-v^xn  differentiation  of  V(y,z)  results  In  a  divergent  series.  The 
electric  field  Is  not  uniquely  detemlned  at  z  >  0  and  z  »  a. 

In  order  to  arrive  at  a  unique  solution  at  all  boundaries,  the 
image  method  Is  employed  to  solve  for  the  z -directed  d-c  electric  fields 
Instead.  A  line  charge  is  placed  at  betveen  two  infinitely 

extended  conducting  plates  idilch  are  perpendicular  to  the  z-axls  and 
sepaieited  by  a  distance  a.  If  both  electrodes  are  kept  at  zero  potential, 
the  potential  relative  to  either  plate  at  a  point  (y,z)  Inside  the  diode 
due  to  the  Influence  of  the  line  charge  may  be  computed  by  linear  addition 
of  the  space -charge  effects  firom  the  line  charge  Itself  and  Its  Infinite 
number  of  mirror  images  shown  in  Fig.  B.l. 


w 

V(y,z)  =  y  in>/(z-2na-z^)2  +  (y-y^)* 


n=-o» 


y  pJre”  )s  +  (y-y  )2 


n=-  «• 


-  y 


pje  /  {z-2a&+z^)  +  (y-yp) 

^  J  (z-2na-z  )2  +  (y-yj® 


-  (z+z^)=  +  (y-y^j)®  1 

.  (z-z^)®  +  (y-y^)®  J  * 


(z+z^-ana)®  +  (y-y^)®-i 
(z-z^-2na)®  +  (y-y^)®  J 


(z+z^-2na)2  +  (y-y^)*  t 

(z-z^-2na)®  +  (y-y^)®  J  ' 
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HFIHITBLr  EXSBUSD  FABALIXL  lUCTRQISS. 
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Because  (y-y^)®  «  (2na)®  for  n  0 


V(y,z) 


r  r  +  (y-yp)^ 


+  2 


y 


1  ^  r  *+*«  1 

]-  y 


Us  -00 


-z+2  -aia  -]  >1 


(B.8) 


The  series  on  the  ri£^t-hand  side  are  uniformly  convergent  as  shown 
In  Section  B.I.l.  As  an  approximation  the  summations  may  be  replaced 
by  Integrals  In  which  n  Is  regarded  as  a  continuous  Integration  variable. 
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But, 


n®T  ^ 

nd.  n>i 


J  ^  ^  [2na-(2+z^)]  lnt2na-(*4«^)] 

n>i 


n=^r 

-  [2iia-(z+z^)]  j- -  i  [2na-(z-z^)]  ln[aaa-(z-z^) ]  -  [2tia-(z-z^) ]j- 
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-  [2Ta-(z-z^)]  in  [2Ta-(z-Zp)]  +  [2ra-(z-z^)] 


+  [2a-(z-zJ]  in[2a-(z-z^)]  -  [2a-(z-z^) ]j- 


=  7  T[21'a-(z+z  )  ]  in[2Ta-(z+z  )]  -  [2ra-(z-z  )]  int2Ta-(z-z^)  ] 

S  I  o  o  o  o 

-  [2a-(z+z^)]  in[2a-(z+Zp)]  +  [2a-(z-z^)]  ln[2a-(z-z^) .  (B.IO) 


Similarly, 
n«?r 


r  r  2na+(z+z  )  -i  r 

nci 


-  [2Ta+<z-z  )]  ln[2Ta+(a-*  )]  -  [2a4(z+z  )]  In[2a4<z+zj] 


+  [2a+(z-Zjj)3  in[2a+(z-z^)]j.  .  (B.ll) 
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The  last  tvo  terns  combine  to  give 


11m 

T-*<» 


n»T 

{s  /■■ 

n=i 


2aa-(z+z^) 

2na-(z-z  ) 
o 


d(2na) 


2na+z+z^  n 
2na+z  -z^  J 


Hiw 

T-»«» 


+  (z+z^)  in 


+  (z-z^)  in 


[ 

[ 


i  {at.  in  [ 


2Ta+(z+z^)  1 
2Ta-(z+z  )  J 


2Ta-{z“Z^) 
2Ta+(z 


:v  1 
-*-) . 


■I-  2a  in 


=  2  in 


(2a)®-(z+z^)^  J  * 


[2Ta-(z+z^)][2Ta+(z+Zp)l  -i 

[2ra-(z-ZQ)]t2r*4<z-Zo)  ]  J 


+ 


2a-(z+z^) 

2a4{z+z^) 


+  (z-z^)  in 


2b+(z-z^) 

2a“(z-z  ) 
o 


[ 2a-( z -z^ ) ] [ 2a+( z -z^ ) ] 
[  2a-(  z +2^ )  ]  [  2a+(  z +z^)  ] 


[  2a+( z -z^ ) ] [ 2a- ( z+z^ ) ] 
[Sa-Cz-Zg) ] [2a+(z+z^) ] 


] 

] 

]} 

■] 


Zq  r  [2a-(z+z^)][2a-(z-ZQ)]  - 
a  L  [2a+(z+z^)][2a+(z-Zjj)]  . 


The  potentled  at  (y,z)  Is  therefore 


V(y,x) 


(z4-a^)«  •«•  (y-yp)^ 

(*-*^)*  +  (y-y^)® 


+ 


2 

a 


in 


p  (2a)«  -  (i-*q)® 

*  “  L  (»)•  -  (iwo)* 

[2a+(*-*Q)][2a-(s+z^)]  n 

[2a-(z-*o)]l2a+(x+i^)  ]  J 


*o  .  r  t2a-(x+x^)l  [2a-(x-Xp)]  -| 

'‘"a”  L  [2a+(x+z^)]  [2a+(*“*o)l  J  J 


(B.15) 


For  x  +  z  <  0.2b 
o~ 


(2a)2 


in 


(x-xJ2 


(2a)®  -  (x+z^,)' 


c-^y-c-^j  ■ 


hzz  xz 

’  7f  ’ 


i-ln 


[ 


[2a+(z-z^)  ]  [2a-(x-w^)  ] 
t2a-(z-z^)][2a+(x-»*^)] 


] 


2.  in 


r  2x  T  "fiZZ^ 

-in  1  -  •»  “T“  >  (B,15) 

ft  ft  • 
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4  in 


[  2a-(  z  +  z^)]  [  2a-(  z  -z^)  1 
[  2a+(  2 + z-  )]  [2a+(  z  -z  - )  ] 


z  r 


z+z 
_ 0 

2  r 

=  In  • 

L  - 

a 

a 

*  *n 

a  L 

'  ‘J 

2zz 

_ g 

a2 


!I!he  last  three  terms  add  up  to  "be  -(jaz^/a®) . 


(B.l6) 


P.  r  r  (z+z  )2  +  (y-y  )®  -i  5zz  ■> 


and 


1:^1  <  0-12  . 

Since  a  majority  of  the  electrons  inside  the  diode  at  any  moment  are 
within  0.2a  frcmi  the  cathode,  the  simplified  equation  may  help  to  save 
-vch  computing  time. 

The  d-c  electric  field  at  any  point  (y,z)  may  be  derived  by  taking 
termrwlse  differentiation  of  the  potential  curve.  E(y,z)  due  to  a  line 
change  at  (y^,z^)  is  then 


E(y,z) 


-  I 


z-z  -2na 

_ o _ 

(z-2na-z^)^  +  (y-y^)* 


z-fz  -2na 
o 

(z  +  Zg-2na)®  +  {y-yo)‘* 


(B.l8^ 
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For  n  e^pproachee  ±  « 


E  (y  -)  =  I  r  - i - i _ +  r  —i _ i— 1  1 

n'-'"  2«6jj  I  I  z-z^-2na  z+Zg-anaJ  |_  z-z^-tfina  z+z^+anaj  J 


z-z 


Z-l-Z 


^  (2-z^)®  -  (2na)^  (*+“0^*  "  (2na)® 


}• 


(B.I9) 


Each  term  la  increasing  with  l/(2na)2j  as  a  result,  E(y,z)  is  a 
unlfonnly  convergent  series,  validating  the  term-wise  differentiation 
procedure.  It  also  allows  the  replacement  of  part  of  the  summation  by 
an  integration  as  an  approximation, treating  E^(y,z)  as  a  continuous  func¬ 
tion  of  n. 


E(y,z) 


r _ =^-^=^0 

2««o  "I  +  (y-yg)*  (z+Zq)®  +  (y-y©)® 


+ 


(B.20) 
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E(y,2)  = 


-(z+Zq) 


E(y,z) 


P 

iite 


i_  r _ !2c 

®0  ^  + 


z+z. 


(y-yj®  ’  (z+*  )®  +  (y-Yo)' 


1  r  I  - .  _ +  _i _ i—l 

2  L  z-z  +2a  z+z^+2a  z-z^-2a  z+z^-2aj 


n=-i 


7 


n=-“ 


aaa  z+z^-2Jja 
o  o 


i 


dn 


n=*o 


n=i 


z-z  -2na  z+z^ 
o  o 


1 _ ‘ 

-2iia 
o  J 


dn 


P 

2ite 


®o  )®  +  (y-y^)®  ^  (z-®o)®  "  (2a)®] 


1  1  1  1  r  2na-(z+z^)Y 

(z+zj=  +  (y-yj"  (z+z^)®  -  (2a)®J  2?  J 


n=-« 


1  r  2na-iz^^)  -| 

2a  [  2nB-(z-z  )  J  J  ' 


2nB-(z-z  ) 
o 


n=i 


P 

2ne 


(y-y  )®  (z-z  ) 
o'  ■  ''o'  o 


=  -  (2.)'] 


(3.21) 


(z+z  )  - - 

°  L  (z+z„)® 


(y-yo)®  (*+*0^ 


®  -  (2a)®] 


.  r  (2a+z+z  )  [2a-(z-z  )]  -i  >> 


-178- 


For  Z'>’Zq<  0.2a,  8(%ie  of  the  terms  la  Eq.  B.22  may  be  dropped. 


E(y,z)  = 


z-z 


z+z. 


2«€^  )  (z-z  +  (y-y^)®  (z+O®  +  (y-yo)^ 


In 


/  Z+Z\  /  z-z  \ 

/  Z  -Z  \  /  z  +Z  \ 


2“ol 


z-z 


ZH-Z 


(z-Zq)^  +  (y-y^)' 


(z+Zq)  +  (y-yp)‘ 


1 


z+z 
_ o 

a 


z-z  n 
_ o 

a 


p,  r  2+z  z  'V 

?  _i_  J  - 2 -  .  - 2 -  +  _2l  (B  23) 

2«e^  I  (z-z^)2  +  (y-y^)^  (z+z^)^  +  (y-y^)^  a^J  ’  ^ 

B.I.l  Uniform  Convergence  of  the  Series  Representation  of  the 
Potential  Function  V.  Uniform  convergence  of  the  Infinite  series  rep¬ 
resenting  the  potential  V  may  be  shown  by  taking  the  partial  sum  of  the 

series  for  Isorge  n.  The  nth  term  V  of  the  series  Is 

n 


■179 


(B.a4) 
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3ince  Is  proportional  to  l/n*  for  a  given  e,  there  exists  an  n^  such 

that  for  >  n^,  for  all  n'  >  n  and  a  >  0 
o 


I 


<  «  > 


(B.25) 


The  series  representing  V  is  uniformly  convergent. 

B.I.2  Location  of  the  Effective  Centroids  for  the  Electric  Field 
and  the  Space -Charge  Potential  Distribution  Calculations .  In  the  vicinity 
of  the  cathode,  the  electric  field  experienced  by  a  cheorged  psurticle  at 
(y,z)  due  to  a  line  cheurge  located  at  (y^,z^)  is  approximately 


E(y,x) 


2«€_ 


{ 


z-z. 


(z-z^)  +  (y-yg)' 


z+z 

_ o _ 

(z+z  f  +  (y-y 

O  0 


p,(y  ,z  )  2z 
~  "^o  o  I 


2nc 


o 

Z®  -  2^ 


(B.26) 


for  y  =  y^.  Now  if  the  electric  line  charges  are  uniformly  distributed 
inside  a  compartment  of  size  Az,  the  weighted  centroid  for  the  electric 
field  computation  at  the  edge  of  the  compartment  will  deviate  fr<»D  the 
geometric  centroid  of  the  compeurtment  considerably  because  of  the  manner 
by  which  the  influences  of  the  electrons  are  weighted. 


N 


I 


£.1 

z 


7" 

Oj 


N 


(B.27) 


where  z  is  the  z -coordinate  of  the  welf^ted  centroid 


-I8l- 


o 

or 


0,066  for  z  =  1 


f  Az 


and  z^j  -  0.1,  0.2,  0.3 


for  z  >  z 


...  0.9 


—  Az  for  z  <  z^ 


(B.28) 


if  z  -  z  <  Az . 
o 

The  weighting  factor  for  the  potential  computations  are  also 
evaluated  according  to  the  average  effect.  From  Eq,.  B.17,  if  y  =  y^, 

o  '■  o  ^ 

Similar  nvnnerical  approximation  will  yield  an  effective  value  of  z^  for 
the  potential  distribution  calculations. 

"z  =  0,U5  Az  for  z  <  z„  , 

O  0 

T  =  0.55  Az  for  z  >  z  ,  (B.50' 

o  0 

where  z  -  z^  <  Az . 


B.II  Dimensions  of  the  Two-Dimensional  Diode 

B.II.l  Physical  Dimensions ♦  The  size  and  d-c  characteristics  c 
the  diode  employed  in  the  Monte  Carlo  method  of  noise  analysis  are 
identical  to  those  used  in  the  density  function  method. 

Cathode  anode  distance  =  0.01852  cm 


V  .  =  10  volts 

anode 


T  =  1160.6*  Kelvin 

e 
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^IMckness  of  one  layer  of  electron  ‘beam  ^  *  6»35  x  10"^  cm 
(emlBBion  current  density)  ■  -1,5  amp/cm^. 

B.II.2  Choice  of  Unit  Electron  Beam  Area.  Theoretically,  a 
sample  Eo:ea  may  he  selected  to  represent  each  layer  of  the  electron  beam. 
For  an  ergodic  random  process,  the  length  of  the  time  period  over  which 
statistical  averages  are  taken  is  inversely  proportional  to  the  sample 
area,  yet  the  computing  time  required  for  the  d-c  electric  fields  goes 
up  as  the  square  of  the  number  of  electrons  Inside  the  diode.  As  a 
result  the  sample  area  should  be  kept  to  a  minimum  so  as  to  save  tiiM 
and,  on  the  other  hand,  the  average  number  of  electrons  emitted  per  time 
Interval  per  unit  sample  area  must  be  maintained  large  enough  so  that 
Poisson's  distribution  for  electron  emission  may  be  applied  effectively. 
Sample  ax'sa  AA  is  chosen  here  to  be  the  square  of  a  beam  layer  thickness 

AA  =  (Ay)2  .  (B.51) 

For  At  =  2  X  10  the  number  of  electrons  emitted  per  time  interval 
per  unit  area  is 


n  =  7»5509  electrons.  (B.32) 

B.II.3  Choice  of  the  Sizes  of  Compe.rtmentB .  The  sizes  of  the 
compartments  Azj  are  chosen  subject  to  the  following  restrictions : 

1.  AZj  must  be  large  enough  so  that  each  ccmpartment  (l,j)  con¬ 
tains  enough  electrons  to  make  the  fluctuations  from  the  mean  small. 

2.  It  is  also  necessary  to  limit  the  sizes  of  the  compartments 

to  avoid  unnecessary  error  in  the  computation  of  the  d-c  electric  fields, 
lumping  of  all  charges  at  the  centroid  of  a  compartment  is  valid  only 
if  the  sepaz^tiois  between  the  electrons  are  small. 
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3>  AZj  should  he  small  enou{^  to  assure  that  Identical  a-c  electric 

fields  are  experienced  hy  all  the  electrons  Inside  the  same  compartment. 

4.  A  final  condition  Imposed  upon  the  choice  of  Az's  la  that 

z  >  Az  >  V  At 
min  avg 

in  the  vicinity  of  the  potential  minimum. 

As  shown  in  Fig.  B.2  there  are  a  totaLL  of  twenty -four  compartments 
in  each  layer  of  the  two-dimensional  electron  team  model  and  their  dimen¬ 
sions  are  listed  in  the  following  table. 


Table  B.l 
Compartment  Size 


_ Az _ 

1.5  X  10"'*  cm 

6.0  X  10“*  cm 

1.5  X  10"^  cm 

0.02  X  10"^  cm 


_ Range _ 

0  <  z  <  1.5  X  10"®  cm 

1.5  X  10-®  <  z  <  4.5  X  10"®  cm 

4.5  X  10"®  cm  <  z  <  1.65  x  10*^  cm 
1.65  X  10"^  cm  <  z  <  1.852  X  10"^  cm 


Number  of 
Compartments 

10 

5 

8 

1 


B.III  Velocity  of  Electrons 

Fluctuations  of  the  velocities  at  seme  preselected  z-planes  are 
recorded  throughout  the  Monte  Carlo  calculations .  During  a  time  interval 
between  t^  and  t^^  +  At,  the  velocity  of  the  /th  electron  at  the  Zj^-plane 
and  layer  y^^^  is  defined  as  ■v^(yjj^,Zj^,ti5,) . 

For  crossing  to  occur,  it  Is  necessary  to  have 

. 

^  \  ’ 

where  Z|(tj^)  is  the  position  of  the  Ith  electron  at  Usvially, 

Z|(t^)  >  Zjj.jL  If  the  compartment  sizes  are  adequately  chosen,  i.e.. 
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n(UYER)  12  3  4  5 

IanoJe 


'//////////. 


CATHODE 


FIG.  B.2  COMPARTMENTS  IN  A  TVfO-DIMENSl(»iAL  DIODE  FOR  EIECTRIC 


FIEID  CCMPUTATICHJS 
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electrons  do  not  cross  more  than  one  coopertront  hoiindary  during  \mit 
time  Interval  At.  Then^  the  crossing  time  t^-f  St  la  evaluated  by  means 
of  the  equation  of  motion  such  that 


^n  ■  ^^*1^  =  V^l^  •  Bt  -  (h  Bta 


(B.33) 


or 


8t 


hE 


d-c 


Since  ^  Is  negative,  the  lower  sign  will  have  to  be  used  In  order  to 
bring  forth  a  positive  6t  such  that  0  <  Bt  <  At.  The  velocity  Is  there¬ 
fore 


(B.55) 


by  definition. 


B.IV  General  Harmonic  Analysis 

B.IV.l  Auto-Correlatlons  and  Cross -Correlations .  For  stationary 
random  process  with  a  sample  function  x(t),  the  time  auto-correlation 
function  R  (t)  Is  defined  such  that 

XX  ' 


R  (  t) 

XX  ' 


Urn 
T-»  <«> 


_1_ 

2T 


x(t)  x(t+T)dt  , 


(B.36) 


where  x  Is  the  time  interval  between  any  two  samples  taken  for  the  randcn 

veuriable  x(t).  This  Is,  according  to  the  ergodlc  theorem,  equal  to  the 

ensemble  correlation  function  R  ,  (x)  which  is  the  expected  value  of 

*1  *2 
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«0  00 

“  J  J  Xi(t)  Xg(t+T)  fg(Xj^,tjx^t-K)(lXj^dx^  ,  (B.37) 

-00  -00 

where  fg(x^,t;Xg,t+T)  le  the  probability  density  of  finding  x  between 

X,  and  X  -idx,  at  t  and  between  x  and  x  4dx  at  t+r. 

Ill  222 

For  x(t)  assuming  only  a  discrete  set  of  values  corresponding  to  a 
discrete  set  of  t's  equally  spaced,  the  random  process  is  called  a  ran¬ 
dom  time  series  and  its  auto -correlation  can  be  rewritten  as 


R  (s) 
xx'  ' 


a 

Ii“  y  x(t)  x(t+sAt)  .  (B.58) 

n-* » 


t=-N 


Similarly,  the  cross -correlation  function  between  two  stationary  rajqdcsn 
time  series  x(t)  and  y(t)  with  a  discrete  peirameter  t  is 


R  (e) 
xy'  ' 


Si 

=  llm  y  x(t)  y(t+8At)  .  (B.39) 

N-*  00  Aj 


t=-N 


Some  Interesting  characteristics  of  the  foregoing  equations  may  be 
revealed  by  replacing  s  with  -s. 


"xx'-"' 


=  Urn 

N-»  CO 


n 

y  x(t)  x(t-8At)  , 


(B.40) 


t=-N 


V 


s  llm 

N-^  CO 


a 

Z  y(t-8At)  . 


t  =  -N 


(B.Ul) 
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Now  if  t*  X  t-8At, 


Urn 
N-*  « 


1 

SN-t-l 


N-BAt 

.1 

t'x-N-sAt 


x(t’+8At)  x(t) 


R  (s) 

XX  ' 


(B.42) 


N-sAt 

Rj^C-s)  *  llm  y  x(t*+8At)  y(t) 

t'=-N-sAt 

=  R^(8)  .  (B.UJ) 

JA 

The  auto-correlation  functions  are  therefore  even  functions  of  time 
while  the  cross -correlation  functions  are  not. 

B.IV.2  Self -Power  Spectral  Density  and  Cross -Power  Spectral 
Density  of  Ergodic  Rtuidom  Processes .  For  a  random  time  series  x(t)  such 
that 

T 

lixa  ^  [  lx(t)|2dt  <  00 

T-» « 

there  exists  a  Fourier  transform  of  the  correlation  function  R(t)  such 
that 

00 

s(f)  =  ^  J  R(t)  e'^^'^^dT  ,  (B.lt5) 

-00 

where  S(f )  a  the  power  spectral  density  function  of  the  sample  function 
x(t)  in  the  frequency  domain,  and 


(0 


=  2xf 
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Slmllax  to  this,  the  power  spectral  density  of  a  discrete  parameter 
random  process  may  he  written  as  follows: 

R(t)  =  R(t)  6{T-sAt)  .  (B.U6) 

Therefore 

S(f)  .  t 

8=-q 

Whenever  R(8)  Is  an  even  function  of  s,  the  imaginary  term  drops  out. 

.P-i 

S(f)even  °  2n  i  /  ^  a»At  +  R(0)  +  R(p)  cos  oxiAtl  (B.48) 

s=i 


R(s)  [cos  ujsAt  -  J  sin  a»At] 


(B.47) 


or  when  R(8)  is  an  odd  function  of  s  such  eis  the  cross -correlation  func¬ 
tion 


S(f) 


odd 


=  ^  ^  [RjjyCs)  Byx(s)]  I  (cos  a»At  -  J  sin  cusAt) 


8=-q 


^  [R^yCs)  +  Rj^Cs)]  cos  a»At 


s=i 


+  R^(0)  +  I  [R^y(q)  -h  Ry^(q)]  cos  cuqAtJ  -  J  ||  |  )_  [K^^{s) 

e=x 


+  2  ox^tj-  .  (B.U9) 


APPENDIX  C.  Klisnc  rC^.’ZF.  n:  iu*  isLfiCTnor?  SEAM 


By  means  of  the  transmission-line  analog^^>  the  existence  of 
two  wave  solutions  in  a  small-signal  modulated  electron  beam  may  be 
Illustrated.  They  are  designated  as  the  fast  space-charge  wave  and 
the  slow  space-charge  wave  depending  on  their  respective  phase  velocity 
compared  to  the  average  velocity  z  of  the  electrons.  The  kinetic  poten¬ 
tial  and  the  convection  current  density  i^^  may  be  expressed  in  terms 
of  the  two  wave  amplitudes,  such  that®®, 


JP  z  -JP  z  -1 

e  +  u  e  e 


(C.l) 


and 


-  i[v 


JP  .  -»  z  1  -JP  . 

^  -  u  e  ^  e 


(C.2) 


where  u^  is  the  amplitude  of  the  fast  space-charge  wave  whose  propaga¬ 
tion  constant  is  P  -  P  , 
e  p 

u  is  the  amplitude  of  the  slow  space-charge  wave  whose  propaga¬ 
tion  constant  is  8  +  P  , 
e  p 

V  *  2V (p p/Pg)p  "the  characteristic  impedance  of  the 
electron  beam, 

Pe  = 

Pp  -  “p/z,  and 

(u^  >  the  radian  plasma  frequency  corresponding  to  the 
concentratlcxi  of  charge  in  the  electron  beam. 

From  the  definition  of  the  noise  parameter  n  in  Eq.  2.31, 
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n  -  |Re  [i*V^] 

-  5  Il*x  ^ 

■  ^  -  l“.l'j  •  (0-3) 

The  right-hand  side  of  Eq.  C.3  represents  the  real  kinetic  power  density 
in  the  electron  beam.  The  two  terms  ju^|®/2W  and  |u  |®/2W  may  be 
regarded  as  the  absolute  value  of  the  kinetic  power  density  carried  by 
the  fast  space-charge  wave  and  the  slow  space-charge  wave  respectively. 
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E 

cathode 

®d-c 
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Fourier  coefficient  for  the  nth  space  harmonic. 
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dimensional  diode. 
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s  0  otherwise. 
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partment  (m,n),  volts/m,  ’ 

z -directed  space  charge  electric  field  at  the  boimdary 
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F^(y,z,z), 
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h 
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Density  function  In  a  tvo -dimensional  space. 

Density  function  in  a  tvo 'dimensional  space  for  y  <■  0. 

Density  function  vrltten  in  terms  of  excess  energy 
parameter  W. 
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D-c  component  of  F(y,z,z,t). 

Amplitude  of  the  a-c  component  of  F^(y,z,z,t). 


Minimum  noise  figure  for  a  beam  layer  at  and  a  distance 
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Minimum  noise  figure  for  the  electron  beam  at  a  distance 
I  from  the  cathode. 
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Weighting  factor  for  the  response  part  of  the  a-c  density 
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homogeneous  part  of  the  density  function. 
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A-c  component  of  the  x-dlrected  magnetic  field  at  (y,z), 
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Average  z -directed  convection  current  density  in  the 
beam,  mp/vfi. 


io(y). 


Average  z -directed  convection  current  density  at  y  or 
amp/m® . 


i 


An  integer. 


z -directed  convection  current  density,  amp/m®. 


ij^(y  jt.)  z-directed  a-c  convection  current  density  at  the  mth  layer 
®  and  nth  plane  during  the  time  Interval  between  t^  and 

tj^+At,  amp/m^. 

1  (z  ,t  )  z-directed  a-c  convection  current  density  at  the  nth  plane 
^  during  the  time  interval  between  t^  and  t^+At,  amp/m^. 


iR 


z -directed  a-c  convection  current  density  at  per 

vinit  bandwidth  associated  with  a  unit  homogeneous  or 
response  current  input  at  the  ith  velocity  class  of  the 
mth  layer,  amp/m^. 


i  )  z-dlrected  a-c  convection  current  density  at  |  per  unit 

^  bandwidth  associated  with  a  unit  homogeneous  or  response 

^  current  input  at  the  ith  velocity  class  of  the  mth  layer, 

amp/m^ . 


^  z-directed  d-c  current  density,  amp/m®. 

J,  k,  t,  m  Integers. 


m 


Mass  of  an  electron,  kg. 


N 


Number  of  electrons  in  a  compartment. 


N(y  number  of  electrons  crossing  the  nth  plane  at  the  mth 

™  layer  during  the  time  interval  between  t^^  euad  t^+At. 

N  (y  ,z  )  Average  number  of  electrons  crossing  the  nth  plane  at 

°  ™  the  mth  layer  per  unit  time  interval  At . 


n 


An  Integer. 


P(5,5,w) 

Po(w) 


Dimensionless  density  function  parameter  defined  as 
P(?>l>w)  =  [a  (Ig/e)]"!  F(y,z,W). 

D-c  component  of  the  dimensionless  density  function  param¬ 
eter  P(5,e,w,t). 


Pj^(?>5>w)  A-c  component  of  the  density  function  parameter  P(C>l>w). 
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iR 


P 

Q„(f) 

Q^Cf) 

<1 


R(t) 


s 


R(t) 

"xxf’' 


Rj») 

V“’ 


»x  X  <’> 
1  s 


Dimensionless  a-c  density  function  parameters  at  (Cj»0  asso¬ 
ciated  with  the  kth  velocity  class  Induced  by  a  vmlt  homo¬ 
geneous  or  response  current  Input  at  the  1th  velocity  class  of 
the  mth  layer.  "±"  sign  stands  for  the  forward  going  and 
backward  going  electrons  respectively. 


The  response  part  of  the  dimensionless  a-c  density  function 
parameter  defined  such  that 


Pi( 

-  aP  (w)-j  I 

*  r77T’“'‘'^J’'’VW 


-je. 

26(w-Wj^)e 


Total  ninnber  of  computer  time  Intervals  employed  in  the 
Monte  Carlo  calculations. 

The  Fourier  transform  of  D  (t),  a  spectral  density  weighting 
factor.  ° 

The  Fourier  transform  of  a  modified  spectral  density 

weighting  factor. 

The  number  of  time  Intervals  ^t  Is  defined  such  that  two 
samples  of  a  random  variable  taken  at  a  time  interval 
greater  than  qAt  apart  are  completely  uncorrelated. 

Random  numbers. 


A  correlation  function  between  samples  of  a  random  process 
taken  at  a  time  interval  t  apart. 

Average  value  of  R(t)  between  t  -  At/2  and  t  +  At/2. 

Time  auto-correlatlon  function  for  a  stationary  randan 
process  with  a  sample  function  x(t) . 

Same  as  R  (sAt). 

xx'  ' 

Time  cross -correlation  function  for  two  stationary  random 
processes  with  sample  functions  x(t)  and  y(t^  respectively. 

Auto-correlatlon  function  for  a  random  process  with 
sample  functions  x. 


S  A  noise  parameter  defined  by  Haus. 

Sj^(f)  Power  spectrsd  density  of  the  convection  current  density 

fluctuations . 


Sg(f)  Power  spectral  density  of  the  eqvilvalent  velocity  fluc¬ 

tuations  created  by  shot  effect  from  a  thezmlcmic  cathode. 
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B^(t) 

Sg(t) 

s^M 

hM 

S(b) 

S(w+Tlg) 

s(5j,6) 

S(|) 

s(z^,n) 

s 

s 

T 

T 

t 

U_ 

u 

% 

V(i,J,m,n) 


V(1,J) 


Power  spectral  density,  Fourier  transfonn  of  R(  t)  . 

Power  spectral  density,  Fourier  txansforo  of  R( t) . 

Power  spectral  density  corresponding  to  R(t);  s  Is  a 
ccnplex  number. 

Power  spectral  density  corresponding  to  a  set  of  delta 
functions  equally  spaced  in  the  time  domain. 

S(s)  Is  defined  as  S(s)  =  {l/2«i)  ^  Sp(s-u)Sg(u)du. 

A  imit  step  function,  S(w  +  n  )  ■=  0  for  w  <  t)  . 

c  c 

=  1  for  W  a  T1  . 

'c 

Noise  parsmeter  at  ( , | ) . 

Noise  parsmeter  at  |. 

Noise  parameter  at  (y  ,z  )  and  frequencies  f  »  n/2^t. 

ID  u 

Noise  parameter  at  and  frequencies  f  =  n/2qAt. 

Complex  frequency,  cps. 

An  Integer. 

Ambient  temperature  In  degrees  K. 

A  dummy  veurlable  employed  during  the  evaluation  of  certsdn 
Improper  lntegred.s. 

Cathode  temperature  In  degrees  K. 

Time  measured  from  some  reference  point,  sec. 

At  the  end  of  the  1th  time  Interval. 

The  amplitude  of  the  fast  space -charge  wave. 

The  amplitude  of  the  slow  space -charge  wave. 

Velocity  of  electrons  In  the  beam,  m/sec. 

Average  velocity  of  the  electrons  In  the  beam,  m/sec. 

Potential  depression  at  the  boundsuT'  between  compartment 
(l,J)  and  compartment  (i+l,J)  due  to  the  presence  of 
space  charge  p  at  the  centroid  of  compartment  (m,n), 
volts. 

Potential  depression  at  the  boundary  between  compartment 
(1,J)  and  compartment  (1+1,J)  due  to  space  charge,  volts. 
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V(y,z')  Potential  at  (y,z)  relative  to  the  cathode  potential^ volts. 


V'(y) 


anode 


'd-c 


a-c 


iR 


Potential  at  the  potential  minimum  of  the  beam  layer 
centered  at  y,  volts. 

D-c  potential  at  the  anode  (applied  voltage). 

D-c  potential,  volts. 

A-c  voltage  across  the  diode,  volts. 

Chu's  kinetic  potential,  volts. 

z -directed  kinetic  potential  at  (C.,|)  per  unit  bandwidth 
associated  with  a  unit  homogeneous'^or  response  current 
Input  at  the  velocity  class  of  the  mth  layer,  volts. 


)  z -directed  kinetic  potential  at  |  per  unit  bandwidth 
associated  with  a  imlt  hcmogeneous  or  response  current 
Input  at  the  w^  velocity  class  of  the  mth  layer,  volts, 

v^  Initial  velocity  of  electrons,  m/sec, 

v^Cy  ,z  )  Mean  square  velocity  of  electrons  at  (y  ,z  ). 
o  m  n  m  n 


‘lc> 


W 


Initial  velocity  of  the  electrons  which  possess  barely 
enough  energy  to  reach  the  potential  minimum,  m/sec. 

Amplitude  of  a-c  velocity,  m/sec. 

Initial  velocity  of  the  1th  electron,  m/sec. 

Velocity  of  the  Ith  electron  crossing  the  nth  plane  at 
the  mth  layer  at  t^^,  m/sec. 


A 


Excess  energy  parameter  defined  as  W  =  z 


z2  - 


2h[V(y,z)-V'(y)  ]. 


W 

W 


Electron  beam  characteristic  Impedance  defined  as 
(aV^/l^)(Pp/Pg),  ohms. 

Average  excess  energy  parameter. 


W 


Dimensionless  excess  energy  parameter  corresponding  to  W 

w  W, 


Wj^  Dimensionless  excess  energy  parameter  of  the  kth  velocity 

class . 

X(ti,w)  A  switch  defined  such  that  X(ti,w.  )  -  1  If  w  <  0,  w+t)>0.02 

B  0  otherwise. 


X 


Real  space  coordinate,  meters. 
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X 

Y(ti,w) 


y 

y 

z 

Z 

Z 

•  • 

z 

Z^i\) 

2' 

z' 

a 

P 

P' 

Pp 


6(w-w^) 


ot 

At 

AA 

Az 


Time  derivative  of  x;  x-directed  velocity,  m/sec, 

A  switch  defined  such  that 

Y(ti,w)  =  1  lfw>0  or  w<0  auid  w+tj  =  0.02 
=  0  otherwise . 

Real  space  coordinate,  meters, 
y-dlrected  velocity,  m/sec. 

R-f  Impedance  of  the  diode,  ohms. 

Real  space  coordinate,  meters, 
z-directed  velocity,  m/sec. 
z-dlrected  acceleration,  m/sec. 

z-directed  space  coordinate  of  the  ith  electron  at  t^ 
meters . 

Space  coordinate,  a  dummy  variable,  meters. 

z-directed  velocity  of  electrons  at  a  layer  y  =  y',  m/sec. 

A  constant  defined  as  a  =  m/21tT  . 

c 

A  constant  defined  as  p  ^  2n^^*  (e|l  I/e 

0  0 

Phase  constant  defined  as  p'  ^  co/z,  m 

Phase  constant  for  a  single  velocity  electron  beam, 

Pg  =  <o/z,  m  . 

Plasma  phase  constant  defined  as,  m  Pp  = 

Delta  function  defined  as  b(j,m)  ^1  if  J  =  m 

=  0  if  J  /  m. 

Delta  function  defined  as 

Wi+€ 

[  6(w-w^)  dw  =  1,  6(w-Wj^)  =  0  for  w  /  w^^ 

=  »  for  w  =  Wj^. 

Incremental  change  in  time. 

-12 

Unit  time  interval  At  =  2  x  10  sec. 

Sample  emlsslcxi  area  ofsnelectz^  beam  layer, 

Incremental  change  in  z,  meters. 
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s 

m  u 

9(y,z,s) 


Aityl) 

a(|) 

V 

I 

i 

t  • 

I 

n(y  ,z  ,n) 
•'m  n'  ' 

n( 5j,l) 

n(i) 


Incremental  change  in  phase  angle. 

Dielectric  constant  of  the  free  space. 

Dimensionless  real  space  coordinate  parameter. 

Time  derivative  of  5* 

Potential  parameter. 

Potential  parameter  at  the  cathode. 

Potential  parameter  defined  as  0  =  2oaV'/m. 

Cross -power  spectral  density  of  the  convection  current 

density  fluctuations  and  the  kinetic  potential  per  unit 

bandwidth  at  (y  >z  )  and  frequency  f  =  D/2(^t. 
m  n 

Cross -correlation  function  between  the  convection  current 
density  fluctuations  and  the  kinetic  potential  at  (y,z). 

Phase  shift  of  the  original  disturbance  associated  with 
the  ith  velocity  class. 

6 

+  A  f  d£  __ 

9.  s  ±a  /  -  -  a  —  .  The  negative  sign  is  for  the 

backward  going  electrons. 

Imaginary  part  of  6(y^,z^,fi) . 

Imaginary  part  of  the  cross -power  spectral  density  of  the 
current  fluctuations  and  the  kinetic  potential  at 

Imaginary  part  of  the  cross -power  spectral  density  of  the 
current  fluctuatiozis  and  the  kinetic  potential  at  |. 

Frequency,  dummy  variable. 

Dimensionless  real  space  coordinate  parameter. 

Time  derivative  of  | . 

Second  time  derivative  of  |. 

Real  peurt  of  8(y  ,z  . 

n  u 

Real  part  of  the  cross -power  spectral  density  of  the 
convection  current  density  fluctuations  and  tge  kinetic 
potential  at 

Real  part  of  the  cross -power  spec^^nO. 
convection  current  densl^ 

potential  at  J.  " 
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Pl 

II 

Pi 

Pin,n 

[pv](y,z) 

T 

Hiyi) 

♦(5) 


♦(6) 
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Space -charge  denelty  In  regions  I  and  II  respectively 
coul./a*. 

Space  charge  per  unit  length  in  the  x-dlrectlon,  coal./afi. 

Space  charge  per  unit  length  in  the  x-dlrectlon  at  the 
centroid  of  ccopartoent  (n^n),  coul./m. 

z -directed  convection  current  density  In  coopartment 
(y,z),  an®/o®. 

Time  measured  from  seme  reference  point,  sec. 

Self -power  spectral  density  of  the  kinetic  potential 
at  <uid  frequency  f  =  0/2gdt. 

Auto-correlatlon  function  of  the  kinetic  potential  at 
(y>z). 

Self -power  spectreLl  density  of  the  kinetic  potential  at 

Self-ixjwer  spectral  denelty  of  the  kinetic  potential 
at 

Self -power  spectral  density  of  the  convection  current 

density  fluctuations  at  (y  ,z  )  and  frequency  f  *  0/2cAt. 

n  n 

Auto-conelation  of  the  convection  current  density 

fluctuations  at  (y  ,z  ). 

'm'  n 

Self -power  spectral  denelty  of  the  convection  current 
density  at  ( 5j,6) . 

Self -power  spectral  density  of  the  convection  current 
density  at  | . 

Frequency  parameter  defined  as  n  =  n/2q^t,  sec."^ 

Radian  frequency  defined  as  cu  =  2nf. 
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